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Abacca it 


Measurements have been made of the heat transfer rates through tur- 
bulent and transitional boundary layers on an isothermal, convexly 
curved wall and downstream flat plate. The effect of convex curvature 
on the fully turbulent boundary layer was a reduction of the loca? Stan- 
ton numbers 20-50% below those predicted for a flat wall under the same 
circumstances. The recovery of the heat transfer rates on the down- 
stream flat wall was extremely slow. After 60 cm of recovery length, 
the Stanton number was still typically 15-20% below the flat-wall pre- 
dicted value. 


Various effects important in the modeling of curved flows were 
studied separately. These are: (1) the effect of initial boundary layer 
thickness, (2) the effect of freestream velocity, (3) the effect of 
freestream acceleration, (4) the effect of unheated starting length, and 
(5) the effect of the maturity of the br-undary layer. Regardless of the 
initial state, curvature eventually forced the boundary layer into an 


asymptotic curved condition characterized by St a Re 


„ -1 


The slope, 


minus one, is believed to be significant. The only exception to this is 

the. case of strong acceleration with curvature where the exponent on 

Re* was found to be -2. This study showed that variations of the 

2 

initial parameters affected only the rate at which the asymptotic state 
was approached. Strong convex curvature is an organizer of the heat 
transfer data, so long as the boundary layer contains sufficient turbu- 
lence. The data on the flat walls upstream of the curved section varied 
considerably from case to case, but within the curve the data were simi- 
lar. The effect of curvature on the mean temperature profile was a 

strengthening of the wake region; there appears to be no lasting effect 
on the structure of the inner portion of the boundary layer. Freestream 
acceleration, on the other hand, has been shown to affect the inner re- 
gion. The present data show that acceleration and convex curvature are 
additive in their effects on Stanton number. 


The effect of convex curvature on laminar and early transitional 
boundary layers, where the turbulent mixing was minimal, was to delay 
and retard the transition process. When curvature was removed, transi- 
tion appeared to progress more rapidly. 
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The physical Interpretation given to these observations is that 

convex curvature compresses the scale of the turbulent boundary layer. 

Tills yields an inner region of nearly constant thickness and mass flux 

with vigorous turbulent mixing activity. Beyond this is an outer layer, 

scaled on the boundary layer thickness, which contains decaying rosidual 

turbulence from the boundary layer that existed upstream of the curve. 

A bent transfer model ean be devised, consistent with the above, where 

the overall heat transfer rate is conduction-limited. This would result 
"1 

in the Re, relationship of the asymptotic curved boundary layer. 

When curvature is removed, the inner layer begins to grow and consumes 
the outer residue layer* 

An existing curvature-prediction model was tested against this 
broad heat transfer data base to determine whore it could appropriately 
be used for heat transfer predictions. 
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Chapter 1 
INTRODUCTION 


Turbulent boundary layers on convexly curved walls ore encountered 
in many engineering applications; tike iorward part of a blunt body, 
leading edges oi air intakes, blade passages of turbomachinery , aircraft 
wings, and rocket nozzles. Curved boundary layers with high heat trans- 
fer rates are encountered on gas turbine blades, where accurate predic- 
tion oi and design ior turbine blade heat loads are critical to the 
reliability and efficiency oi modern high- per lorinance engines. Hi ere is 
ample evidence that curvature affects heat transfer. Some of this ef- 
fect lias been attributed to the extra rates ot strain associated with 
streamwise curvature which significantly a fleet the structure oi turbu- 
lent boundary layers [2,33] and the heat transfer rates [13,32). 

Hie primary objective ot this experiment was to measure the effect 
of convex curvature on the heat transfer rate over a representative 
domain oi initial and boundary conditions. Hie work was undertaken as 
part oi an ongoing series ot projects at Stanford University sponsored 
by NASA-Levis Labs. The motivating problem was the need to understand 
the mechanisms and accurately predict the heat transfer rates on gas 
turbine blades. 

Wien the curvature project began, the state, of the art tor design 
heat-load calculations was the computer code STAN* [1]. This code 
solves the partial differential equations which govern transport of 
thermal energy and momentum in boundary layers. Hie quality of the pre- 
dictions it makes Is dependent upon the applicability oi the Reynolds 
stress model It, uses. Though the stress modeling of STAN 5 is quite 
simple, it is a trusted program for use within Its data base, because 
its Reynolds stress model is supported with empirical input from twelve 
years oi careful experimentation at Stanford. STANS accurately accounts 
for streamwise acceleration or deceleration and/or transpiration blowing 
or suction. Presently, however, it does not account for streamwise cur- 
vature. Hie purpose of the Stanford curvature program is to add curva- 
ture to its useful domain. 
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I .1 Previous Research 

A comprehensive survey oi the literature on curvature effects, 
prior to 1972, was given by Bradshaw [2] * Ikj pointed out that the 
effect of curvature is about ten times as strong as one would predict 
from a thin shear layer, eddy-viscosity model by simply adding the extra 
rate of strain OV/Ox to the existing strain Held. 

Experimental wor 1 on the effects oi curvature dates back to the 
days of Ludwig Pra One early study by Wilcken [3], a student of 
Prandtl, was the f documented study oi curvature effects where the 
facility was designed to keep secondary ilow effects to a minimum. At 
that time, Prandtl had put forth a stability argument for curvature, 
Wilcken was to test it. 

"Because the various parts oi the boundary layer are 
variously affected by centrifugal force in the presence of a 
curved surface, a concave surface produces a tendency to 
force the fast parts of the flow toward the surface and the 
slow parts away from it. This tendency lavors the exchange 
of the alow layers next to the surface with the faster ones 
on the inside of the flow. IVius , it reinforces the already 
existing turbulent exchange procedure, The contrary is the 
cose for the convex laces. Here the centrifugal force has a 
stabilizing effect, reducing the turbulence. 

Prandtl also had developed the early form of a mixing length model, then 
called "free path," tor a flat plate, Wilcken, though admitting some 
secondary ilow influence, found that curvature significantly affected 
the "free, path lengths," Ifc stated, "Boundary layer events on curved 
surfaces should be ascribed more importance than has generally been the 
case up to the present «" Curved flow research continued under Prandtl at 
the Kaiser Wilhelm Institute for Flow Research with the studies of Vtendt 
[4] and Schmidbauer ['>] . The general conclusion was that, even for weak 
curvature, the boundary layer hydrodynamics are significantly affected. 

The next documented study was in a curved channel, where the flow 
could become fully developed. This study by Wattendorf [b] in 1934 
showed that the full}. developed flow was significantly influenced by 
curvature, although the overall pressure drop was not. Wattendorf also 
found, from his mean velocity and wall static pressures measurements, 
that the power, n, in the powi»r-law velocity profile, U + = cy^^ n , 
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decreases with stronger convex curvature and Increases tor stronger 
concave curvature. Ills descriptor of the strength oi curvature was the 
parameter v/Ru , where R was positive for concave curvature. 

In 1937 , Clauser and Clauser [7] investigated the effect oi 
curvature on the transition from laminar to turbulent ilow. This ia 
the ttrst time that hot-wire anemometry had been applied to curved 
flows. They found that stabilizing convex curvature increased Lhe 
critical Reynolds number and delayed transition. Their observations 
showed that, on typical airioils of the time, convex curvature might 
double the critical Reynolds number. Hans Ltepmann [8] extended the 
study oi the effect of curvature on transition to very weak curvature 
CO < 6., \ 0.001). In this range it was found that transition was not 

La 

affected by convex curvature but was affected by concave curvature. 
Tills indicated that the process of transition was different for the 
concave wall than lor the flat or convex wall.. 

in 1955, Frank Kreitli [9] performed a clever heat transfer test 
that showed quite conclusively that the heat transport from a concave 
wall was considerably more than from a convex wall. For his channel 
ilow, he concluded that curvature effects scaled on U/r, the forced 
vortex parameter. No measurements of local heat transfer rates were 
taken. Hie hydrodynamics of a fully developed, curved, turbulent chan- 
nel flow similar to the Krelth facility were studied in detail by 
Eskinazi and Yell [10]. Using hot-wire anemometry, they measured the 
downstream development of profiles of streamwise velocity fluctuations 
and, for the fully developed flow, measured the profiles of u' , v' , 
and u T v T • They stated that one of the most important iniluencee of 
curvature Is on the v' -production term -u'v’ (U/r), They noted that 

near the convex wall u’ T vJ was positive, indicating a suppression 

* ‘ 2 2 

(negative production) oi v' . Their spectral measurements of u' and 

v' 2 showed that the decrease in turbulent mixing activity was largest 
in the low-wave-number range. 

The first curved-flow heat transfer test with wall-measured heat 
flux data was that of Schneider and Wbde [LI] . With plug-type heat flow 
transducers they measured local heat fluxes on a convex wall that were 
50% of the predicted flat-wall values and considerably less than would 
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be predicted by the model ol Krelth [9j. llielr tunnel aspect ratio 
was 1.0, so considerable. contamination by secondary flow was prob- 
able . 

V. C. Patel [12] studied the hydrodynamics ot the flow through a 
90 p curved duct similar to the facility used in the Schneider and Wide 
study, except that the aspect ratio was '>.(). ho attempt was made to 

separate streamwise acceleration and deceleration from the other effects 
oi streamwise curvature, and there was some contusion about whether the 
change in the mean velocity protiles was due to curvature or the local 
acceleration and deceleration within the bend. Patel measured only mean 
quantities, which, even with the increased aspect ratio, may have been 
influenced by secondary flows. lie noLed that curvature affects the 

shape factor and, therefore, the rate of entrainment. 

Also In 19t)8, Thoraann [13], at the Aeronautical, Research Institute 
of Sweden, made, detailed local heat transfer measurements on surfaces 

that were straight and eonvexly and eoncavely curved. The. Ihomann 

study was performed in a wind tunnel with a treestream Mach number of 
2.5, uniform static pressure on the tost wall, and a boundary layer 
thickness-to-radius of curvature, ratio of 6^/R n U.02. The resultant 
effect of curvature was found to be an increase in convective heat flux 
of ^ 20% for the concave case and a decrease of « 15% for the convex 
case with respect to the flat-wall case (see Fig, 1-1). Since the 
Thomann study was done in a supersonic freestream, significant compres- 
sible effects were present. As discussed by Bradshaw [2], compression 

or dilation produces strong extra rate-of-strain effects in the boundary 
layer which may alter the turbulent transport process in much the same 
fashion as does curvature. 

In I9b9, Bradshaw [14] discussed the analogy between streamwise 
curvature and buoyancy in turbulent shear flow and introduced a modifi- 
cation of the Richardson number used In meteorological work to curved 
and rotating flow computation. The gradient Richardson number for 

curved flows is written; * 2S(1+S) where S is the stability 

parameter, S * (U/R)/(3U/3n) , positive for the convex wall and nega- 
tive for the, concave wall. He then proposed [14,15] that the ftonin- 

Obouhkov formula for the correlation of the apparent mixing length with 


4 



small buoyancy effects H/SL q - 1 - 8Ri could be used to model the 
effects of weak curvature. Ibis approach met with considerable success. 
In fact,, the value of the constant 6 could be inferred by analogy from 
meteorological experiments in stably and unstably stratified boundary 
layers. It is generally agreed to be the order 10. 

So and Mel lor [16,17,18] published results from a very detailed hy- 
drodynamic experiment on curved-wall boundary layers . In their experi- 
ment, the ratio of boundary layer thickness to radius of curvature ratio 
was 6^/R * 0.07, and the aspect ratio was ~ 8.0. Because of an 
imaginative design which employed wall jots, secondary flows were kept 
acceptably small. Profiles of all the Reynolds stresses were measured. 
On the convex wall it was found that the turbulent shear stress was 
"turned off" in the outer half of the boundary layer. Over the concave 
wall, they found evidence of a stationary system of longitudinal vor- 
tices, analogous to those formed between rotating cylinders. Wall shear 
stress was inferred from a Clauser plot technique, but the turbulent 
shear stress profile was not measured close enough to the wall to enable 
them to check the wall value by extrapolation. Ihe Clauser plot tech- 
nique has since become an accepted method for inferring wall shear 
stress and has been verified by later experimenters, e .g . , Gill-' [35] • 

Kilts and Joubert [19] of the University of Melbourne measured 
profiles of mean velocity in boundary layers of a curved duct. they 
noted that the width of the logarithmic zone was curvature-dependent. 
Convex curvature caused the velocity profile to become wake-like at a 
lower value of n and the opposite for concave. They searched for, 
but could not find, similarity of the mean flow, as represented by a 
defect type law for either the fully developed, curved, turbulent 
channel flow or the outer regions of convex wall boundary layers . On 
the concave surface, Ellis and Joubert found evidence of Taylor- 
Gortler vortices. 

IXjring the 1970s, Bradshaw undertook a series of experiments on 
curvature effects. The first experiment was one of very weak curvature 
( ^ 9 95 early results were presented by Meroney [20] and 
the final results by Hoffmann [21] . Mean and turbulence data were taken 
in a curved duct of 30° bend and aspect ratio of 6.0. No attempt was 
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made to separate the weak, viecolcrution/deeulerat ion et tecta . The data 
constated ot mean quantities and the important higher-order quantities 
through order 4. Spanwise variations of skin irictlon due to weak 
Taylor-Gortier cellular activity were observed on the concave wall, they 
measured variations ot 20% in local skin-irict ion coefficient. Utet 
experiments on curvature at the Imperial College included the recovery 
from an "impulse" of curvature in a curved duct (Sraits, Young, and 
Bradshaw [22]) and in an axisyminetric flare (,SmiLs, Eaton, and Bradshaw 
[ 2 J ] ) , 'lhe second case showed the combined effects of sustained lateral 
divegenee and recovery from curvature. It was found that Taylor- 
Gortier cells were established in the curved duct but not in the 
axisymmetrie flare, indicating an effect of lateral divergence. In both 
studies, mean data and fluctuating data through order 4 were taken. 
During this same time, Castro and Bradshaw [24] were studying the eon- 
vexly curved, free-mixing layer. 

Effects oi mild curvature were studied by Rauaprian and Shlvapra- 
shad [23,26,27,28]. Their test facility was very similar to the Hoii- 
mann and Bradshaw facility, with the noLable exception of an aspect 
ratio of 2.3. They did not observe Taylor-Gortier cells on the concave 
wall, contrary to the findings of Hoffmann and Bradshaw. The authors 
admit that this difference may be attributed to significant secondary 
flows . 

Detailed measurements of the hydrodynamics of a curved channel flow 
were made by Hunt and Joubert [29] in 197 9 . The data for this study 
were taken in an apparatus in which a duct downstream of the nozzle 
could be straight or bent into a large-radius curve. The channel was 
shallow, allowing the convex and concave boundary layers to merge about 
40% of the distance, around the 45° bend. The curve started at the exit 
of the nozzle, so the. ratio of boundary layer thickness to radius of 
curvature was extremely small < ^ 99 “ 0,0 °^ aC the 8tart of curva- 
ture, then grew rapidly around the bend. The aspect ratio was large, 
13.2 to 1.0, to minimize secondary flow effects. Roll-cell vortices of 
the Taylor-Gortier type were observed in the channel just off the con- 
cave wall. These cells resulted in variations in skin friction across 
the span of ± 4%. Data taken included mean velocity profiles and 
profiles of the important Reynolds stresses. It was found that the 
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mensured Reynolds shear stress extrapolated to the convex wall agreed 
with skin iriction values found with the Ciauser plot technique. 
Profiles at the end ol the test region were iully developed in the mean 
quantities and nearly iully developed in the turbulence quantities* 

In 197/, Hr ini oh and Graham (30] measured mean velocity and 
temperature protiles and wall heat transier rates in a curved channel 
similar to the one in which Eskinazi and Yeh [10] took hydrodynamic 
data. The heat transier facility, however, had an aspect ratio of 0.0, 
and considerable secondary flow influence was observed. Nevertheless, 
they noted an increase in the heat transfer rate on the outer (.concave) 
wall and a slight decrease on the inner (convex) wall. Their tempera- 
ture profiles reached an asymptotic shape characterized by a skew ol the 
peak toward the concave wall, giving a much steeper temperature gradient 
on the concave than on the convex wall. Temperature profiles were not 
plotted in inner coordinates, so it Is not known whether they showed a 
log- Law relationship. 

R. M. 0. So [31] continued his curvature work beyond the hydrody- 
namic studies discussed earlier [lb, 17, 18] with an analytical prediction 
ol the eliect ol curvature on the turbulent Prandtl number. He devel- 
oped the rather simple relationship: 

1 - } 8 ' Ri(l - 

l M _ _ ay" 

‘h i -i 8 '(/Mu - i) ' +i v ) 

where y is the value ol the turbulent Prandtl number of a correspond- 
ing plane flow. The value 8' “ b .0 has been found to give the best 
correlations with curved-flow data, swirling-flow data, and meteorologi- 
cal data [85,bb,57]. This equation gives an increase of ior 

> 0 (convex) and Pr <. 1 . Ibis ratio will decrease if, either one of the 
two inequalities is reversed. For air flowing over a convex wall, the 
diffusivity of thermal energy decreases faster than the eddy diffusiv- 
ity. Since the Prandtl number for air is nearly unity, the deviation 
from Reynolds analogy ( e >i/ e h ” 1) i 8 predicted by the above equation 
to be small. 
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Recently, Muyle, Blair, and Kopper [32] measured local heat trans- 
fer rates In a curved duct which had an aspect ratio oi 4.2 5* Their 
study was with low-velocity air (M * 0.06 as opposed to the high- 
velocity study oi Thomarm [13]), in which the test wall static pressure 
and temperature were uniform. The ratio of boundary layer thickness to 
radius of curvature at the start oi the curve was ~ 0.01. Their find- 
ings (see Fig. 1-2), confirming those of Thomann, showed a dectease in 
the wall heat flux of * 20 % lor the convex wail and an increase oi 

« 33 % for the concave wall relative to the flat-wall heat flux predic- 
tion of Reynolds, Kays, and Kline [33]. 


Previous heat transfer studies dealt oniy with the curved region, 
none reported results of recovery downstream of the curve. The present 
investigation is the first time the effects on heat transfer of both the 
introduction of and recovery from curvature have been investigate*;, 
when looking at the problem of gas Lurbine blade cooling or many other 
applications where curvature effects on heat transfer are significant, 

f t Ku » 1 n »*niv 'atm » v * o f ^ Anrv f p r*ii At'f’pp i Al 1 

UUL* I tut: I. Vjt, x.vuq w* u v,iwu ft t " '"*• h- j 

regions of weak or no curvature, so the recovery effects are as perti- 
nent as the effects of the introduction of curvature. 


The present heat transfer results are from a program studying both 
the heat transfer and fluid mechanics of a convexly curved flow. Since 
heat transport Is by turbulent motion, the hydrodynamic study is an 
essential input to the heat transfer study. Detailed results of the 
hydrodynamic study have been reported by Gillis and Johnston [34,35]. 
Their work was performed on essentially the same tunnel configuration as 
the present study (see Fig. 2-1) and with boundary layer thlckness-to- 
radius of curvature ratios of 0.1U and 0.05, two of the cases pre- 
sented herein. The Gillis and Johnston results are summarized in Figs. 
1-3 through 1-7. Fig. 1-3 shows the effect of convex curvature on skin 
friction, as calculated using a Clauser plot. The initial response to 
the introduction of curvature was fast, and the behavior was seemingly 
near-asymptotic. At the end of the curved region (- 16 boundary layer 
thicknesses downstream from the start of curvature for the <■> 
0.10 case), the skin friction reached a value ~ 30% less than the 
flat-wall predicted value. In the recover> region, the return of the 
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■kin friction wan alow. Altar ~ 20 boundary layer thickncaaea Iron 
the atart of the recovery region, the akin friction waa atlll ** 20-2)1 
below the flat-wall value. They alao found that, though the boundary 
layer thlckneaa at the beginning of curvature differed by a factor of 
two for the two experlaenta, the variation of Cj/2 with diatance 
within the curved region waa aurprlslngly similar tor the two caaea . 
Fig. 1-4, a typical plot of mean velocity profllea In wall coordlnatea, 
shows that the log region ahortena within the curved section but la 
always d lscerniblle . The effect of convex curvature on profiles of the 

— 5 "7 

turbulence quantities (u' , q , and u 1 v’ ) la shown on Figs. 1-5, 1- 
6, and 1-7, respectively. Within the curved region, the turbulence 
activity declined quickly throughout the boundary layer, but especially 
in the outer regions, to a near aelf-slmllar profile. 'Hie recovery to a 
normal flat-wall profile was slow and appears to be by propagation from 
the wall. The slow growth of the turbulence profiles in the recovery 
region may be responsible for the slow return of the skin friction to 
flat-wall ’alues. Fig. 1-7 shows that, at the beginning of curvature, 
the shear stress in the outer portion ot the boundary layer changes sign 
for a short distance. This term is linked to the production of TKE a 
hi ;n-rcviTNdl indicates a negative production of TKE in this region. 

Clllls and Johnston identified a parameter called the "shear layer 
thickness," the n-dlstance of the extrapolated shear stress profile (see 
Fig. 1-8) labeled " ’6^ ." They proposed a model by which to visualize 
the strongly curved boundary layer: a two-layer model shown in Fig. 1-9. 
The Inner layer Is characterized by non-zero TKE and shear stress, a 
typical boundary layer, but constrained in thickness to 6^ , which 
remains essentially constant within the curve, The outer layer is char- 
acterized by non-zero TKE but essentially zero shear stress. The outer 
layer is the decaying residue of the upstream boundary layer. At the 
beginning ot the recovery region, the restriction on due to curva- 
ture is lifted and the inner layer grows within what remains of the 
outer layer. The results of the present heat transfer study are consis- 
tent with this model. 

The hydrodynamic study of Clllls and Johnston is part of an over- 
all program involving heat transfer rates on smooth, convexly curved 
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mrficti. Tlit* iri'kviit experiments build upon their study of the hydro- 
dynamic processes, Invest lgat lng the heat transfer processes through 
structurally similar, but heated, boundary layers. 

1.2 Objectives 

The final objective of the Stanford curvature program Is to add 
curvature to the list of effects that can be appropriately predicted 
with turbulent boundary layer prediction programs, e .g . , STAN5 ( l J . 
With this goal In mind, the objectives ot the present study were to: 

1. Accurately measure wall heat transfer rates from a smooth, convexly 
curved wall and dot. ist ream flat recovery plate over a large enough 
domain ot carefully controlled Initial and boundary conditions so 
that the eventual model can be used with minimal extrapolation. 

2. bit Id upon the understanding of the curved boundary layer gained 
during the dills and Johnston study (33) and others listed in 
Section 1.1. For this study a series ot 15 runs were made with 
differing Initial and boundary conditions to learn about the sen- 
sitivity of the curvature effect to various parameters thought to 
be important. This will aid in developing a prediction model based 
on maximum physical insight . 

3. Teat the preliminary prediction model proposed by Gillis and John- 
ston (35], making changes where appropriate. 

4. Construct a heat transfer facility that Is sufficiently flexible 
to accomplish the above objectives and future objectives of the 
overall program, specifically, convex curvature with discrete Jet 
injection, which simulates modern gas turbine f /.lm-coollng geom- 
etries . 


1 .3 The Experim ent. 

in the following experiment, thermal and hydrodynamic boundary 

layers were grown on a flat preplate, then were introduced to a 90®, 45 

cm radius of curvature convex wall, followed by a recovery wall. The 

wall temperature was maintained uniform, and the static pressure or. the 

test wall was controlled to a prescribed function of downstram distance, 

2 

either uniform or constant K (v/ll • dU /ds) . Detailed wall heat 

pw pw 
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flux data aa wall aa profllea of velocity and temperature ware mea- 
sured In the developing, curved, and recovery regions. The Initial and 
boundary conditions were varied to determine sensitivity to various 
parameters. The 13 cases In the present study Investigated the sensi- 
tivity of the effect of convex curvature on heat transfer t?: 
a) Initial boundary layer thickness, 
b; Mag nltude of free-atream velocity. 

c) Free-stream acceleration. 

d) Location of the beginning of heating with respect to the beginning 
of curvature or the beginning of recovery. 

e) Maturity of the momentum boundary layer at the beginning of curva- 
ture . 

One other parameter which should be varied syatemat leal ly la the 
radius of curvature of the test wall. Itie cost and complexity of a test 
facility with such flexibility precluded this entry. 


11 



STANTON NUMBER, St(T*/T,J 


•7 - 



■ ■ i 


1 1 

o 

6 

•5 - 

* 

" » m 

' w 

— ji „ 

* „ A 



'I 

» ■ "O 

1 a 

0 

k 

» 

I 

2 *> <3 

* — — <4 

au. , 

\ 

> ^ 

* A 

•1 



0-9 1 

"o H l 

• 

2 1 

u 

1 1 


^ / 


Fig. 1-1. Stanton number versus streumwlse distance, M *• 2 . 
T l( -T w - 77.5 k , froa Thn— nn [13) • 



REYNOLDS NUMBER, Re* 

Fig. 1-2. Stanton number versus x-Reynolds number, from 
Mayle, Blair, and K opper [32]. 
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Skin friction versus streamwise distance, Sgg/R 
0.10, from Glllls and Johnston [35]. 
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Fig. 1-4. Moan velocity profiles In wall coordinates, 5q t) /R 
0.10, from Gill is and Johnston [35], 
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Fig. 1-5. Streamwlse turbulence Intensity profiles versus 

distance In the streamwise direction, >S 99 /r - 0.10, 
from Glllis and .lolmston [35]. 
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Fig. 1-7. Reynolds shear stress profiles versus distance in 
the streamwise direction, ^q/R ■ 0.10, from 
Cl 11 is and Johnston [35] . 
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Fig. 1-8. Reynolds shear stress profile showing " r s i". 1 rom 
Gill is and Johnston [15]. 



Fig. 1-S . Two-layer curved boundary layer structure, from 
Gillis and Johnston [35]. 
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Chapter 2 

THE EXPERIMENTAL APPARATUS 


The test facility used In the present heat transfer study was used 
in the hydrodynamic study of Glllls and Johnston [35] • The hydrody- 
namic aspects of the tunnel, discussed In detail In Ref. 35, are siasmar- 
lzed below. Systems added to the facility to broaden the hydrodynaaic 
operating domain and systems specific to the heat transfer measurements 
are discussed in detail. 

2 . 1 General Description 

A curved-wall heat transfer facility has been constructed (Fig. 
2-i) that allows development of a smooth, flat-wall, heated boundary 
layer upstream of a 90° bend of 45 cm radius of curvature followed by a 
flat-wall recovery. There were five circuits to the facility: 

(1) The main loop: from the main fan, through the return 

ducting, oblique header, heat exchanger and screen pack 
and contraction nozzle combination, then into the test 
regions and back to the main fan via a plenum box. 

(2) The charging loop: discharge flow to the room via louvres 
and slots then return to the tunnel via the filter box and 
charging blower. 

(3) The suction loop: suction from the preplate, reinjected to 
the plenixn box via the suction fan. 

(4) The cooling water loop for the heat exchanger: a 3031 

capacity water tank, water supply, and discharge lines to 
and from the tank, for circulation and for make-up. 

(5) The hot water loop: heated the preplate and recovery 

walls, using two temperature-controlled water heaters. 


*Cases 2 and 3 of Ref. 35. 
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2.2 Hydrod ynam ic Aip«cti of the Facility 

The main flow waa driven by a fan which delivered air Co an oblique 
header Chat turned Che flow Into a heat exchanger. The flow paeaed 
through Che exchanger, a screen pack, and an 11:1 contraction nozzle 
before entering the tunnel teat region. Oetalla of the acreen pack and 
nozzle can be found In Ref. 44. Flow exited the teat region Into a 
plenum box, which aupplled the main fan. Hie tunnel velocity waa 
controlled by changing pulleya and belts on the fan and motor, and could 
be varied from 3.5 m/a to 26 m/ s . 

The developing region waa 16.5 cm by 56 cm In cross section and ~ 
200 cm long. The outer wall was straight and adjustable by pivoting 
about Its upstream edge ao that the static pressure could be made nearly 
uniform, allowing the growth of a normal flat-wall (non-accelerated) 
boundary layer. A 15 cm section of the test wall In the preplate 
beginning ~ 84 cm upstream of the start of curvature was constructed 
with ~ 2000-1/16 inch diameter holes uniformly spaced across the span 
and connected to an auxiliary fan (Fig. 2-la) allowing boundary layer 
auction. Tills system extended the operating domain to thinner boundary 
layers and very low momentixo thickness Reynolds number laminar and 
transitional boundary layers. When a thick, fully turbulent, two- 
dimensional boundary layer was desired at the beginning of curvature, 
the boundary layer was tripped just downstream of the nozzle. When the 
suction fan was operating to produce thinner boundary layers, the trip 
was located downstream of the suction holes. 

Within the curved region, the flexible outer (concave) wall was 
adjusted so that the static pressure on the convex test wall followed 
the desired function of streamwise distance, usually uniform. When the 
static pressure on the test wall was uniform, there was no streamwise 
acceleration of the Inner region of the boundary layer. Cases with 

acceleration were set up by trial and error until a nearly constant K 
was achieved within the curve. Twenty-six 0.025-inch (0.6 mm) diameter 

static pressure tap holes distributed in both the streamwise and span- 
wise directions were used for these adjustments. Hie tunnel was main- 
tained slightly above ambient pressure with the charging blower (Fig. 
2-la), which took air from the room via a filter box. Because the 
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tunnel was pressur 1 re J , separation of the boundary layer on the concave 
wall, upon entering t he curved region, could be prevented by discharge 
o! the boundary layer fluid through a series of seven apunwlae louvres. 
Secondary flows in and downstream of the bend were reduced by peeling 
off the sidewall boundary levers rpon entry to the bend, by using side- 
wall slots within the bend, and by installing boundary layer fences on 
the convex surface near the side walls beyond the heated portion of the 
span. The secondary flow control was developed by J. C. Gil 11s. The 
details of the evolution and a more thorough discussion of the final 
design are presented In Ref. 35. 

The recovery region was a straight tunnel of dimensions 15 cm by 53 
cm and approximately 125 cm long. Boundary layer fences continued down- 
stream of the curve for the first 60 cm of recovery length. 

The outer walls were constructed so that profile data could be 
taken at five stations In the developing region, six stations within the 
curve, and six stations within the recovery region. Stations typically 
had seven spanwise positions for checking two-dimensional ity . 

2.3 Heat Tra n sfer Aspects of the Fa c ilit y 

Tiie test wall was constructed of copper and segmented and Instru- 
mented so that the local heat flux could be measured on the preplate, 
the convex ly curved wall and the recovery wall. 

The preplate was divided Into 48 segments, each 2.5 cm long in 
the streamwlse direction and each Instrumented with an embedded lron- 
constantan thermocouple. The last 24 segments were heated with circu- 
lating hot water (Fig. 2-lb) and were Instrumented tor direct heat flux 
measurement. Each segment consisted of a 3 mm thick copper plate backed 
with 3-0.5 mm thick bakelite sheets, the center one containing a silver- 
constantan thermopile at the centerline location. The thermopile signal 
was correlated with heat flux by calibration (see Appendix C for cali- 
bration details). Ihe heated water circulated through a 1.2 * 2.5 cm 
copper waveguide channel behind the bakelite sheets. This wall was con- 
structed In 1960 and had been used In the studies of Kefs. 43 through 
47. Although the circulating water was isothermal, a film buildup over 
the years prevented the test wall from being perfectly Isothermal. 


21 



Small segment-to-segment temperature differences existed, typically 
smaller than 0.2‘C, which resulted in scatter of the Stanton number data 
of approximately 2X. This scatter was not considered in the uncertainty 
analysis. The plate-to-plate streamwlse conductance for each gap within 
the preplate was measured during the heat flux meter calibration. These 
values were then used in the data-reduct ion program to correct for small 
plate-to-plate heat flows. 

The convex wall was constructed of b mm thick copper stock seg- 
mented into S cm lengths in the streamwlse direction. Lach segment was 
electrically heated, allowing steady-state measurement of the spanwlse- 
averaged wall heat flux by energy-balance techniques. Ualu reduction 
Included correction for: plate-to-plate heat conduction, losses to the 
support ussembly, and radiation losses. The plate-to-plate conductances 
were calculated. Utp conductances at the ends of the plates and tor the 
preplale/curve and curve/ recovery wall gaps were measured during prelim- 
inary tests where the temperature drop across the gap under study was 
made artificially large and the other heat flow paths were Insulated or 
controlled to zero AT. The power delivered to each plate, less correc- 
tions for other losses, was presumed to be conducted across the gap with 
the measured AT. clap- conductance uncertainties were Incorporated into 
the uncertainty analysis, they are typically small contributors to the 
overall uncertainty. The 14 segments of the curved wall were supported 
by ten c ire umt erent lal phenolic ribs. The ends of the ribs were held by 
a large aluminum frame. The curve was cut by first adding aluminum ribs 
(see Fig. 2-2) tor additional support then turning the entire assembly 
while cutting the 90° arc with a single-point cutting tool. After the 
machining operation, the alumliuiu support ribs and tide rail spacers 
were removed. The center used tor the machining operation was part of 
the frame and was the center of rotation of an arm used for traversing 
the various boundary layer probes. The heating elements tor each curved 
wall segment were two 4b cm lengths of AWJ F28 chromel w're embedded 
into parallel heater grooves and epoxy-bonded. At one end of each 
plate, the two embedded wires were connected by a large copper bus bar 
and at the other end were connected to the output terminal of a variable 
transformer. The overall heater resistance was about 8 11. 
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To have wel l-stablllzed power to the plates, the building power waa 
passed through two voltage regulators In aeries. The output oi the sec- 
ond regulator was connected to the Input of a step-down transformer. 
Tills reduced the voltage to typically ~ 40 volts AC. Next In line were 
powerstats that controlled each test plate power. This arrangement al- 
lowed Individual control of the power over a wide range with accuracy. 
All electrical power cables were enclosed Inside conduit to minimize In- 
terference with Instrumentation cables. A swlt hlng arrangement permit- 
ted the Insertion of a precision wattmeter intc each circuit. 

Constructed Into the test facility but not used for the present 
study was a system for Injecting 1 cm diameter discrete Jets ot air 
through the last 13 plates of the curved wall at a 10 Injection angle. 
Hiesi will be used in future curved-wall, film-cooling experiments [4H| 
which are a continuation ot the Stanford tlat-wall studies [43,44 r 45|. 
For the present study, the Injection holes were filled with balsa wood 
plugs before the surfa'e was sanded and polished. 

K.ich plate had three 1 ron-consi antc.n thermocouples distributed 
across the span tor redundancy and to detect spanwlse variation of 
temperature. For most cases, the output ot the three thermocouples was 
nearly the same, Indicating two-dimensionality, but lor cases where 
transition occurred within .he curved section, the convective heat 
transter coefficient was s lgr.l f leunt ly higher near the ends ot the 
plates where corner secondary flows promoted earlier transition. For 
these cases, a simple spanwlse heat-conduction model was necessary that 
could be used to estimate the centerline heat flux from the measured 
spanwlse average heat flux. This model was tested by comparison to 
measurements ot centerline heat flux using stlck-on heat flux meters 
(see Section 2.5.c). Although this correction was small for cases where 
transition was complete at the beginning of the curve. It was used In 
reducing the data of all the heat transfer runs. the uncertainty of 
this correction was Included In the overall uncertainties given in 
Appendix (1. The embedded thermocouples were laid Into milled slots so 
that there was a long contact length (L/D > 20) to minimize conduction 
error, following the techniques of Moffat [49 J. Hu* alimlnvmi trame sec- 
tion near the ends of the copper plates was heated with hot circulating 


23 



witter to minimize end losses. Hie copper tubes (or ducting this heating 
water were crushed against the aluminum with the side-rail blocks shown 
in Fig. 2-2. The side-rail, fraae, and tube assembly was tilled with 
hlgh-conductlvlty grease betore crushing. Hie large support d r m was 
heated with patch heaters that were controlled to a specified tempera- 
ture. 

The recovery wall was the same as the preplate wall, so that the 
test surface was symmetric about the center of .he curve. 

To minimize the radiation heat transler from the heated copper 
walls, they were sanded with progressively finer sandpaper then polished 
with commercial copper polish. After polishing, the surface was shiny 
enough that details of the surroundings could be examined In the reflec- 
tion on the surface. The surface was regularly polished to remove oxide 
buildup. With these precautions, the surface emisslvlty was held to an 
estimated 0.05 to 0.15. 

2.4 Inwt rumrui .it lon SuI lfydrodyn.ua toasuroaeiit* 

Mean velocity profile measurements were taken using a total pres- 
sure probe and wall static pressure ports. Hie outside diameter of the 
total pressure probe was 0.7 mm. The wall-static pressure and total-to- 
statlc pressure differences were read from a Statham model PM-97 trans- 
ducer calibrated to assure linearity to within ± 0.251 of full-scale 
output. Output was read with a Hewlett-Packard Integrating digital 
voltmeter model 2401C using a 33-second Integration time. Before each 
test, a two- point recal ibrat lon against a Comb 1st micromanometer was 
made. Pressure differences less than approximately 0.25 mm of water 

were read directly off the Combist micromanometer. In the curved 
region, the static pressure was read at the wall. Hie local velocity 
was then calculated from the formula: 

U - ~ (Pt(n) - Ps( n) ) 1/2 

where Pt(n) was the measured local total pressure and Ps(n), the 
local static pressure, was calculated *rom Psw as 
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Flow angle* were measured with a Conrad probe constructed with 45" 
bevels and an Indicator to measure relative angles. Cbnvergcnce angles 
In the boundary layer were referenced to the potential core flow 
direction at the same streamwlse and spanwlse position. Hie uncertainty 
ot the angle measurement was an estimated 1°. The differential pressure 
across the two tubes of the probe was sensed with the Statham PM-97 
transducer and read with the Hewlett-Packard Model 2401C Integrating 
digital vol tmeter . 


2.5 Instrumentation for lk*at Transfer Measurements 

a. Temperature Measurements 

All temperature measurements were taken with lron-constantan 
thermocouples, except for the boundary layer traversing thermocouple, 
which was ehroinel-constantan . Samples of the curved section lron- 
constantan thermocouples and the traversing chrome 1-const antan thermo- 
couple were calibrated against a Hewlett-Packard quartz Thermometer , 
which was accurate to within 0.02"C. The samples of lron-constantan 
showed the same calibration curve, at the beginning of the data-taklng 
and near the end, to within 0.08°C. The calibration curves and the 
uncertainty of calibration were Incorporated into the data- reduc t Ion 
program . 

The embedded thermocouples In the preplate and recovery wall were 
calibrated at operating temperature by referencing them to a calibrated 
sample of lron-constantan wire. This was done by taping the reference 
wire onto the copper plate section over the thermocouple to be calibra- 
ted, then pressing a 2.5 cm thick by 15 cm wide by 53 cm high layer of 
styrofoam over the thermocouples so that they were centered. Outside 
the styrofoam was a full-area patch heater which was controlled so that 
the temperature difference across the styrofoam was less than ~ 0.1C. 
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Alter about 20 minutes, the entire aaaeably reached theraal equilibrium, 
with no temperature gradient acroaa the copper and air gap separating 
the two thermocouple Junction*. The calibration tor each thermocouple 
In the preplate and recovery wall waa then taken and Incorporated Into 
the data-reductlon program. 

All the thermocouple wlrea were brought together Into an Isothermal 
zone at the back of the console panel, where they were connected to 
rotary thermocouple selector switches leading to the display o! the 
signal through a Hewlett-Packard digital voltmeter, Model 3465B. to 
avoid temperature gradients along the thermocouple wires, they were In- 
serted Into 1/8-ln. polyflo tubing for a portion of their length and a 
metal cable tray thermally Insulated on the Inside for the remainder. 
The Isothermal zone box was lined inside with 0 .8 mm copper sheets to 
conduct away local hot spots and Insulated on the outside with aluminum 
foil-backed rock wool Insulation. An aluminum tree-standing radiation 
shield was Inserted between the Isothermal zone box and power equipment 
In the lab. One diagnostic thermopile was Installed from corner to 
corner within the zone box to test tor thermal gradients. The output ot 
this thermopile was typically less than 3-4 mV and never larger than 9 
mV (~* O.l'C). All the iron-constantan thermocouples shared the same 
Ice bath reference Junction. 

Temperature traverses were made with a chromel-constantan thermo- 
couple probe constructed with 0.0b mm wire and designed for minimum 
conduction error (50]. This probe was calibrated against the quartz 
thermometer and measured temperature of the flowing stream with an esti- 
mated uncertainty of ± 0.0BC. Processing of the traverse data included 
correction tor viscous heating, following Moffat [51], and for the 
effects on fluid properties of temperature and humidity. The output was 
read on a Hewlett-Packard integrating digital voltmeter Model 2401 °C. 

b . Lm bedd ed Heat Plux Me ter s 

Sui face heat flux for each preplate or recovery wall section was 
measured with a heat flux meter installed between bakelite lamlnat .8 
behind each 0.3 mm thick copper plate segment. Lach meter Is 5 cm wide 
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and 0.4 am thick and la wound with multiple sllver-constantan thern. 0 - 
piles to measure temperature difference across Its thickness. 

The embedded heat flux meter calibration was done with a specially 
constructed calibration heater. Power to the heater was measured with a 
Sensitive Research Ga 1 vanometrlc type precision wattmeter. After cali- 
bration, the heat flux meter constant was known to an estimated uncer- 
tainty of ± 3X. Details of the calibration and the uncertainty estimate 
are given In Appendix C. Hie heat flux meter output was connected to a 
selector switch via shielded Instrument wire. The output was read on a 
Hewlett-Packard digital voltmeter model 3465B. 

c . Stlck-On Heat Flux Meter s 

Stlck-on heat flux meters were used lor spanwlse uniformity studies 
and for coarse verification of the Installed instrumentation. They were 
mul tl- Junction copper-constant an thermopiles laminated Into 0.08 mm 
thick sheets ot Mylar and were attached with double-stick plastic tape. 
Since they were detachable, they could be moved from place to place 
during a test. They measured local heat flux with an estimated uncer- 
tainty of 5-10X. Though they were calibrated with the same procedure as 
with the embedded meters ( App . C) , the calibration changed slightly with 
each sticking and removing cycle--hence the larger uncertainty. 

d . Electric P ower 

Power delivered to each plate of the curved section was measured 
with a single Sensitive Research galvanometric type precision AC watt- 
meter by switching it Into the power circuit of each plate, one at a 
time. The wattmeter was calibrated in a DC mode against a standard 
to within 0.05 watt. The wattmeter calibration was Incorporated Into 
the data-reductlon program. A correction was made within the data- 
reductlon program for the wattmeter Insertion losses. This correction 
Is discussed In detail In Ref. 44. 
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2.b i^ua lit li' Alton ot the Facility — Hydrodynamics 

A detailed discussion ot the hydrodynamic qualification ot the I low 
wun presented by Gil 11s and Johnaton (3b], important point s are summa- 
rized be'ow. 

The potential (low In the developing region waa found to be unltorm 
in velocity to typically i 0 . i *> percent, and the level ot turbulence 

intensity 

V u l2 /u. 

waa typically less than O.SX. The apanwlae variation ot the moment um 
thlckneaa at the beginning ot curvature waa typically leaa than 1 bi. 

Vlierever there la atreaallne curvature, there la a cross -at ream 
pressure gradient. It the flow la In a tunnel ot finite dimensions, 
tills pressure gradient will accelerate boundary layer tluld on the side 
walls t rom the concave surface to the convex. Streamlines on the convex 
surface then tend to converge toward the centerline as the secondary 

flows continue their spiral. Converging streamlines are inevitable; the 
designer of the curved tunnel must make the convergence acceptably 

small. Convergence angles, measured with a Conrad probe, were 2" over 
ttie central span of 13 cm and less than b* over a span ol 2b cm. 

Streamline convergence became perceptible alter about bO* ol curvature, 
and continued to grow throughout the remainder ot the curve, then per- 
sisted down tlie recovery wall. An Integral nu>mcntum balance using 

curved flow definitions ot the Integral parameters suggested by iion.tml 

| b2 ] was presented In Kef. 37 and showed closure to within bX lor the 
baseline case. The Influence ot secondary flow Is seen predominantly in 
the wake region ot the boundary layer and Is expected not to influence 
significantly the wall-measured data or extrapolations from tin* log 

region to the wall. 

2.7 Qualification ot the Facility— Mnat Ttanafar 

Hie flow In the developing region was tound to be uniform In 
temperature to within O.ObC. It was monitored continuously during each 
run and controlled constant to within O.ObC. Temperature control ot tiie 
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tunnel dir was achieved by dumping d portion ot the water circulating 
through the heat exchanger and replenishing with make-up water t rom a 
cold-water supply main. Hie spanwise uniformity ot the heated boundary 
layer at the Introduction of curvature waa checked by meaaurlng the 
local heat flux with itlck-on heat- 1 lux metera and by measuring the 
local enthalpy thlckneaa at various spanwise stations. Hie enthalpy 
thickness non-uni iorml t y was less than ± bZ , and the local heal- flux 
non-unitormlty was of the order ot , or less than the uncertainty of the 
meter . 

Within the curved region, thermocouples distributed across the span 
of each plate indicated that spanwise non-uniformity of the wall tem- 
perature was typically less than 0.07C and as low as 0.03C when preplate 
wall suction was applied. Spanwise uniformity of the heat flux within 
the bend was tested with the stlck-on heat-flux meters. In the first 
half of the bend, the heat flux was uniform to within the uncertainty ol 
the sllck-on heat tlux meter. In the last halt ot curvature, a spiral 
vortex driven by secondary flow spilling over the top ot each boundary 
layer fence was observed by using a tuit on a wand. This vortex was 
centered about 1 cm inboard of the fence and about 5 cm away from the 
wall. Measurements of the local heat tlux using the stlck-on heat-flux 
meters (Fig. 2-J) indicated that these vortices were augmenting t he heat 
transfer over a 2 cm span at the end of each segment. These vortices 
are not believed to have influenced t lie average heat flux significantly, 
what little Influence they may have had rhould have been accounted for 
by the spanwise conduction correction model incorporated Into the data- 
reduction program. In the afterplate, the fences were located Inboard 4 
cm from the walls to catch this vortex, break it up, and contain It be- 
tween the fence and the wall. The corner vortices are not believed to 
have Influenced the afterplate data significantly, since these data were 
taken with heat flux meters embedded at the wall centerspan. A conduc- 
tion analysis on the wall showed that the 3 mm thick copper plate sepa- 
rating tiie heat flux meter from the air averaged the heat flux over an 
effective span of only about 8 cm of the 40 cm fence-to-f ence total 
span . 
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Since two different techniques were being used to measure th* 
average heat flux, It seemed appropriate to measure, as a means of 
qualification, the heat flux on the preplate, curve, and recovery walls 
with a common Instrument, 1 .e . , a stlck-on heat flux meter. During a 
representative run, thl stlck-on heat flux meter was moved progress- 
ively from the last of the preplate to the first of the curve, the last 
ot the curve, and the first of the recovery wall. It was found that the 
built-in Instruments and stick-on meters agreed to within their uncer- 
tainties . 

Hie effect of small secondary flows on the wall-measured data or 
data found by extrapolation from the log region to the wall is expected 
to be minimal. An earlier design of the boundary layer fences allowed 
convergence angles In the recovery region that were about 50Z more than 
those for the final design. Stanton number data taken before and alter 
the modification were the same to within their uncertainty, indicating 
no significant sensitivity to small ^?condary flows. Secondary flows 
slightly Influence the wake regions of the thermal and hydrodynamic 
boundary layers, however, making an Integral energy or momentum balance 
along the tunnel centerline difficult. 

Further qualification of the Stanton number data was made by com- 
parison to the moan temperature profiles for the baseline cases. These 
profiles extended Into the conduction layer to n + values as low as 
4.0. Hiey followed the curve T + ■ Pr n + (see Fig. 3-3), indicating 
accurate values of wall and fluid temperatures, Stanton number, and skin 
friction. Also, local heat flux values were calculated from the 
temperature profiles as 



where T was taken a distance n away from the wall near n* ■ 6.5. 
This local heat flux was then used to calculate r ne local Stanton 
number, which was within 5Z of the wall-measured values for all the 
stations except the last station of the recovery region, where the 
gradient-calculated value was ~ 91 higher. 
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Hie Stanton ninhfr data were repeat able within their uncertainty 
band. The batu-llne case wmh the laat case run ot the 13 caaea discussed 
herein. Hie first qualified caae to be run waa oatenalbly the same 
case, but was completed eight months prior to the baseline case. Al- 
though the tunnel had been set up In many different coni lgurat Ions In 
the Interim, the Stanton number data for the two cases were the same to 
within 3Z. 


2.8 Hie Knerg^r Balance 

A tlnal all-inclusive test of the hydrodynamic Mnd heat transfer 
mcuHurcment a is the energy balance. Hi1b was performed by taking de- 
tailed wall measurements and velocity and temperature profile measure- 
ments while the rig was carefully controlled to the same steady-state 
conditions. Hie wall measurements taken at each wall segment are shown 
in Fig. 2. A. Hie energy balance data are tabulated in Appendix G as 

case 100779 and are also plotted as Fig. 2.!) in St versus Re. coor- 

2 

dlnates. Table 2.1 shows the excess 1 lowing energy in the boundary 
layer per meter ol span lor various at Teamwise locations: 



puli-l^) dn 



Hits was calculated for each of the spanwlse positions by: (1) inte- 
grating the velocity and temperature protlles and (2) Integrating tne 
energy Integral equation in the streamwlse direction using wall-measured 
values of heat flux. It all measurements were perfectly certain and the 
streamlines were perfectly parallel, the agreement between the two meth- 
ods would be exact . 

Hie overall closure was 93Z. This Is reasonable, the closure for 
the momentum balance for the same case [33] was 93Z, which Is considered 
good [37]. Closure tor the first 80Z of the curve is better than 93X, 
indicating secondary flow effects only on the profiles taken downstream 
of the curve. This secondary flow effect Is expected to be on the wake 
region of the profiles and is expected to have no significant eftect on 
the wall-measured Stanton numbers. 
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Tabic 2.1 
Energy Balance 


s <€■) 

-35 

10 

41 

61 

119 

Flowing 






Energy (KJ/aec-a) 


13' 

52° 

78* 


Froa Prof Ilea 

0.216 

0.475 

0.648 

0.720 

1.027 

Froa Wall Meaaureaente 

0.227 

0.506 

0.663 

0.746 

0.981 

Closure (Z) 


93 

99 

97 

93 
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Fig. 2-la. Plan view of the facility. 
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Fig. 2-lc. Location of measuring stations. 
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Fig. 2-3. 


Spanwlse variation In local heat flux near the end 
of the curvature. 
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Chapter 3 

THE EXPERIMENTAL RESULTS 


The following chapter flrat dlscuaaea a repreaentat 1 ve heat trans- 
fer run In aoae detail, then presents results of invest lgat Ions of sepa- 
rate effects: cases similar to the baseline esse but with varying values 
of a particular parameter. 

3.1 The Bas eli ne Case 

The baseline case Is a fully turbulent boundary layer responding to 
the Introduction of, then withdrawal of, convex curvature subject to 
uniform wall static pressure and teaperature. Data from this case are 
shown in Figs. 3-1 through 3-3. 

3.1.1 Stanto n Numb er versus Strea aw lse L ' stance 

Figure 3-1 shows the local Stanton numb, r versus distance In the 

streamwlse direction for a ratio of boundary layer thlckness-to-radlus 

of curvature of 6gg/K * 0.10. The wall static pressure coefficient on 

the test surface, C , was held constant to within ± 0.02 and the wall- 

P 

to-f reestream temperature difference, nominally 17*C, was held constant 
to within ± 0.6*C. The Reynolds number based on momentum thickness at 
the beginning of curvature was 4173, and the shape factor was 1.41, In- 
dicating a mature turbulent boundary layer. The uncertainty of the 
Stanton number data was typically 3.32, 4.82, and 3.52 for the upstream 
developing region, the curved region, and the recovery region, respec- 
tively. These typical values are shown In Fig. 3-1. The 952 certainty 
band for each data point Is listed with the data In Appendix G for all 
the cases discussed. The uncertainty analysis Is presented in Appendix 
A. Small differences In teaperature from segment to segment Introduce 
variations in local Stanton number which appear as scatter. These vari- 
ations are not considered In evaluating the uncerta.nty Intervals. 

Figure 3-1 shows that the effect of curvature on the heat transfer 
rate Is quite dramatic. The wall heat flux decreases ~ 152 when the 
curvature Is Introduced and continues to decrease within the bend so 
that, at the end of the 90* curve, the Stanton number is ~ 35-402 below 
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Che value that would be expected on a flat wall (daahed line). The 
recovery of the heat tranafer rate* on the flat wall downatruam of the 
bend la extremely alow. After bO cm of recovery length, the Stanton 
number la atlll ~ 1 3X below the flat-plate expected value. Rile beha- 
vior la the same aa that of the akin friction ahown . n Fig. 1-3. 

If Reynolda analogy held, the behavior of the turbulent tranaport 
of thermal energy (v't 1 ) would be similar to that of the Reynolda 
shear stress (-u'v n ) measured by Gillis and Johnaton (35), Fig. 1-7. 
The shear stress profile at the beginning of curvature la dramatically 
altered. The shear stress la essentially shut oil In the outer part of 
the boundary layer and reduced In the Inner layer, resulting in reduced 
wall shear and heat flux and a steepening ot the velocity and tempera- 
ture profiles at the wall. Note also, 1 rom the shear-stress profiles 
(Fig. 1-7), that the turbulent mixing activity returns very slowly to 
flat-plate values In the recovery region — hence the slow recovery of 
skin friction and Stanton number. 


3.1.2 Sta nton Nu mber versus Enthalpy Thickness Re ynol ds Number 
There exists a definite relationship [3b] between the Stanton num- 
ber and the enthalpy thickness Reynolds lumber tor a fully turbulent 
boundary layer on a flat plate with uniform wall temperature and static 
prcssur* and no unheated starting-length el feet. This relationship 1s 
shown in Fig. 3-2, along with the data from the present baseline case. 
Hie preplate data Join this correlation as the unheated starting length 
effect disappears, but, when curvature Is Introduced, the data begin to 
drop below the correlation. Shown also in Fig. 3-2 is the equivalent 
correlation for a laminar boundary layer. The slope of the laminar 
boundary layer correlation Is -1.0, which Is the same as the slope for 
the turbulent curved-wall data. This would lead one to think that sta- 
bilizing convex curvature is causing something similar to "relamlnarlza- 
tion,” a name used In connection with strongly accelerated turbulent 
boundary 'ayers (3b). Turbulence measurements by Gillis and Johnston 
[33] and others, show that the curved turbulent boundary layer is still 
very much a turbulent boundary layer, but one In which the production of 
turbulence Is contained to a thinner layer than if the streamlines were 
not curved. Their mean velocity profiles [33] also show no lasting 
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growth of the viscous sublsysr with convex curvature 


rherefore , 


though the St versus Re A ^ tr * ce similar to thst of s laminar 

boundsry layer, the boundary layer Is not laalnsr nor Is the sublayer 
growing significantly within the curve. In the recovery region, Fig. 3- 
2 shows the Stanton niaaber slowly returning to the flat-plate corre- 
lation. After a sufficient downstream distance, the curvature effect 
become* distant history and data should once again fall on the corre- 
lation (similar to the disappearance of the unheated starting length 
effect) . Comparing the data to the flat-plate correlation of St 

versus Re. seems a reasonable way of determining the degree of 
2 

recovery from curvature. The Re A ^ correlation of Fig. 3-2 Is based 
upon local conditions, whereas the Re x correlation of Fig. 3-1 > 

based upon the virtual origin of the turbulent boundary layer, which is 
permanently altered by curvature. Viewing the data In terms of St 


versus 


Re, 


allows one to see the effect of curvature without 


contamination by Reynolds number effects. Figure 3-2 shows that, at the 
end of the recovery section, the local Stanton number Is still ~ 20Z 
below the flat-plate correlation. 


3.1.3 Mean Temp er ature Profiles 

Figure 3-3 shows representative temperature profiles expressed In 
Inner coordinates. One profile Is taken from the preplate region, one 
at the beginning of the curve, four from the curved region, and four 
from the recovery region. The first profile in the recovery region was 
taken only 1.6 cm downstream of the end of the 90° bend. The profiles 
show that the effect of curvature Is to Increase the strength of the 
wake region. This is seen as a steepening of the temperature profiles 
In the outer part of the boundary layer. Curvature apparently does not 
result In a noticeable lasting change of the thickness of the conduction 
layer . 

Since the mean velocity profiles (Fig. 1-4) had a discernible log 
region, both within the curved region and the recovery region, and since 


*There appears to be a small, rapid Increase of the sublayer thick- 
ness, followed by a quick return to normal as curvature Is Introduced 
[3SJ . 
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local Reynolds analogy Is expected to be approximately true , * 
1.0 [31], one would expect the temperature profiles to have a log region 
as well. Figure 3-3 shows that there is Indeed a log-linear region of 
the temperature profiles, though it becomes quite small within the curve 
(as with the mean velocity profiles). "Hie thermal law-of-the-wal 1 ex- 
pression which was used is: 

Pr + 

T - 13.2 Pr ♦ In yy y (Ref. 3b) (1) 

Tills expression was derived tor a normal tlat-wall turbulent boundary 
layer. Since the Stanton number and skin friction coefficients are 
measured Independently, the only unknown available with which to torce 
Eqn. (1) into the data is the average turbulent Prandtl number. Hie 
log- region dashed lines In Fig. 3-3 are 1 rom Eqn. (1), which fits ac- 
ceptably well the log-linear region or each prollle, except the one 
taken 1 .6 cm downstream of the end of the curve. Tills boundary layer 
Is apparently not In equilibrium Hie mean velocity profiles used In 
deducing /2 for Eqn. (1) (and making a recovery correction to the 
Indicated temperature) were measured in the heated boundary layer with a 
pitot tube. Hie turbulent Prandtl numbers used to match Eqn. (1) to 
the data are plotted versus streamwlse distance In Fig. 3-4. Hie intro- 
duction ol curvi.'ure causes a slight Increase in the turbulent Prandtl 
number, as predicted by K.M.C. So [31]. Hie rise in Pr^ when curva- 
ture Is Introduced Is partly due to the disappearance of a weak unheated 
starting length effect and partly due to the Introduction of curvature, 
so that the curvature effect would be slightly smaller thun the total 
effect shown in Fig. 3-4. Shown also are Pr { values found with the 
same technique, but for the case 6^/R “ 0.05. It Is expected thai the 
unheated starting length effect Is Insignificantly small by the first 
point (S ■ -35 cm) for the 6^g/R • 0.05 case and by ~ 30° of cur- 
vature for the ■ 0.10 case. 

Hie effect of the withdrawal of curvature on the turbulent Prandtl 
number Is quite dramatic. ttien the turbulent activity begins to in- 
crease near the wail and propagate outward (as shown in Figs. 1-b and 
1-7), the transport of thermal energy apparently increases faster than 
the transport of momentum. It's not unreasonable to expect a violation 
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of Reynolds analogy with this rapidly changing boundary layer, since the 
two quantities lu the ratio • M /« |( are fundamentally different, one a 
scalar and the other a vector. It la remarkable , however, that this 
violation peralata ao tar downatream. Hie trend la almllar to that Been 
In the strongly decelerating, but attached, turbulent boundary layer and 
la repeated lor the A^/R • 0.05 caae. 

3.1.3 Kner^y Ha lance 

Since both detailed pro! Ilea and detailed wall data were taken for 
thla example caae, all the terma In the energy Integral equation (Appen- 
dix B) are available. 

A check on the data waa made by comparing the growth ot the en- 
thalpy thlckneaa Reynold a number, aa calculated t rom the energy Integral 
equation, to the enthalpy thlckneaa Keynolda lumber a tor each profile. 
The comparison, Klg. 3-5, ahowa that they agree to within ~ *) X. Hu* 
largeat contributor to the lack ol agreement la the Influence ot the 
amall secondary t lows on the wake ot the boundary layer, aa previously 
discussed. Secondary t Iowa convert thermal energy toward the tunnel 
centerline on the convex wall, resulting In higher measured enthalpy 
thicknesses than predicted by Integrating the energy equation along t lie 
centerline streamline. Agreement within la reasonable 13/). Second- 
ary 1 Iowa are not expected to a 1 gn 1 1 1 cant 1 y Influence the wall-measured 
Stanton number or the portions of the profiles between the wall and the 
wake . 

The following sections present several studies made to investigate 
separate parameters. Descriptors ot the baseline caae and cases in t lie 
tol lowing studies are presented In Thble 3-1. 

3.2 The Effect ol the Initial Boundary layer Thickness 

Three cases were compared lor this study: the baseline case (A^/K 
(0-0) - 0.10), a case with (0-0) - 0.05, and a third case with 
Aqq/R (0-0) - 0.02 (see Klg. 3-b) . Hie data upstream ot the curve dif- 
fer slightly for the three cases because they have different virtual 
origins (Reynolds number effect). Within the curve, the data (Klg. 
3-6) are similar tor the three cases. In these coordinates, showing the 
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1 1 r ■ t evidence ot the way In which convex curvature organizes the data. 
The pariuseter Agg/R, labeled the "strength of curvature," has been used 
as a curved-wall flow descriptor (2). Figure 3-b shows that, within the 
range 0.02 * 6gg/K (0 ■ 0) * 0.10 and for the conditions tested, there 
la only a weak dependence on Agg/K (0 - 0) as viewed In these coordi- 
nates . 

If the streamwlse distance were scaled on the radius of curvature 
or shear layer thickness, A g ^, the above conclusion about the effect 
of Agg would again hold. It the streamwlse distance were scaled on 
the boundary layer thickness (either the local value or the value at the 
beginning ot the curve), one would reach a quite different conclusion. 
If the abscissa were S/Agg the case with larger A would have the 
abscissa compressed relative to those with smaller A. One would then 
conclude that, for cases of larger Agg(t)-0)/k, the effect of curvature 
Is more than for cases of smaller Agg(t)»0)/R because the slope of 
St versus S/Agg(t)-0) is steeper. At this point, there is no strong 
reason to choose one scaling length over the others, so the data are 
left In terms of the dimensional quantity ”S". 

Figure 3-7 shows these data recast in the coordinates St versus 
enthalpy thickness Reynolds number. Kventually the Stanton number data 
for all three cases follow a line of -1 slope in these coordinates. 
The data of the Agg/R ■ 0.10 and 0.05 cases rapidly approach the 
-1 sloped line at the start of curvature. Because the cases of Agg/R 
■ 0.10 and 0.05 immediately enter the asymptotic state, they can be 
considered to represent strong curvature. It appears that, once the 
curvature Is sufficiently strong, the parameter Agg/R is of minimal 
further significance. Clllls and Johnston (35] found that A g ^/R " 0.03 
for both their cases (Agg/R - 0.10 and 0.05). Therefore, one could 
expect that. If no other parameter were changed from its test value 
but Agg/R, cases with Agg/R > 0.03 would have their region of non- 
zero u'v* compressed to 0.03 R, whereas cases with Agg/R < 0.03 
would show continued growth of the boundary layer until Agg/R ■ 0.03. 
The data of the Agg/R - 0.02 case Is quite different from the above 
"strong" curvature cases. The asymptotic state is reached much more 
slowly, presumably because the initial boundary layer thickness Is less 
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than the eventual shear layer thickness. The abrupt Introduction of 
curvature results In an immediate decrease In Stanton nmber In this 
weak curvature case similar to that seen for the strong curvature cases, 
however. The mlnlmixn Stanton numbers tor the three cases of Fig. 3.7 
are 3O-40X below the flat-wall correlation. 


Within the 60 cm recovery length, the Stanton data for each case 

returns to within 10-131 of the Re. correlation values. The rate of 

a 2 

return of the Stanton number to the flat-wall turbulent boundary layer 
correlation seems to be proportional to the difference between the local 
Stanton number and the flat-plate Stanton number for the same enthalpy 
thickness Reynolds number. 


3.J The Effect of U 


pw 


Figures 3-8 through 3-10 show the results of a study of the effect 
of Figure 3-8 shows Stanton number versus distance In the stream- 

wise direction for the baseline case, U^ w ■ 13 m/s, and for two other 
values of U pw : ? m/s and The two hlgh-veloclty cases re- 

spond nearly the same (when corrected for the Reynolds number effect), 
while the low-velocity case seems to be approaching closely to the same 
state In the last 30° of the curve. Since U pw appears in the denomi- 
nator of Stanton number, this means that the wall heat flux at any posi- 
tion on the curve Increases almost linearly with U pw » The shear stress 
profiles from Gillls and Johnston (Fig. 1-7) display an asymptotic state 
for 0 > 60°. The cases of higher velocity (112779, 070280 In Fig. 3-8) 
show the Stanton number data rapidly locking into this curved flow 
state. The lower velocity case (011280) shows the data approaching the 
other data more slowly, and about 60° of the curve Is needed for the 
low-velocity data to attain the curved boundary layer condition. In the 
last 30° of bend, one can find a very weak U pw dependence as viewed In 
these coordinates; the data of the 2b m/ s case (112779) is 8X below the 
baseline data (070380) and the U ■ 7 m/s case (011380) data are 14X 

r * 

above . 


The velocity that would exist at 
boundary layer. 


the test wall, 


were there no 
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Figure 3-10 shows a three-dlaenalnal plot In which l) pw ia used aa 
the third axis. Thia allows one to viaualize the surface acre easily 
and provides space for insertion of a fourth case which haa the interme- 
dlate velocity but a boundary layer thickness of 6^/K • 0.05 (see 
Table 3-1). in Fig. 3-10, the low- and high-velocity cases can e^ch be 
compared with an Intermediate case of nearly the saae boundary layer 
thickness at the beginning of the curve, separating the effect of 6^^/R 
from the effect of 

The following paragraph exemplifies the difficulty in viewing data 
in these coordinates. The Stanton number data at any given S-locatlon 
within the last 30“ of bend are nearly the same for the three cases, the 
lower-velocity cases yield lower Re x values and therefore higher pre- 
dicted flat-wall Stanton numbers. The flat-wall Stanton lumber plane is 
shown in Fig. 3-10 as a nearly horizontal plane. The lower-velocity 
cases would therefore show a stronger effect of curvature in terms of 
St/St o . The minima within the curve are 50X, 42X, and 39X below the 
flat-plate correlations for U ■ 7, 15, and 2b m/s , respectively. 
If, in the asymptotic curved boundary layer, the Stanton number is de- 
pendent upon only local conditions, as now appears to be true, the his- 
tory of the development of the boundary layer upstream of the bend is 
not important. Comparisons to flat-wall values in these coordinates, as 
done above and commonly done with curvature studies, are not appropri- 
ate. The Reynolds number effect muMt be separated from the curvature 
effect . 

Figure 3-9 shows that the data within the curved region, for all 
three cases, eventually follow a -1 relationship with Ke^ . Careful 
examination of the introduction of curvature shows that, for the low- 
velocity case, the data approach this -1 slope in the first 30-40X of 
the curve. In the mid-velocity case, the data quickly reach the -1 
correlation (as observed in the previous study, Fig. 3-7) at the begin- 
ning of the curve. The data of the highest-veloc i ty ca e drop below 
the -1 correlation, then approach this correlation in the first 20-30° 
of curvature. This seems to be the principal effect of U p W : t0 c ^ an 8 e 

the shape of the St versus Re^ curve without changing its even- 

2 

tual -1 slope. 
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Figure 3-8 or 3-10 shows that the increase In Stanton masher In the 
recovery region Is faster In the lower-velocity cases. The reason for 
this becomes evident when the data are viewed in the coordinates of Fig. 
3-9. The rate of return of Stanton number within the recovery region 
appears to be proportional to the distance from the flat-plate correla- 
tion (a first-order lag in these coord lnates) . 

3.4 Effect of Freestream Acceleration 

Figures 3-11 and 3-12 show the combined effects of streamwlse cur- 
vature and acceleration. The curved region Is the only region where 
there is freestream acceleration, K * 0.0 for the preplate and recov- 
ery regions. In their discussion of the accelerating boundary layer, 
Kays and Moffat [38] point out that, for the Asymptotic accelerated 
boundary layer, the moment van thickness Reynolds number reaches a 
constant value while the enthalpy thickness Reynolds number continues to 
Increase. In the following study, the momentum thickness Reynolds 
number at the start of curvature was within 10% of the value they gave 
for the asymptotic state corresponding to the acceleration parameter 
K. Therefore, the momentum bnoundary layer was expected to lock into 
the asymptotic acceleration state at the start of curvature. 

Three cases are compared for this study: K ■ 0.0, 0.57 x 10~ b , 
and 1.25 x 10 b . The boundary layer thlckness-to-radius of curvature 
ratio for all three cases was about 0.05. The K ■ 1.25 x 10 b case 
of Fig. 3-12 shows a different Stanton number trace within the curved 
region than seen in the preceding figures or in the case of milder 
acceleration. There apparently is an acceleration effect which is as 
important as the curvature effect. 

A considerable amount of work has been done on accelerating bound- 
ary layers on flat walls, (e.g., Ref. 38). The following is a very 
brief review: Examples of the effect of acceleration on smooth, flat- 
wall boundary layers [38] for K • 0.57 x 10 -b and K - 1.45 x I0~ b arc 
shown In Figs. 3-13 and 3-14, respectively. Figure 3-13 shows that 
acceleration with K ■ 0.57 x 10 b does not signlflcntly influence the 
heat transfer rates. Figure 3-14 shows that acceleration with K - 
1.45 x 10 -b results in a dramatic dropoff of the Stanton number. 
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This Is predominant ly s result of the acceleration's thickening the 
viscous sublayer or conduction layer, that la, it la predominantly an 
Inner-region effect. Figure 3-14 shows also that, for the strong 
acceleration cases studied, the slope of the Stanton number curve in 
tnese coordinates Is between -0.23 and -1.0. 

Hie combined case of curvature and acceleration (Fig. 3-12), shows 
a combination of two effects which act on different parts of the bound- 
ary layer, acceleration on the inner and curvature on the outer. The 

• lope of Stanton lumber versus enthalpy thickness Reynolds number for 
the case of strong curvature with strong acceleration of K - 1.23 * 
10 b (Fig. 3-12) is approximately -2.0, which shows that both effects 
are significant and the two effects are additive. Hie case of K ■ 0.37 

* 10 _b (Fig. 3-12) shows an eventual -1 slope, but there appears to 
be an effect of this lesser acceleration at the beginning of curvature. 

A careful review of the data for this study shows an Interaction of 
the effects of curvature anc‘. acceleration. Figure 3-13 shows that, on 
the flat wall, acceleration with K ■ 0.57 * 10 b should not show an 
effect of acceleration. Figure 3-12 shows that K ■ 0.57 * 10 b 
results In an effect of acceleration at the Introduction of curvature 
that apparently disappears as the expected -1 slope appears. Figure 
3-14 for the flat wall shows that K ■ 1.45 * 10 b results in a -0.5 
slope; a K value of 2.5 * 10 -b would be required to have a -l .0 
slope. Figure 3-12 shows that K ■ 1.25 * 10 b results In a -2 
slope. It appears, therefore, that the curved boundary layer is more 
sensitive to the growth of the conduction layer than is the flat-wall 
boundary layer. This Is similar to the Interaction between acceleration 
and wall suction [38 j. 

3 .5 Effect of Un heated Starting Le ngth 

Figures 3-15 and 3-lb show the effect of the location of the begin- 
ning of heating with respect to (1) the virtual origin of the turbulent 
momentum boundary layer and (2) the beginning and ending of curvature. 

The first comparison is to show that the unheated starting length 
Is not Influencing the studies already discussed. Cases 070280 and 
060480 have heating beginning ~ 63 cm upstream of the beginning of cur- 
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vature. Case 070280 haa heating beginning 134 cm downstream of Che 
virtual origin of Che turbulent momentum boundary layer, whereat* caae 
060480 haa heating beginning 24 cm upatream. foe two caaea are there- 
fore quite different In unheated atartlng length. Figure 3-16 ahowa 
that, even for caae 070280, where the unheated atartlng length la moat 
aevere, the effect of the unheated atartlng length dlaappeara quickly 
and la gone before curvature la Introduced, the Stanton number data are 
on the correlation. Dlfference8 In the Stanton number data from caae to 
caae for all the caaea dlacuased thua far, therefore, cannot be attrib- 
uted to the unheated atartlng length eftect. 

The second comparison la to show the effect ot starting a thermal 
boundary layer Inside of a well-developed turbulent momentum boundary 
layer In the vicinity of the curved section. 

Case 113079 of Figs. 3-13 and 3-16 shows the effect of curvature 
where the thermal boundary layer began at the start ot curvature. Wien 
curvature is introduced, the production of turbulence is confined to the 
Inner part of the boundary layer, effectively reducing the scale of the 
layer. In the cases where heating begins within the curve, it is rea- 
sonable to assume that much of the added thermal energy Is contained 
within the Inner (active) region. If this is so, the thermal energy Is 
not allowed to mix effectively to the outer region and the Stanton num- 
ber drops fast at first, then quickly levels off. Wien viewed in the 
coordinates of Fig. 3-16, the data of the two cases where heating began 
upstream of the curve (cases 070280 and 060480) eventually reach the 
presumed asymptotic state characterized by the line of -1 slope. In 
the case where heating began at the beginning of curvature, the data 
again attain this -1 sloped line as well. Using the above scenario, 
this condition was eventually reached when the inner layer received a 
sufficient amount of thermal energy to establish an asymptotic tempera- 
ture profile. Fig. 3-16 shows that the recovery for this case is 
similar to that of cases 070280 and 060480 characterized by a rate of 
return proportional to the difference between the local Stanton number 
and the correlation. 

Case 030280 had the same hydrodymamic boundary layer as 113079, but 
the heating began at the beginning of the recovery section. In this 
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case, the thermal boundary layer began to grow a» the shear layer of the 
Glllls and Johnston model began to recover to flat-wall conditions. The 
data finally attained a trajectory where its return to the correlation 
of Fig. 3-16 was at about the same rate as seen for the other cases, and 
was slower than would be expected In a normal flat-wall boundary layer. 

There Is evidence from this study that the state of the flow 
outside of the shear layer Is not Important In the recovery process. A 
comparison of cases 070280 and 030280 shows that they recover similarly 
after the unheated starting length effect of case 030280 disappears. 
However, the distributions of the thermal energy within the recovering 
layers for the two cases are quite different. Case 070280, In which 
heating began within the preplate region, would have a considerable 
amount of thermal energy as residue outside ”^ a l"* Case 030280, In 
which heating began at the beginning of the recovery section, would have 
minimal thermal energy beyond Since they recover in a similar 
fashion, this residual thermal energy must not be Instrumental In the 
recovery process. 


3.6 The Eff ect of t he Mat urity of the Momentum Boundary Layer 

The results of a study showing the effect of the maturity of the 
momentim boundary layer at the start of curvature are shown in Figs. 
3-17 through 3-21. The descriptor to Identify the separate cases is the 
momentum thickness Reynolds number at the start of curvature. * For this 
study, Re^ ( b "0) varies from 222, Indicating a lamlnar-to-early tran- 
sitional boundary layer, to the baseline case of 4173 Indicating a fully 
turbulent boundary layer. 


Figure 3-17 shows that, if the boundary layer Is sufficiently tur- 
bulent, (cases 070280, 011380, 010680, and 012480), the convex curvature 
becomes an effective organizer of the data. Though the preplate data 
for the four cases are radically different, the traces of Stanton number 
versus streamwlse distance within the curved region (Fig. 3-17) are 


A simple model has been devised for extrapolating the measured 
starting profile data to the beginning of the curve. This model is 
presented as Appendix D. 
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remarkably similar. A plausible explanation of this can be made by 
using the "shear layer" model proposed by mills and Johnston (33]. 11 
the boundary layer Is sufficiently turbulent, the thin layer (~ 301 of 
the original boundary layer), which will become the shear layer, has 
well established turbulent mixing activity. Mien curvature Is Intro- 
duced, the turbulence production Is confined to this Inner region, and 
the remainder of the boundary layer Is of minor Importance. A mature 
turbulent boundary layer would have large-scale turbulent mixing. In 
the asymptotic curved turbulent boundary layer, the scales become adjus- 
ted to the shear layer thickness. Within the shear layer, then, the 
asymptotic states of both transitional and fully turbulent cases are 
essentially the same. The history is contained In the outer layer, 
which apparently has a minor Influence. iltese late transitional and 
turbulent cases, Fig. 3-18, display the slope -1, considered to be the 
asymptotic curved turbulent boundary layer state. 

Figure 3-18 shows that the recovery data for cases of late transi- 
tional or turbulent boundary layers are similar to those of the mature 
turbulent boundary layer. 

In contrast to the above four casus are those In which the boundary 
layer Is laminar or early transitional at the start of curvature. Figs. 
3-19 and 3-20, characterized by the nearness of the data to the laminar 
flat-plate correlation. For the nearly laminar cases (e.g., 062880 and 
062380), the small amount of turbulent activity that may be present Is 
not sufficient to force transition to continue In the face of stabiliz- 
ing curvature. The data remain near the laminar correlation until the 
layer has matured: the convex curvature has therefore delayed transi- 
tion. Mien transition finally begins, within the bend, It appears to be 
retarded by curvature. Figure 3-20 shows the transition, as seen by an 
Increase of Stanton number toward the turbulent correlation, to be slow 
until the stabilizing convex curvature Is removed and then to progress 
quickly . 

A notable difference from the two cases discussed above Is case 
060980, a case with extremely low velocity and momentum thickness Rey- 
nolds numbers. In this case, although there appears to be some effect 
of curvature, particularly at the introduction, the curvature does not 


32 



seem to retard the transition as strongly as It did In the other casus 
and transition Is nearly complete by the end of the curved region. 


Hie laminar and early transitional cases in Fig. 1-20 show that the 
curvature did not prevent transition. The Introduction of the stabilis- 
ing curvature resulted In a momentary return toward the laminar correla- 
tion, but transition effects dominated the trajectory within the curved 
region. At the very end of the curve, cases Ub0980 showed a trend simi- 
lar to that observed In fully turbulent curved boundary layers (the -1 
slope). This may be due to the growth of the boundary layer, as dis- 
cussed above, In conjunction with the emergence of turbulence activity 
as the dominant transport mechanism. Hie uncertainty tor these transi- 
tional data was as large as ± 10-131 (see App . 0 tor details). 


Figure 3-20 shows the recovery data for lamlnar-to-ear ly-transl- 
tlonal boundary layers. Hie data return quickly to the turbulent flat- 
wall boundary layer correlation. Hie most significant effect of curva- 
ture was the delay and retardation ot transition. 


Figure 3-21 shows the combined plot of Figs. 3-17 and 3-19 as a 
three-dimensional view of Stanton number versus sire, unwise distance and 
momentum thickness Reynolds number at the beginning ot the '•urve Re^ 
(b«0). Hie figure shows three distinct regimes. In the first, exempli- 
fied by the most distant curve, transition dominates the trajectory and 
curvature Is only a mild retardant. In the second, the next two curves 
in order of Increasing Re, (0*0), curvature delays and retards the 

I 

transition process and, when curvature Is removed, the boundary layer 
returns to "normal" quickly. In the third regime, there is the typical 
response of a mature turbulent boundary layer to convex curvature. 


53 




•rrurjiMwi'.r pittan<'| (cm,' 

ri(. J-l. Stanton number vinui »ttf«n%Ur distance for the 
baseline c» i# (070280) - Turbulent <*9<)/R( A "0) ■ 
0.10, U pv - 14.8 m/s , K » 0.0. 
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rig. 3-2. Stanton number versus enthalpy thickness Reynolds 
number for the baseline case (070280) - Turbulent, 
■ 0.10, • 14.8 m/s , K • 0.0. 
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Fig* 3-5. Comparison of Re^ from the energy integral 

equation and Re^ - from the measured profiles for 
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Fig. 3-13. Tho effect of mild f ree-st ream acceleration (from 
Kays and Moffat [38]) - Stanton number versus 
enthalpy thickness Reynolds number. 



Stanton 

numfcer 


/ u -*: 

Enthaloy thickness *eyno ds rurter, I — 


Fig. 3—14 . The effect of strong free-stream acceleration (from 
Kays and Moffat [38]) - Stanton number versus 
enthalpy thickness Heynolds number. 
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FIr. 3-15. The effect of unheated starting length - Stanton 
number versus Mn.un.ltr distance. 
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FIr. 3-16. The effect of urheated starting length - Stanton 
number versus enthalpy thickness Reynolds number. 
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Fig. 3-17. The effect of the nwturltv of the momentum boundary 
layer - Stanton number versue strearawise distance 
for the late-transit lonil and turbulent cases. 
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Fig. 3-18. The effect of the maturity of the momentum boundary 
layer - Stanton number versus enthalpy thickness 
Reynolds number for the lat e-t rans It lonal and 
turbulent cases. 
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for the laminar and early-transit lonal cases. 
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3-20. The effect of the maturity of the momentum boundary 
layer - Stanton number versus enthalpy thickness 
Reynolds number for the laminar and early- 
transitional cases. 
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Fig. 3-21. The effect of the maturity of the momentum boundary 
layer - Stanton number versus enthalpy thickness 
Reynolds number and momentum thickness Reynolds 
number at the start of curvature. 





Chapter 4 

HEAT TRANSFER PREDICTIONS 


On*? of the objective# of the curvature project waa to develop a 
turbulence model for curved-wall heat tranafer predictions that could be 
uaed for practical engineering calculation#. Thla la, In fact, the moat 
uaeful form In which to convey the result* of these curvature experi- 
ments • 


4.1 State of the Ar^t 

Prior to the beginning of the curved-wall studies at Stanford, the 
moat prevalent curvature model was one which stemmed from a suggestion 
by bradshaw [14] and was formulated and teated by Johnston and hide 
(41)]. Recently, Gil 11a and Johnston (35], with the benefit of their 
recovery data, proposed a second model consistent with their "shear 
layer" Interpretation ot the curved turbulent boundary layer. both 
models are discussed below. They modified the distribution of the min- 
ing length across the boundary layer to account for streamwlse curva- 
ture . 

Other, higher-order , models have been developed by Launder et al. 
[ 34 ] , Irwin and Smith [39], and Gibson [33]. These models are quite 
good tor mild curvature. Although higher-order models give more detail 
of the turbulence structure, they are generally more difficult to run 
and less useful as engineering tools. In order to make the results of 
the modeling effort most useful to Industry, It was deemed desirable to 
keep the computation scheme slmple--hence the mixing- length approach. 

The mixing-length hypothesis uses an empirically prescribed length 
scale to relate the mean velocity gradient to the Reynolds shear stress, 
through the formula: 


- uv 


2 3U 
T n 



(4-1) 


A typical distribution of mixing length "1" for a flat-wall boundary 
layer Is: 
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I • «n |^l - exp(-n*/A*) j (4-2) 

the modified van Driest model (36) for the Inner portion of the boundary 
layer, where x, the Karman conatant , la 0.41. For the outer portion 
of the boundary layer, the typical flat-wall mixing length distribution 

has 

t - 0.085 fi 99 (4-3) 

The smaller of (4-2) or (4-3) is used. The parameter A + Identifies 
the thickness of the viscous sublayer and Is affected by acceleration/ 
deceleration or blowlng/suct Ion . The value tor a flat-wall boundary 
layer with uniform freestream velocity and no transpiration is A* ■ 25. 

The curvature model suggested by Bradshaw (14) modified this flat- 
wall distribution for curvature according to the relationship: 

l - i o (l - UR1) , 

where t 0 Is the local flat-wall mixing length and Ke Is the local 
gradient Richardson number defined by 

Ki - 2S(1 + S) , 

where S, the stability parameter, Is 



The parameter p Is an empirical parameter generally of the order 10. 
The effective radius of curvature, ^ett* cofflel1 from a first-order lag 
model proposed by Johnston and Fide (41): 

d(l/R ff ) 
eft 

dx 

Ro is the wall radius of curvature. This mixing-length model was In- 
vestigated in detail by Johnston and Elde (40,41). Their recommenda- 
tion for the value of the constant 6 was 6.0. Figure 4-1 shows its 
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performance for the baseline case. The prediction was Initialized with 
starting profiles taken 35 cm upstream of the beginning of the curve. 
The test wall was isothermal and U was uniform. No correction for 

p W 

the effect of curvature on turbulent Prandtl number was made. The re- 
sponse to the Introduction of curvature was predicted quite well, but, 
after the curve, Cj/2 and St recovered far too quickly. If the tur- 
bulent Prandtl number were decreased In the recovery region according to 
the curve of Fig. 3-4 , the recovery process would be even faster. This 
model was devised before recovery data had been taken. Since It is 
merely an empirical fit of curved and rotating wall mixing lengths, It 
is not surprising that the model's performance In the recovery region 
was poor. 

G1 1 1 1 8 and Johnston [35J improved the recovery prediction by devis- 
ing a model in which the main effect of curvature was to confine the 
turbulent motions close to the wall and allow the previously existing 
large-scale motions to decay. In the recovery region, they modeled the 
regrowth of the inner boundary layer as a normal boundary layer growing 
on a flat plate within a velocity gradient. Their model follows .* 

The mixing length in the outer portion of the boundary layer was 
modeled as: 


£ « 



6 


* 

si 


where 6^ is the width of the active shear layer and 6^ is the dis- 
placement thickness of the boundary layer integrated out to the shear 
layer thickness ^ 8 i • The constant of proportionality was found to be 
0.10. The shear layer thickness in the flat-wall preplate or recovery 
wall van found by extrapolating the linear portion of the shear stress 
profile to uv • 0, as shown in Fig. 1-8. Within the curve, Gillls and 
Johnston use a proposal from Gibson (42), that there is a critical value 
of the stability parameter, S, defined as: 


S - U/R 
b TIJ7TH 


See Ref. 35 for more details. 
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above which Che shear stress cannot sustain itself. Gibson found that 
this critical value was approximately 0.17. Gillls and Johnston found 
that a critical value of 0.11 optimized the prediction of their data. 

The model used in the following predictions of the heat transfer 
data is essentially the sam* model as presented in Ref. 35, with the 
exceptions of two changes made to the estimation of 6 g ^ within the 
curve . 

To describe the first change, it is necessary to define the follow- 
ing two par ame t e r s : 

Y The n-distance away from the wall where S ■ S 

crit tru 

The n-distance away from the wall where the extrapo- 
lated shear stress equals zero. 

In the Gillls and Johnston model, 

£ - 0.10 ( 6 ,- 6 .) 

si si 

where 6 , ■ Y . in the developing or recovery regions and 6 , ■ 

Slot SI 

Y in the curved region. With their model it was found that the 

crit 

switchover at the beginning of curvature quickly altered the shear 
stress profile driving Y gl down to a small fraction of Y orit . But, 
as the change in the shear stress profile began to alter the mean veloc- 
ity profile, Y 8 ^ began to rebound, passing up through and finally well 
above Y This rebound eventually affected the mean velocity and 

temperature profiles, causing the calculated values of Cf/2 and St 
to increase in a manner much like recovery from curvature, though still 
within the curve. The critical stability argua^nf states that, in the 
asymptotic curved state, the extrapolated shear stress. Y g |, equals 
the n-distance where S ■ S crlt , ^ cr it* 11111 * n ^is model, Y g j could 
be well above Y cr j t , especially, it was found, when 6^^/R was small. 
The change in the model was to construct a "controller” that would not 

allow Y , to grow much above Y . _ . The new model calculated the 
si crit 

6 s i to be used in computing the curved-region mixing length profile as: 

6 sl " Y crit wnen Y sl < Y crit 
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Y 

4 .1 • V crit (‘-O' 2 -o( T ~ - 1-0 ) 

' ' crit " 

wh *" v .l > v ctlt but > 0-33 • 

It Is a stiff "controller" that drives 6 , down when Y , > Y, . 

8* 8 1 C r 11 

trying to force Y^j ■ ^ cr ^ t for the asymptotic curved boundary layer. 
Note that when Y g ^ < Y cr j t the new instruction has no effect. 

The second change to the model was developed when trying to predict 
the three cases of differing 6 99 /R. It was found that the optimum 
value of S cr j t was a t unction o f 6 99 /r. Fitting the optimum values 
for cases of <5 99 /K ■ 0.10, 0.05, and 0.02 with a smooth curve gave; 

< 4 99 /R)l/3 
S crlt " 5 

Since the 6 99 /R used in the model Is the local value and not the value 
at the start of curvature, is continually updated as the bound- 

ary layer thickens. This correlation gives the correct trend with 
6 99 /R tor the cases studied. It is not known, however, that it would 
do equally well it 6 99 /R were changed from case to case by changing 
R while holding 6 99 (0-O) constant. 

In the following sections, predictions are made of many of the heat 
transfer runs presented earlier. In these predictions the turbulent 
Prandtl number war. not changed as a result of curvature. Figure 3-4 
shows that, within the curve, the change of turbulent Prandtl number is 
small and in a direction that would decrease the Stanton numbers 
slightly below the following predictions. 4 mild Improvement to the 
model could be made, then, by incorporating a turbulent Prandtl number 
correction for the curved boundary layer. One possibility is the model 
developed by R.M.C. So [31] discussed in Chapter 1. 

Figure 3-4 shows that there seems to be an effect of recovery on 
the turbulent Prandtl number. The following predictions, which make no 
such correction, are expected to show too slow a recovery. At pres- 
ent, no model exists for making this correction. Profiles of local 
v T t T /u 1 "v T across the recovering boundary layer would be valuable in 
developing a model that accurately accounts for this effect. 
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The present model can be quite temperamental. The calculation of 
Y b1 la done by extrapolation of a linear portion ot the shear stress 
profile to the wall. Alter some distance downstream ot the beginning ot 
curvature, the shear stress profile develops Inflections that can cause 
Y h j to Jump a considerable distance 1 rom step to step. This usually 
results In variations In the slope ot St versus Re^ (or s) , os- 
cillating about what appears to be a reasonable value. Uccasionally the 
program will extrapolate to an extremely small Y r1 due to an inflection 
very near the wall. When this happens, the solution goes irreparably 
out of control. Because ol this, the program can be quite frustrating 
and expensive to run. An improvement might be to calculate the shear 
layer thickness as an Integral thickness; 


Y , Is a function 
si 


o 

of / y dn , 

J w 


where n Is the n-distance when the shear stress goes to zero. Thin 


would stabilize the model considerably 


The tollowlng sections demonstrate how well the above model with 
the above changes predicts the heat transfer cases discussed in chapter 
3. 


4.2 Predicting the l.ttect ol Curvature tor the Baseline Case 

Figures 4-2 and 4-3 show the prediction ot the baseline case. The 
Introduction ot curvature rapidly decreases the skin friction and Stan- 
ton number. For this case the effect Is strong enough to overshoe t and 
oscillate tor some time before settling down to the asymptotic curved 
state. The prediction of both the skin friction and the Stanton number 
within the curve Is quite good. Note that a Fr t correction would de- 
crease St slightly near the end of the curve. The recovery is quite 
accurately predicted. Note that the Pr t correction would increase the 
predicted St In the recovery region. Tb try to assess Just how much 
the correction would increase St, the ratio of St/(Cj/2) was calcu- 
lated from the profile data for this case t rom Chapter 3, then multi- 
plied by the predicted C^/2 to get an estimated Stanton number with a 
recovery correction on Pr^. This corrected Stanton number Is shown in 
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Figs. 4-2 and 4-3. With this correction, the prediction and the data 
are quite close. Ideally, the prediction curve of Fig. 4-3 would follow 
the data with a -1 slope. Hie strong Introduction results in a 
steeper slope. As the shear layer thickens, Y g ^ continues to grow, 
and the slope decreases. When Y g ^ becomes larger than Y crit , the 
"controller" decreases 6 g ^, trying to drive it to Y cr it* ®een as a 
downward trerd near the end of the curve. Runs made with sustained 
curvature show that the model predicts a continued decrease In St with 
nearly a -1.0 slope lr these coordinates. The shape of the recovery 
Is approximately correct, characterized by a rate of return roughly pro- 
portional to the distance from the flat-wall correlation. 

Hie following sections discuss the predictions of cases with dif- 
fering values of a single parameter. Table 4-1 shows the case numbers 
and figure numbers that constitute each separate-effect study. 

Table 4-1 
Prediction Runs 


Effect 

— 

Case 

Parameter 

Value 

— 

Figure #' s 

699 / R 

070280 

0.10 

4-2, 4-3 


022680 

0.05 

4-4, 4-5 


060480 

0.02 

4-6, 4-7 

U pw 

112779 

26.4 m/s 

4-8, 4-9 

070280 

14.8 m/s 

4-2, 4-3 


011380 

7 .0 m/s 

4-10, 4-11 

K 

022680 

0 

4-4, 4-5 


051080 

1.25 * 10 " 6 

4-12, 4-13 

Unheated 

022680 

Base 

4-4, 4-5 

Starting 

113079 

A 2 /R(0-O) - 0.0 

4-14, 4-15 

Length 

030280 



A 2 /R(0-9O°) - 0.0 

. 

4-16, 4-17 
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4.3 Prediction of Cfr ses with Dl f faring Init ial Bo undary Layer Thlckneaa 

Figures 4-2 through 4-7 show the ability of the program to model 

cates of differing boundary layer thickness. The "predictions" are 

good, even for the thin boundary layer case where 6gg/K ■ 0.02. The 

term "prediction" might be misleading, especially In this study, because 

values of S , _ were chosen that would optimire the comparison of cal- 

c r 1 1 

culated and measured Stanton number and skin friction. No adjustment to 
the model has been made for subsequent studies; therefore, sections 4-4 
through 4-b can, more honestly, be called predictions 

Note that the preplate prediction of Stanton number (Figs. 4-3, 4- 
5, and 4-7) is progressively lower than the correlation for thinner 
boundary layer cases. "nils is probably due to the Reynolds number ef- 
fect. In Case 060480, where <5^/R(t>-0) - 0.02 (Figs. 4-b, 4-7), the 

momentum thickness Reynolds number was 4b4 at S * -35 cm and estimated 
to be 1724 at the start of curvature, indicating a late-transitional and 

early turbulent boundary layer. The turbulence model was developed for 

a fully turbulent boundary layer and apparently underpredicts the pre- 
plate data for the thin boundary layer cases. Hie modeling is expected 
to Improve as the boundary layer thickens downstream. The predicted 

Re* at the end of the curve for Case 0b0480 was 3200. 

6 2 

The predicted recovery of the thinner boundary layer cases appears 
to be slightly too rapid. Case 022680, where 6g^/R(t»“0) - 0.05 (Figs. 
•-4, 4-5), shows a faster recovery of skin friction than seen in the 
measured values and a recovery of Stanton number similar to that seen in 
the data. Figure 3-4 shows a decrease in Pr t within the recovery re- 
gion for this case, however, so that a Stanton number prediction with 
a Pr t correction would show too strong a recovery. 

4.4 Prediction of Cases with Differing U^ w 

Figures 4-2, 4-3, and 4-8 through 4-11 show predictions of cases 
with differing U . There seems to be a slight U effect that is 

r" r " 

not modeled. Hiis effect apparently begins at the introduction of cur- 
vature. Since the present model is based mainly on data from the asymp- 
totic state, it is not surprising that there are some deficiencies near 
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the beginning of curvature. One might expect that the strength of the 
introduction of curvature might scale as 



4.5 Pred iction o f Ca ses with Free atre am Acceleration 

Figures 4-12 and 4-13 show the prediction of the combined effects 
of strong tree-stream acceleration and curvature. The prediction models 
the -2 slope (Fig. 4-13) quite well In the early region. Apparently, 
however, the shear stress profile has developed an Inflection, and the 
extrapolation technique Is jumping from one 6^ to another quite dif- 
ferent one. This is seen both near the end of the curve and In the 
recovery section. The model proposed in Section 4.1, where 4 j is 
calculated by integration might correct this problem, and It Is expected 
that the prediction of the acceleration case would be quite good. 

4.6 Prediction of Cases with Di fferin g Unheate d St a rtin g Lengths 

Figures 4-14 and 4-15 shw«w the prediction of the case In which 
heating begins at the beginning of curvature inside a mature momentum 
boundary layer. Figures 4-lb and 4-17 show the prediction of a similar 
case in which heating begins at the start of the recovery region. The 
figures show an excellent prediction of the two cases. 
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STANTON NUMBER 


0 0040 


0 00 10 


0 OO-'O 


0 0010 


0 0000 


Fig. 



STKEAMWISE DISTANCE («'M ) 

4-1. Prediction of the baseline case with the .Johnston 
and Elde model. 
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b 

£ 

X 

* 


s 

5 


fc 


OOCIO 



STREAMKlSE DISTANCE (CM I 

Fig. 4-4. Prediction of Case 02268C, f^/Rf^O) * 0.15: 

Stanton number versus »t teamwise distance. 



ENTHALPY THICKNESS REYNOLDS NUMHER 

Fig. 4-5. Prediction of Cate 022W0, • 0.05: 

Stanton nvaiber vrrtui entlialpv thickness Reynolds 
nuaher . 
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STREAM VISE DISTANCE (CM) 

Fig. 4-6. Prediction of Case 0404B0, ; 99 /R(A-n) - 0.02: 

Stanton number versu» streamwlse distance. 



ENTHALPY THICKNESS REYNOLDS NUMBER 

Fig. 4-7. Prediction of Case 0S0480, fl 99 /R(6«0) - 0.02: 

Stanton number veraus enthalpy thickness Reynolds 
number . 
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S7AKT0N Ni'HBlfc 



STK»*yWIM t'lhTAS I (CM) 

FI*. 4-R. Prediction of C*»e 112779, ’ pw ■ 26 m/a : Stanton 

number vrraui atrramvlar distance. 
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:.n fie iseis 



S'l h l> M • , t 0I5TAN » (CM | 

rig. 4-10. IT ed let Ion of Case 011580, t' • 7 »/»; Stanton 
number versus at ream* lee distance. 




INTMAIM TMICKNFSS RlYNoUri ROMMBK 
rij. 4-11. Prediction of Case 011580. - 7 ■/■: Stanton 

number versus enthalpy thickness Reynolds number. 
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STANTON N MPfR 



STN»AM»ISI UlSTANi * (CM) 

Fig. 4-12. Predict Ion of Caee 0510*0, P • 1.25 * 10* 6 : 
Stanton nuolrr vtruJi atteanrwl»e dlstanie. 



number versus enthalpy thickness Reynold* number. 
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STKKAMWIhK DISTANCE (CM) 

Fig. 4-16. Prediction of Cane 050280, A 2 /R(0«90°> - 0.0: 

Stanton number versus streamwise distance. 



ENTHALPY THICKNESS REYNOLDS NUMHER 

Fig. 4-17. Prediction of Case 030280, Ai/Rf -90") - 0.0: 

Stanton number versus enthalpy thickness Reynolds 
number. 
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Chapter 5 

CONCLUSIONS AND RECOMMENDATIONS 


Tills study presents results of experiments on the effect of convex 
curvature on heat transfer through turbulent and transitional boundary 
layers. For Lite baseline case, the effect of curvature was separated 
from other effects such us streamwise accleration, which usually ae- 
companies curvature in practice. Miere there were combined effects, 
they were carefully controlled, i.e. , where streamwise acceleration 
witli curvature was to be tested, acceleration was with a constant K - 
v/u' 


pw 


dU /ds. 
pw 


Hie results of the 15 cases studied here provide considerable in- 
sight into the effects of curvature on turbulent boundary layer heat 
transfer. Hiey complement the detailed hydrodynamic measurements of the 
previous investigations, e.g., Gillis and Johnston [35], So and Mellor 
[lb] , etc. Hie studies of separate effects isolated several Important 
parameters and revealed some new characteristics of curved boundary lay" 
ers (e.g., the -1 slope in Stanton number versus enthalpy thickness 
Reynolds number coordinates). Also, they checked the robustness oi 
hypotheses, like the "shear layer hypothesis," derived in single-run 
experiments by subjecting them to a broad range of conditions. Most 
important to the modeler and, eventually, the designer, these tests 
provide a broad data base with which to verify future prediction models. 


5.1 Cbncluslons 

Hie main conclusions are as follows : 

1. Hie effect of convex curvature on heat transfer rates is signifi- 
cant for all the. cases studied. Reductions of Stanton number of 20-50% 
below the expected flat-wall values were observed. Hie effect of convex 
curvature is seen immediately at the beginning of curvature and contin- 
ues to grow throughout the curve. Hie recovery of the heat transfer 
rate on a downstream flat wall is extremely slow. After 60 cm of recov- 
ery length, (15-20 boundary-layer thicknesses), the local Stanton number 
may still be 15-20% below what would be expected without upstream curva- 
ture, at the same enthalpy thickness Reynolds number. 
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2 . 11k* boundary layer eventually roue ho a an asymptotic curved oondi- 

ta ^ 

tlon characterized by St « Re^ . Ttie ~l slope Is considered to be 
significant. It may be indicating that the overall heat transfer from 
the wall has become conduction-dominated as in a laminar boundary layer 
or a strongly accelerated boundary layer. 

3. Separate effects studies on this one curved surface show that the 
boundary layer thickness, free -stream velocity, and maturity oi the 
boundary layer affect the response to the introduction of curvature and 
the rate at which the asymptotic state is approached, but not the even- 
tual asymptotic state characterized by the =1 slope. 

4. Studies of combined stream wise acceleration and streamwise curva- 
ture effects show that they both can be significant and that they seem 
to be additive. Cbnvex curvature seems to increase the boundary layer’s 
sensitivity to acceleration, as has previously been observed when accel- 
eration and wall suction were combined. 

5 . Studies of the effect of convex curvature on laminar and transi- 
tional boundary layers show that this stabilizing curvature delays tran- 
sition from laminar to turbulent . Mien the curved boundary layer is 
near-laminar or early-transitional, the effect of the curved region is 
soon "forgotten," and the recovery data rapidly return to flat-wall 
expected values. Mien the curved boundary layer Is late-transitional, 
recovery is slow, similar to that of a mature turbulent boundary layer. 

b. Turbulent Prandtl number deduced from the thermal law-of- the- wall 
or calculated as St/fC^/i) increases slightly (10-lb*) within the. 
curved region and decreases rapidly in the recovery region. At the end 
of the recovery section, the log-law deduced Pr t is ~ 0.50 and con- 
tinuing to decrease. 

7. The "shear layer" model proposed by tiillis and Johnston [35J seems 
to be an effective model with which to view the effect of convex curva- 
ture on turbulent boundary layers. 

8. Hie curvature prediction model presented by Gillis and Johnston 
(35] and modified herein predicts the present heat transfer studies 
reasonably well, although it must still be considered "developmental." 
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5.2 Recommendations for Future Mirk on the Convex Problem 

Tills program has answered several important questions about the 
curved wall convective heat transfer problem. It dealt with the scaling 
and sensitivity questions and provided a broad data base to build upon. 

Ai Important next step is to study in detail the difference between 
the effects of curvature on the diffusivity of heat and momentum. This 
involves taking profiles of the shear stress u T v T and the turbulent 
transport of heat V't' to give profiles of Pr t . Development of a 
four-wire anemometer probe at Stanford is now nearly complete, and such 
data should be available soon. Of the four wires, three are for measur- 
ing the three components of velocity and the fourth is a temperature 
sensor. Prediction of the recovery of Stanton number will require de- 
velopment of a turbulent Prandtl model, which, in turn, will require 
that we know more about the processes than the present data can reveal. 

Another important separate effect study not done in the present 
experiment would be of the effect of radius of curvature. It has been 
assumed, but not shown, that the appropriate parameter for the strength 
of curvature is fiqq/R. This study has indicated that affects 
only the rate at which the asymptotic state is approached. The ques- 
tions that remain are: (1) If the radius is changed but dqq/RCO-O) is 
held constant, does the approach trajectory to the asymptotic state 
change? (2) Does the nature of the asymptotic state itself change? 

Other important questions remain about the curved and recovery 
regions. These are: If the curved wall were allowed to continue indefi- 
nitely (without influence of secondary flow) , would the conclusion that 

-1 

the asymptotic state is represented by St « Re. remain true? Or, if 

2 

the recovery wall were much longer, what would be the continued trajec- 
tory of Pr^ versus s (Fig. 3-4)? 

Still unanswered are questions about the combined effects of convex 
curvature and transpiration, or surface roughness. These are key ques- 
tions to the designer. A hydrodynainically smooth wall may become tran- 
sitionally rough when it becomes convexly curved and the scale of the 
boundary layer decreases. It has been shown that convex curvature or 
wall suction increases sensitivity to streamwise acceleration. How 
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sensitive to acceleration would a curved boundary layer with wall suc- 
tion become? 

Another Important question to gas turbine blade designers is: How 
does convex curvature influence the jet-freestream interaction for film- 
cooled blades? it most likely will be learned that the optimum blowing 
rate is curvature-dependent* 


* 


Research activity now under way at Stanford. 
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Appendix A 

THE UNCERTAINTY ANALYSIS 


Processing the data included the calculation of each reduced 
quantity and the uncertainty of that quantity (output in Appendix G) . 
The uncertainty can be interpreted as the band above or below the stated 
value within which it is believed the value would lie, with 20:1 odds, 
or 9b% of the time. The uncertainty tells how closely each stated value 
listed in Appendix G and diseussed in the body of this report can be 
believed. The uncertainty analysis also tells the designer or data** 
taker how much effort should be given to the measurement of any parties 
ula** quantity. The uncertainty analysis was brought on-line early in 
the design phase of this project and used in that capacity to guide the 
design. 

One. note: the uncertainty interval applies to any one data point. 
Mum many data points are taken over a region where there are no discon- 
tinuities, i.e., the preplate, the curved wall, or the recovery wall, 
the mean of the data is known with more certainty than any one data 
point; integration or averaging tends to decrease the uncertainty. 
Though the uncertainty presented is typically 3-5 a, the uncertainty of 
the mean of the data is much smaller. 

Tine uncertainty analysis built into the Stanton number data- 
reduction program was one that, as a first step, calculated all the 
sensitivity coefficients, 


S 


P' 


3St 

3P 


6P' 


the influence on the Stanton number uncertainty of the uncertainty of a 
particular parameter, P* . This is done by first rerunning the data- 
redu>’tion with that parameter changed by, typically, 1 X, to get the par- 
tial derivative, then multiplying by the uncertainty of P' . After the 
sensitivity coefficients was calculated, the overall uncertainty was 
calculated using Kline and MeCllntock [A-1J : 


6 St 



A-l 



The following table show* the parameters included in the uncer- 
tainty analysis and their estimated uncertainty.* 

Table A-l 

Input to the Uncertainty Malyaia 


Parameter 

Uncertainty 

Embedded, calibrated thermocouple output 

3 Mv 

Calibration constant for thermocouples in pre- 
and after-plate 

30% 

Plate-plate gap heat conductance 

10% 

Heat flux meter signal 

20 Mv 

Wattmeter reading 

0.05 w 

Heat flux meter calibration constant 

3% 

Plate area 

0.46 ram^ 

Anbient pressure 

7.6 mm Iig 

Relative humidity 

1 % 

Free-streara temperature 

3 mv 

Dynamic pressure difference 

0.13 mm H2O 

Emissivlty 

0.02 

Wattmeter calibration constant 

0.015 w 

Plate heater electrical resistance 

0.015 ft 
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(Table A-l (cont.) 


Plate Ifeater and delivery line electrical 

resistance 0.015 0 

Plate heater and delivery line electrical 

resistance, including transformer 0.015 ft 

Resistance of the ammeter circuitry 0.005 0 

Tunnel static pressure 0.13 mm 1^0 

Spanwise heat transfer correction 20% 

Support section temperature 4 Mv 

Recovery factor of free-stream thermocouple 0.04 

Thermocouple calibration constant 0.10 

Saturation pressure of water 0.14 atm 

Saturation density of water 0.0016 kg/nr* 

Inductance of ammeter circuitry 0.003 ft 

Resistance of voltmeter 100 ft 

Resistance between wattmeter and transformer 0.005 ft 

Preplate-curved section gap conductance 10% 

Curved section/afterplate gap conductance 10% 



Tin* uncertainty analysis was written into both the profile data- 
reduetion program (listed In App. K) and the Stanton number-reduction 
program (listed in App. F) . The total uncertainty of each quantity is 
given, along with the reduced quantity In Appendix G. hore detailed 
output giving the maxima and minima of the sensitivity coefficients 
within the developing region, curved region, or recovery region for each 
input to the uncertainty analysis is also available. Tills is the de- 
tailed intormulion used in the rig design. At example of this is given 
on the following page for the baseline case. 


Re ferenoca 

A-l. Kline, S. J., and F » A. MeClintoek, "Describing Uncertainties in 
Single -Sample Experiments," Moo ban i calEng live ering, Jan. 1953. 
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Appendix B 

Curved-Will Integral Parameters and the Momentum 
and Energy Integral Equations 


Wien there is streamwise curvature, it becomes necessary to alter 
the definitions of the integral parameters somewhat and, hence, the form 
for the momentum and energy integral equations. The following curved- 
wall definitions : the displacement and momentum thickness presented by 
Honami {52,37] . . » sed in the present study: 


bisp lac erne n t_ Th i c k nes s of Curved- Wa ll B o und ar y layers 

s i 

Moment um Thickness of Curved- Will Bound ary La yer s 
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Wiere it is assumed that Up * Up W /(l+kn), the integral momentum 


equation becomes [52,37]; 
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The definition of the curved-wall enthalpy thickness used in 
following study was derived by starting with the basic equation: 


the 



p U (i -i ) dn 
00 p w 00 



pUU-iJ dn 


( B-4) 


After substituting the estimate Up * Up W /(l+kn) and integrating: 
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The simplest form of the energy integral equation is: 
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Uniting in terms of the Stanton number defined as 

St - q*7[o U (i -i )] 

‘o 1 » pw w 00 * 

Eqn» (B~b) (assuming constant free- stream enthalpy and density) becomes: 
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Taking the derivative, Eqn . (IH>) becomes 
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Equating the two terms that (B-7) and (B-8) have in common leads to: 
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on the flat wall, k « 0, and (B-9) reduces to 
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Equation (B-9) is the integral energy equation used in the present 
study. Note that the appropriate equation for high-speed flows would 
have accounted for changes in and p^ . 



In the Stanton number data-reduetion program, Eqn. (B-9) is inte- 
grated over a wall strearawise distance 6s to determine the change in 
enthalpy tnickness as; 
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Appendix C 

Calibration of the Heat Flux Meters 


Calibration of the heat flux meters required designing a heater 
that would asure a precisely known and carefully controlled one- 
dimensional heat flow through the heat flux meter, This device is shown 
in Fig. C-l. 

Hie heater consisted of a unifora-power-dlstribution patch heater 
that was 15.2 cm wide by 45.7 cm high, backed by a 2.5 cm thick piece of 
styrofoam. On the outside of the styrofoam was another patch heater, 
and beyond that another 2.5 cm thick piece of styrofoam. The power to 
the heating pad sandwiched between the two styrofoam sheets was con- 
trolled so that the temperature difference between the two heaters was 
small (typically less than 0.1°C), preventing any back-loss from the 
heater adjacent to the copper plates. The input signal to the control- 
ler came from an iron-constantan thermopile. 

The power supplied to the heater adjacent to the copper plates was 
adjusted to give a representative heat flux. The input power was regu- 
lated AC power controlled with a transformer and was measured with a 
galvauometric-type precision AC wattmeter. The temperature of the cop- 
per plates was controlled to a representative value with the built-in 
controller on the hot-water circulating system. After one-half hour, 
the entire assembly reached steady-state, and the average heat flux of 
the heating pad was equal to the local heat flux near the center of the 
assembly . 

The heating pad covered six copper plate segments. It was found 
that heat flux meter calibration constants, K, measured with the plate 
at the two center positions agreed with one another to within less than 
2%, and usually less than IX. 

th 

Writing down the energy balance for the i plate and dividing by 
the plate area gives: 


K i HF i + s i-i<VW + wvw 


where 



q£ ■ the measured heating pad power divided by pad area, with 
small corrections for wattmeter insertion and back- 
conduction; 

HFj * the output signal of the heat flux meter; 

S i-1* S i+1 " 8 a P conductance between the 1 plate and the 

upstream or downstream plate. 

Mien the test plate was in one oi the two center positions, there 
were either two or three guard plates on each side, and the neighboring 
plates were at nearly the same temperature. Therefore , only a very 
small correction was needed in the above equation for strearawise conduc- 
tion. Tills correction was made using estimated values oi S. The only 
unknown, then, was the desired calibration constant, . The measured 

calibration constant was then input to the Stanton number data-reduction 
program along with a correction for temperature, as suggested by the 
manufacturer, 

K i " Kref,!* 1 " (T r 80°F)/700°F) 

Once the calibration constant was found, the plate was moved to one 
side of the heater and the above test was rerun- -then to the other side. 
Mien on the side of the heater, plate-plate temperature differences were 
larger and the streamwise conduction terras became significant. The two 
S-values could then be calculated from the two tests, since there were 
two equations and two unknowns. Hie calibration device was designed for 
very accurate calibration of the heat flux meters but, unfortunately, 
did not give accurate S-value measurements . The estimated uncertainty 
on the S-values was 10a. Tills terra is not a large contributor to the 
overall uncertainty, however. 

An uncertainty analysis on the heat flux meter calibration using 
representative values was performed with the technique of Kline and 
McClintock [ G— 1 ] . It was found that the uncertainty on K was ~ 
2.7%. A value of 3% was used in the Stanton number data-reduction 
program. 

Refer en ces 

C-l. Kline, S. J., and F. A. McClintock, "Describing Uncertainties in 

Single-Sample Experiments," Mechanical Engineering , Jan. 1953. 
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Appendix D 

Simple Model lor Ordering the 
Transitional Boundary layer Cases 


In the study of the effect of maturity of the momentum boundary 
layer, It was decided to order the eases according to Re^O-G) . 

Initializing profiles were taken 35 cm upstream of curvature, and, in 

the case of a fully turbulent boundary layer at this point, standard 

correlations were used to estimate Re^^ (O^O) . Mien the boundary layer 

£ 

at S « "35 cm was laminar or transitional, the extrapolation became 
more difficult. The following simple model was used for this extrapo- 
lation. 

The following correlations were used for the fully laminar or £u ? ».y 
turbulent boundary layers ; 

Turbul ent : 

Re,r = 0.037 Re 0,8 , 

6 2 o 

C f /2 - 0.0295 Be“°‘ 2 

i o 

C^/2 S dS2/ds 
leading to : 

“*0.2 s 

/ds - 0.0129 Re, 

6 2 

and 

« 2 « 0.037 s 0,8 

' u pw ' 

(s is the distance from turbulent boundary layer virtual origin.) 
Laminar : 

Re o - 0.664 Re 1/2 

S 2 8 

C f /2 » 0.664 Re~ 1/2 
I 8 

C f /2 5 d<$ 2 /ds 

D-l 



leading tos 


d^/d* * 0*441 Re^ 


*2 " 0 




(a is distance from laminar boundary layer virtual origin.) 

It was also assumed, where needed, that transition was always in 
the interval 

400 < Rex < 800 

6 2 

Wien Reg (s » -35 cm) was less than 400, the average Reg^ between 
s - -35 era and the beginning of transition was estimated as 


Re, - (Re. (s * -35 cm) + 400)/2 
°2 2 

from which the average d6 2 /ds was calculated. 

^2 - 1.0 

- 0.441 Re* 

de 5 2 

Since 6 2 (s - -35 cm) and 5 2^ Re <S 2 are availab: l e » tlie 8 “ 

distance where Rex * 400 can be estimated as: 

d 2 

/ d6 2 

s(Re^ ■ 400) ■ -35 cm + A6^ / 

Within the transition region, 


Re. » (400 + 800)/2 
°2 


^0.441 Re“ 1,0 + 0.0129 Re~°’ 25 ) + 2 


and s(Reg^ * 800) is estimated. 
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Hie virtual origin of the turbulent boundary layer can be ealcula™ 


ted from 


s(Rcg » 800) and ^(Reg “ 800) 
u 2 2 

from which Rc^Cu-0) can be calculated, 

ib 

Of course, each case is different, some may be transitional at s " 
-35 cm and others at s » 0 cm, requiring variations on the above 
scheme . 


D-3 



Appendix E 

Listing of the Profile Program 


The following listing is that of the Fortran program used to reduce 
the velocity and temperature^prof ile data. 
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1 

z 

5 
4 

a 

6 
7 
6 


9 

10 

It 


10 

13 


14 

15 

16 
17 


1« 

19 

?.a 

21 

22 

23 

24 


♦MATFIV LIST 

C 

C VELOCITY t TEMPERATURE DATA REDUCTION PROGRAM 

0 REVISE!) JULY 1980 

C 

DIMENSION MLNs 16 ) , AU3t 60 ) ,A'!B( 60 ) ,ZA( 60) 

DIMENSION! U(60J.Tlo0liVl603.T1fc0 1,Z(6O),CF(60) 

DIMENSION TUI 60 1 iDURt 60 UQ 1 60 ) iQSPADl 60 ) , YMl 60 ) 

DATA PGT,NUAX,1« , MA^,AUINf ,AUB,ATfi,ATO, AVISO, U,V,T/3*0, 303*0./ 
DATA Y»Z /12G#0./ 

DATA DC HP . BELT 99 , VPW, VIST » TO/5* 0 . / 

INTEGER ZUM.P'G 
ZLitroo 
c 

C *1# READ PROFILE INFORMATION 

C 

C NO0RT 

C IDT 

G 

C 

C ITG 

C 
C 
C 
C 

C MODE 

C 

C 

C REG 

e 
c 
c 

READ (5 t IO> NPORT , *P7» ITC , MODE , REG 
10 FORMAT (5X6) 

C PGM RETURNS HERE «»* TO TERMINATE PGM, READ BLANK CARO 
IE IMPORT. EG. 0. AND. M8DE.FQ.0) GO TO 670 
C 

C #C» DEAD PUN NUMBER AND r.UN INEO. , A FORMAT COL 1“72 
C 

READ (5.50) (KRNd). 1=1,36) 

50 FORMAT (3SA2) 

C 


DOST NUMBER FOR PROFILE MEASUREMENT 

PRESSURE TRANSDUCER FOR VELOCITY PROFILE 

IPT=1 RGMDIST MICRGfJANGUETER 

IPG 2 EPCGCURE TRANSDUCER 

iHfPMDCOUPLE PROBE FOR TEMPERATURE PROFILE 

STG-’O NO AGCGNPAHilMS TEMPERATURE PROFILE 

XtG-f GOOSE ML CM THERMOCOUPLE PROBE l 6. 5 Mil OFFSET) 

2T6-C EIP3U THERMOCOUPLE PROBE (21.5 MIL OFFSET) 

STC"3 01 AGHKELL PROBE (5.0 MIL OFFSET) 

1 T PE OF DATA REDUCTION REQUIRED 

MODL^O COMPUTATION OF INTCGPAl PARAMETERS, U+.Y+.XVO 
MODE cl COMPUTATION OF INTEGRAL PARAMETERS ONLY 
f LOW REGIME - F GET CALCULATION OF CF/O AND XVO 
HEG=0 TURBULENT 
REG"; LAMINAR 


C *3* 

C 

C 

C 

C 

C 

C 

C 

109 

C 


c 

C #4# 
C 
C 
C 

c 

c 


150 


200 

250 

C 

C *5* 
C 


READ TUNNEL CONDITIONS FOR VELOCITY PROFILE 

TRECCV MAINSTREAM TEMPERATURE, I-C TC (MV) 

PSTAT WALL STATIC PRESSURE (INCHES H2Q) 

PAMD AMOUNT PRESSURE ( INCHES HG CGKRECTEO TO 32 DEG FI 
RH AMBIENT RELATIVE HUMIDITY (PERCENT) 

H 'ROPILE MEASUREMENT LOCATION FROM START OF CURVATURE-IN 

RAD RECIPROCAL OF PA0IUO OF CURVATURE (1/IHCHES) 

READ 15,100) TRECOV, PSTAT, PAM3.RH.X, RAD 
FORMAT (6FJ0.G) 

PHOT TUBE O.D, = 0.020 INCH 

CD=0.025 

OFFSET=0 .014 


READ VELOCITY PROFILE DATA 

Y MICROMETER READING ( INCHES ) 

U TRANSDUCER READING ( INCHES OF WATER) 

READ BLANK CARD AS LAST PROFILE POINT 

DO COO 1=1,50 

READ (5,150) Ym.Utl) 

FORMAT(CFIO.O) 

NMAX=I“1 

IF (UCI).LE.O. ) GO TO 250 

CONTINUE 

CONTINUE 

READ LIMITS FOR COMPUTING SHEAR STRESS FROM CLAUSER PLOT 



25 

26 

27 


26 


29 

30 

31 

32 

33 

34 

35 


36 

37 


38 

39 

40 

41 


42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


C 

C 

C 

C 

C 

C 

C 

C 

C 


C *6* 

C 

C 

C 

C 

C 

C 

C 

901 

C 

C #7# 

C 

C 

C 

C 


YPLO LOWER LIMIT OF LOC REGION 

SUGGESTED VALUE 75 FOR HI RED ELS FLOWS 

22 FOR VERY LOW REDEL2 FLOWS 
40 FOR CURVE OR RECOVERY 
YPHI UPPER LIMIT OF LOG REGION 

SUGGESTED VALUE 125 TO 200 FOR HI RE0EL2 FLOWS 
62 FOR VERY LOW REDEL2 FLOWS 
90 FOR CURVE OR RECOVERY 


READ (5,150) YPLO, YPHI 
IF (ITC.EW.O) GO TO 370 

READ TUNNEL CONDITIONS FOR TEMPERATURE PROFILE 
READ(5»901 ) PAM , PSTA , PON »RL ,TO , STANT 

PAM AMBIENT PRESSURE (INCHES HG CORRECTED TO 32 DEG F) 
PSTA WALL STATIC PRESSURE (INCHES H20) 

PCN MAINSTREAM TOT PRES-KALL STATIC PRES (IN H20) 

RL AMBIENT RELATIVE HUMIDITY (PERCENT) 

TO PLATE TEMPERATURE, I-C TC (MV) 

STANT LOCAL STANTON NUMBER 
FORMAT(7F10,0) 

READ TEMPERATURE PROFILE DATA 
2 MICROMETER READING (INCHES) 

T THERMOCOUPLE READING (MV) 

READ BLANK CARD AS LAST PROFILE POINT 


DO 360 1=1 ,50 

READ (5,151 ) 2(1), T(I) 

151 FORMAT (2F10.0) 

NZMAX=I-1 

IF (T(I).LE.O.) CO TO 370 
360 CONTINUE 
370 CONTINUE 
C 

C WRITE ALL RAW DATA 

WRITE (6,5) (NRN(I), 1=1 ,36) 

5 FORMAT (1 HI ,/9X, 'CURVATURE PROJECT #H# NASA-NAG-3-3 ##* CONVEX 

1 TEST SECTION'/I OX, 'VELOCITY AND TEMPERATURE PROFILE RAW DATA',/10X 
2,36A2) 

WRITE (6, 6) NPORT.IPT.ITC, MODE, REG 

6 FORMAT (/10X, 'NP0RT=',I3,' IPT=',I3, ’ ITC=',I3,' MODE='»I3, 

7' REG=',I3) 

WRITE (6,8) 

8 FORMAT (/10X, 'NPORT' ,5X» ' IPT( PRESSURE ) ' ,4X» *ITC( THERMOCOUPLE) MOD 
IE' , 14X, 'REGIME '/I OX, ' 1=BOTTOM 1 =COMBIST MICRO 0=NO TEMP PROFILE 
5 0=IMTEGRAL, + ,XVO' 

2,2X» ' 0=TURBULENT'/1 OX, '4=CENTER 2=TRANSDUCER 1=GOOSC NECK' 
3,7X, ' 1 =INTEGRAL ONLY' ,3X, ' 1= LAMINAR '/I OX, *7=T0P' ,22X, '2=ELBOW , 
4/37X, '3=BLACKWELL PROBE' ) 

WRITE (6,7)RAD 

7 FORMAT (/I OX, 'INVERSE RADIUS OF CURVATURE ' ,F6.4 ,'( 1/INCHES) ' ) 

WRITE (6,12) 

1? FORMAT (/I OX, 'VELOCITY PROFILE RAW DATA') 

WRITE (6,15) TRECOV,PSTAT,PAMB,RH,X 
15 FORMAT 1/10X, 'TRECOV(MV)=' ,F6.3,' PSTAT(2N.H20)=' ,F6.3, 

1 ' PAMB( IN.HG,32F ) = ' ,F6 ,2 » ' RH(PCT) = ' ,F4.0, ' X( IN, )= ' ,F7.2 ) 

IF (IPT.EQ.1 ) WRITE (6,20) 

IF (IPT.EQ.2) WRITE (6,21 ) 

21 FORMAT ( 1 0X, • Y( IN ).,U( INCHES OF WATER) «') 

23 FORMAT ( 1 0X, ' Y( IN . ) , U( INCHES OF COMBIST FLUID-S.G,=.8154 *') 

WRITE (6,25) ( Y( I ) »U( I ) , I=1,NMAX) 

25 FORMAT ( 9X,F6.4,F8.3,4X,F6.4,F8.3,4X,F6.4,F8.3,4X,F6.4,F8.3,4X 
1 ,F6.4,F8.3) 

IF (ITC.NE.O) WRITE (6,30) 

30 FORMAT (/I OX, ' TEMPERATURE PROFILE RAW DATA') 

IF (ITC.NE.O) WRITE (6,35) PAM, PSTA, PON, RL, TO 
35 FORMAT (/10X, 'PAM(IN.HG,32F)=' ,F6.2, ' PST-WALL( 1N.H20)=' ,F6.3, 

1 ' PDN( IN.H20 )= ' >F6.3, ' RL(PCT)=' ,F4.0, ' TWALL(MV)=' ,F6.3) 

IF (ITC.NE.O) WRITE (6,40) 

40 FORMAT ( 1 0X, *2( IN. ) , T(MV)« ') 

IF (ITC.NE.O) WRITE (6,25) (Z(I),T(I), I=1,NZMAX) 
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ORIGINAL PAGL Its 
OF P(X)R QUALITY 



C 


61 



iniPT.ru. 11 CAl 10-0.8 156 

62 



TRITOV TCURtCOV) 

63 

C 


iniPT.ru. 2) CAIXP"! .0 

ri)YH DYNAMIC PRESSURE IN INCHES OF HCO 

66 



rDTN UlNMAXJKCALIP 

65 



CAT l TlNNl L 1 PDYN , RAD , YC NMAX ) , PSTAT , PAM8 , TRECQV , RH , UINF , TINF , VI5C , C 
1PiRIIPi?,RM*P»PR) 

66 

C 


UPN UINF*l l.O*RA0*Y(NMAX)) 


C 


CORRECTION TO MICROMETER READING 


e 


Al'Pl Y YOUNG C MAAS CORRECTION FOR PROOE 


c 


Oil UP IWHGR 


c 


FIN!) DYNAMIC PRESSURE 

67 

c 


DO 120 11, NMAX 

INI ERPOl ATE TO FIND TOTAL PRESSURE 

68 



A UPADHrm 

69 

c 

120 

UtlKUiD POYN#M ,0»1 .0/CAMA)) 

FIND GRADIENT Of DYNAMIC PRESSURE, 


c 


YOUNT? AND MAAS PARAMETER 

70 



UGRADt 1 ) Q121/Y12) 

71 



HlfNIIAN 1 

72 



DO 123 I 2 ,Nt 

73 


123 

u?RAom iuti+1 j-qn-i n/tYd+i i-yci-i n 

76 



UGRAOl NMAX )-QGR ADC N1 ) 

75 



DO 123 1-1, NMAX 

76 



YMlTI'-ODHRGRADlD/qm 

77 



I? 1 Y!H I l.GT. 0. 1 ) GO TO 126 

7ft 



DELTA' 00*1 .80HYM1I) 

79 



CO TO 126 

eo 


126 

DILYA- 0.1 6*00 

01 


126 

TUU) TINT 

02 


120 

YCI) YCI HDELTA+OFFSET 

03 



II’ (ITC.Nr.O) GO TO 896 

86 



DO 692 1-1, NMAX 

85 



711) 0. 

06 


092 

TUI I )"TINF 

07 



N’MAX-0 

08 



GO TO 312 

09 


096 

IF mC.EU. 1 ) OFFS~0.006S 

90 



IF CITC.ru. 2) OEFS' ; O.OC15 

91 



IF (ITC .EQ. 3) OFFS". 0055 

92 



DO 360 I' 1 , N/MAX 

93 




96 



ir ; tNPORT.ru. n ZAm=zm 

95 


330 

TCI)-TtCtTlI)) 

96 



TG'TCC TO ) 

97 



TIN11 ■ TlNZMAX) 

90 



IF INPORT. FU.1) ZLIM-NJMAX 

99 



IF CMODl.tU. 2) GO TO 312 

100 



IF t YC i l.GE.Zl 1 ) ) GO TO 301 

101 



DO 309 I'l.NiiAX 

102 



NT --I 

103 



IF CYm.GC.ZU )) GO TO 300 

10* 


309 

TUC I l~TO+CT( 1 )»TO)«YCI)/ZCn 

105 


300 

KCNT 

106 



GO TO 302 

107 


301 

KC"t 

too 


302 

NC"1 

109 



DO 303 X KC.NNAX 

110 



IF (Y(I).GE.ZCNZMAX)) GO TO 307 

in 



DO 306 K"NCiN/MAX 

112 



If C.NOT.Z(K).GT.Y(D) GO TO 306 

113 



J"K =1 

116 



TUm = CTCJt1 )-Tt J ) )»( Y( I )-Z( J) )/CZC J*1 )-ZCJ))+TCJ) 

115 



JJ=K 

116 



GO TO 308 

117 


306 

CONTINUE 

118 


30/ 

TU(I)-TINZMAX) 

119 



GO TO 303 

120 


300 

NC-JJ 

121 


303 

CONTINUr 



122 312 CONTINUE 

C X LOCATION UNCERTAINTY, INCHES 

123 DX : 0 . I 

C TFMPERATURE UNCERTAINTY, DEG F 

124 DT-0.12 

C DYNAMIC, STATIC PRESSURE UNCERTAINTY, SLANT TUBE MAW., IN. H20 

125 DPMAN'0.010 

C DYNAMIC PRESSURE UNCERTAINTY, PM-5,97 TRANSDUCERS, PCT/100. 

126 DPTR"0.5/1 00 , 

C AMBIENT PRESSURE UNCERTAINTY, IN. HG 

127 DPANCT0.10 

1 CB GC-32 . 1 739 

129 FUNIT3=C1 16.21/33.932/12. 

130 C=2.#Ge#PUNrrSHRM*CAlIB/P 

131 DO 275 1=1, NMAX 

C VELOCITY UNCERTAINTY, RATIO 

132 DURm=0.5*5qRT({OT/fTINF+459.69>)**2+DPTR**2+(DPAMB/PAMB)**2) 

133 U( I )=C*( TU( I J+459.69 )*U( I )-UPW*UPW*RAD*Y( I )*( 2.0+RADMY1 I ) )/( ( 1 .0+R 
1AD*Y(I))**2, ) 

134 Uf I )~S0RTf U( I) ) 

135 2 7 5 CONTINUE 

136 DUINF R-DUR ( NMAX ) 

137 WRITE (6,45) 

130 f»5 FORMAT (/10X,'Y LOS (IN.), VELOCITY UNCERTAINTY RATIO, PCT »•) 

139 WRITE (6,25) ( Y( I) ,DUR(I )*1 00. , Is 1 , NMAX) 

140 CFfist . 

141 IF (MODE.EQ.O) CALL STRESS (NMAX,REP'.U,Y,UPW,VISC,CF2,DUR,YPLO, 
2YPHI.CF) 

(42 890 CALL DEL99 ( 1 ,NMAX,U,UPW,RAD,T,Y,TO,DEL) 

193 CALL DELI 2 (NMAX, RAO, UPW.U.YrTU, DEI ,DE2,DUR,DDE2R) 

166 H=OE1/DE2 

U*5 REM=UPN»DE2/(VI8CK|2. ) 

146 DRr,MR^SQRT(DUINFR**2+DDE2R**2)*100. 

167 XVO=0. 

1-,8 DXVO-O, 

U»9 NEX=1 . 

150 DREXR-0. 

151 DEH2=0. 

152 0DEH2-0. 

153 IF (MODE.NE.O) GO TO 310 

154 IF(REG.EO. 1 ) GO TO 891 

1 55 XVO=DE2*( ( UPW/VISC )**0 . 2 )/( 1 2 . #0 . 037 ) 

156 XVQ"(XVO**l ,E5)*12. 

157 DXVO=XVO*SqRTm.C5*DDE2R)**2+(0.25*DUINFR)**2> 

158 GO TO 093 

159 891 XVO"(DF,2*DE2*UPW)/(5.292*VISC) 

160 B\VO=XVO*SQRT( ( 2 . 0*DDE2R )**2+( OUINFR )**2 ) 

161 893 DXVOR- DXVO/XVO 

162 XVO-X-XVO 

163 IF (REG.EQ.O) REXs(REM/0.037)**t .25 

164 IF (REG.EQ.1) REX=(REM/0.664)**2 

165 DREXR"SQRT( DUIhFR**2+( DX/XVO )**G+( DXVOR )**2 )*1 00. 

166 310 WRITE (6,305) R EX » DREXR , REM » DREMR , XVO , DXVO 

167 305 rORMAT (/I0X, 'REX=' ,E12. 5, 5X, 'UNCERTAINTY IN REX (PCT) =',F6.3,//, 

110X, 'REMS', F6.0.5XV UNCERTAINTY IN REM (PCT) =' ,F6.3,//10X, 

2 'XVO (IN.) =',F6,1 ,5X, 'UNCERTAINTY IN XVO (IN.) =',F6,3) 

168 IF (ITC.EQ.O) GO TO 600 
C 

C ***** TEMPERATURE REDUCTION SECTION ***** 

C 

169 UISC=VISC 

170 ouinfrsdphan/pon 

1 7! CALL DEL99 ( 2,NZMAX,U,UPW,RAD ,T,Z,TO,DELT99) 

C NO PROVISION IS MADE TO MAKE SLIGHT ADJUSTMENT IN VELOCITY 

C PROFILE IF PDN.NE.POYN 

1 72 PDM=PDYN 

173 CALL TUNNEL ( PDN.RAD ,Z(NZMAX) ,PSTA,PAM,TINFF,RL,VINF,TINF,VIST,CP, 
1RHOG,RM,P,PR) 

174 VPW=VINF*(1 ,0+RAD*Z(NZMAX)) 

175 IF ( Z(l).GE.Yd)) GO TO 401 

176 DO 399 1=1 ,NZHAX 

177 NT=I 
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17 * 
I 79 
ISO 
181 
I St 
IS) 
184 
16 $ 
166 

187 

188 
169 

190 

191 

192 
195 

194 

195 

196 

197 

198 

199 


200 

201 

202 

20 ) 

204 

205 

206 
207 
206 
209 


210 

211 

212 

21 ) 

214 

215 

216 

217 

218 

219 

220 
221 
222 
22 ) 

224 

225 

226 


227 


226 

229 

2)0 

231 


IF I2IXI. 62. VI 1)1 60 TO 400 
VI I HUM >*ZII)/TII) 

)99 CONTINUE 

400 KC*NT 

CO TO 402 

401 KC*1 

402 NC*1 

DO 40) I*KC.NZMAX 
IF IZIII.6E.VINMAX)) 60 TO 407 
DO 404 K-HC.NHAX 

IF I .NOT. TIKI. GT. 21 II) GO TO 404 
J*K-1 

VI X t*UI JHIUI JM l-UIJ >1*121 I >-T I J > )/< VI J» I 1-TI Jll 
JJ*K 

00 TO 408 
404 CONTINUE 

407 VI I HUPU/I 1 . 0*RA0*2( 1 1 I 
DURI X HDUINFR 
60 TO 40) 

401 NC* JJ 

IF I NZMAX. GT . NMAX ) DURI I I *DUINFR 
40) CONTINUE 

C OBTAIN STATIC TEMPERATURE PR0FUE5 
C RECOVERY FACTOR FOR TC HIRE NORMAL TO FIOM 
RTC*0 .66 

00 987 1*1, NZMAX 

967 T I IHTII)-RTC*VII)*VII)/I2. *32. 2*776. 26*CP) 

DELT*TO-TINF 

CALL DELH2 I NZMAX ,RA0 ,UPN . V.T ,2 , TO.CP.DEM2 ,DUR .OSUMMR ) 

REM=VPW*DEH C/VIST/I 2 . 

DREMRtSQRTI DUINFR**2»D5UHHR**2 )*1 00 . 

00CM2*DEH2*05UMHR 
WRITE 16,595) REH.DREHR 

595 FORMAT I /I OX , 'REM* *, F6 . 0 , 5X , 'UNCERTAINTY IN REM I PCT 1 ■',F4.I) 

C 

c 

c 

600 WRITE 16,605) INRNII), 1=1,36) 

605 FORMAT! 1H1 ./9X.56A2 ) 

XVOCM s XVO*2 . 54 
DE2CM=DE2*2.54 
0EM2CM=DEH2*2.54 
UPWMS*UPW*0.)04« 

DELCM=0EL*2 . 54 
9E ITCM S 0E LT99*2 . 54 
VISC5I*VISC*0.0929 
RHOGM=RHOG* 16.02 
DE ICM*DE 1*2.54 
VI’WMS s VPW»0 . 3046 
VISTSI=VI5T*0.0929 
XCM=X*2 . 54 

TINFC*5 . *1 TINF-)2 . )/9. 

TOC*5 . *1 TO-)2 . )/9. 

IF (ITC.NI.O) WRITE 16.610) REX ,REM,REH .XVOCM .DE2CM.DEM2CM, 

1 UPWM5 ,DE LCMiDELTCM, VISCSI .DEI CM , VPWM3 >NPORT ,H .VIST SI , XCM , 

2 CF2.T1NFC.RMOGM.TOC 

610 FORMAT I/10X, 'REX * ' , E 1 2 . 5, 7X , 'REM »' .F12.0.24X, 'REH *',F12.0 

l/IOX.'XVO * ' ,F 1 2 . 2 , ' CM DELE *',F12.),' CM* , 21 X, '0EH2 «',F12 

2.),' CM 1 ,/1 OX, 'UPW *' ,F12.2, ' M/S DEL99* ‘ ,F 1 2 . ), ' CM',21X,'DEL 
1T99 = ' , F 1 2 . 5 , ' CM ' , /I OX , 'VI5C s'.EIZ.S,' M2/S DELI *',F12.3,' CM 
4 ' , 2 IX, ‘ UPW *',F12.2,' M/S ',/10X,'PORT *' ,9X,I3,7X, 'H *',F12 
5. 3.24X, 'VISC s'.EIE.S,' M2/S './lOX.'XLOC =',F12.2,' CM CF/2 
6 s ' ,E12.5,24X, ' TINF *•,712.2,' DEG C ' ,/1 OX, 'DENS «',F12.2,' KG/M) 
7' ,43X, 'TPLATE *',F12.2,' DEG C' ) 

IF IITC.NE.O .ANO. STANT . NE . 0 . 0 ) WRITE 16,611) STANT 

611 FORMAT l 7 7X , * STANTON^ ' , E 1 2 . 5 > 

IF (ITC.EQ.O) WRITE 16,615) REX, REM, XVOCM, 0E2CM.UPWMS.DELCM, 

1 VISCSI, DE1CM.NPORT.H, XCM, CF2.RHOGM 

615 FORMAT l/10X,'REX * ' , E 1 2 . 5 , 7X, 'REM *',F12.C, //10 

IX.'XVO *' ,712.2, ' CM DELE »',71*. ),' CM 

2 ZIOX.'UPW *',712.2,' M/S DEL99*' ,712.5, ' CM 

3 ZIOX.'VISC *' ,E12.5, ' M2/S DELI «', 712.3,' CM ', 
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I M 
2)) 


239 

2W 


2 56 

237 

210 

21 » 

ZhO 

*41 

2*3 

29 * 

2«S 

296 

297 

?<*a 

299 

zso 

251 

252 

253 
259 
255 
25 * 

257 
250 

259 

260 
26t 
262 
263 
269 

265 

266 
267 
260 
269 
2 70 

271 

272 

273 

279 

275 

276 

277 

2 70 
279 

200 

201 

202 

203 

209 

205 

206 

207 

200 


9 /I OX f • FO0T »* .9X.I3.7X. M ■ , ,M2.3, 

5 ZIOX.'XLOC ■ ' »F12.2i * C H CF/2 i'iIIMi 

6/I0X.0EN3 •Sni.ti' KG. M3') 

17 IITC.N2.0 .AND. STANT. EQ.C.C) WRITE 16.620) 

620 FORMAT I/IOX.'YICM) Y/DEL Ul M/S > U/UP Y ♦ U»'.6X. 

l'CF/2 TlOtC C) YICH) Y« UtM/S) TIDC6 C) TBAR I -TBAR 

r i 

IF IITC.Nt.0 . AMO . STANT . NE . 0 . 0 ) WRITE 16.621) 

621 FORMAT l/IOX.'VICM) Y/DEL UIFt'S) U/UP U U**.6X, 

l'CF/2 TIOIO C) YICM) V* Ul M/5 1 TlOlG C) TBAR T* 

1 PRT ' ) 

IF IITC.EQ.O) WRITE 16.625) 

625 FORMAT! /I OX. ‘YICM) Y/OEL UIM/3) U/UP Y* U« CF/2'/) 

F1«UPW*5QRTICF2)/VI3C/I2. 

F 2 ■ 1 ./SQRTICF2 )/UPW 
MIN=NMAX 

IF IN2MAX.LT.NMAX) MII1-N/MAX 
MAX* NMAX 

IF I N2MAX . 6T . NMAX ) MAX*NZMAX 
KK*0 

DO 660 X*1 .MAX 
KK*KK*I 

IF 1KK.NE.6) GO TO 602 
KK*1 

WRITE 16.601 ) 

601 FORMAT I 10X) 

602 CONTINUE 
YCH=Y(I>*2.59 
ZCM-ZI I )*2 . 59 
UMS>U( 1)90. 3090 
UM3T *VI I )*0 . 3090 
TK=5.*ITI X )-32. )/9. 

YBAR* V( X )/DEL 

UBAR >U( I )•( 1 . 0*RA0*Y( I ) l/UPM 

YPL»Y(I)*F1 

ZPL*Z( X )*F 1 

UPL«UC X >*F2 

IF IITC.EQ.O) WRITE 16,6951 YCM. YBAR ,UMS ,UBAR , YPL ,UPL ,CFI X ) 

II l XTC.EQ.O) GO TO 660 

TBAR*ITO-TlX))/DElT 

IF I3TANT.EQ.0.0) GO TO 603 

IF IX .GT. NZMAX) GO TO 603 

TPLUS s l 1 SQRTI CF2 ) l/STANT >*TBAR 

PRT*I I TPLU3-1 3.2*PR 1*0.91 )/l A LOG! ZPL/1 3.2)) 

603 CONTINUE 
TM1 *1 .-TBAR 

IF IX. LE. NMAX) TXK«5.*I TUI XT-32. )/*. 

IF I I .GT MIN) GO TO 690 

IF (STANT. EQ. 0.0) WRITE (6.630) YCM, YBAR ,UMS .UBAR , YPL .UPL.CF ( I ) , TK 
2K. ZCM. ZPL. UMST, TK, TBAR, TM1 

IF (STANT. NE .0.0) WRITE 16,630) YCM, YBAR, UNS. UBAR, YPL. UPL, CFI I ),TK 
2K . ZCM , ZPL .UMST . TK , TBAR , TPLUS , PRT 

630 FORMAT I 0X . F 7. 3,F7. 3 , F7. 2 , F0 . 3. F7 . 1 , F0 . 2 , F 1 0 .6 . F0. 2 , SX , F0 . 9 , F0 . 2 , 
2F7.2>F9.2,2F7.3,1F9.3) 

GO TO 660 

690 IF I NMAX. GT. NZMAX) WRITE (6.695) YCM . YBAR ,UMS .UBAR . YPL , UPL , CFI X ) ,T 

2KK 

695 FORMAT I 0X , F 7 . 3 . F7. 3 , F 7. 2 . F0 . 3 ,F7 . 1 . F0 . 2 . F 1 0 .6 , F0 . 2 ) 

IF ( NZMAX. GT. NMAX .AND. STANT . EQ. 0 . 0 ) WRITE 16,650) ZCM, ZPL , UMST ,T 
1K.TBAR.TM1 

IF I NZMAX . GT . NMAX .AND. STANT . NE . 0 . 0 ) WRITE (6.650) ZCM, ZPL UMST.T 
1 K , TBAR .TPLUS , PRT 

650 FORMAT ( 79X, F9. 9 , F0 . 2 , F 7 . 2 . F 9. 2 ,2F 7. 3 . 1 F9. 3 ) 

660 CONTINUE 
GO TO I 
670 CONTINUE 

WRITE (6.600) 

600 FORMAT I INI ) 

RETURN 

END 
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269 

C 

FUNCTION TCIT) 

FUNCTION CONVERTS TEMP FROM IRON-CONSTANT AN 


C 

FREE -STREAM THERMOCOUPLE CALIBRATED 3-60 

290 


TC«S2 . 350*35. 300»T- .905»T*T 

291 


RETURN 

292 


END 


291 FUNCTION T7CITI 

CC THIS FUNCTION CONVERTS TIMP FROM CHROME L -CONSTANT AN MV TO DCS f. 
C BLACKWELL THCRflOCOUPLC PROBE CALIBRATION 3-80 
299 TTC»32. 976*31 .026*T-.59000»T»T 

295 RITURN 

299 END 


297 


296 

299 

300 

301 

302 

303 
309 

305 

306 

307 
306 

309 

310 

31 1 

312 

313 
319 

315 

316 

317 
316 

319 

320 

321 

322 

323 
329 

325 

326 

327 
326 
329 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 


SUBROUTINE TUFMEL I POYN.RAD , YMAX . PSTAT , PAMB . TRECOV .RHUM.UINF , TINF , 
1 VISC (CP iRMOS »RH . P . PM I 

RHOG FREE STREAM DENSITY 

VISC FREE STREAM KINEMATIC VISCOSITY 

CP FREE STREAM SPECIFIC HEAT 

PR FREE STREAM PRANOTl NUMBER 

H FREE STREAM ABSOLUTE HUMIOITY 


THIS ROUTINE COMPUTES THE HIND TUNNEL FLOW CONDITIONS 

UINF FREE STREAM VELOCITY 

TINF FREE STREAM STATIC TEMPERATURE 


SATURATION DATA FROM K AND K 1969 STEAM TABLES 
DIMENSION TEMPI 1 0 1 , PSATI 1 0 1 , RHOSAT! 1 0 1 


OATA TEMP/ 
90.0. 
DATA PSAT/ 


1 


I 


100.627. 
DATA RHOSAT/ 


I 


90.. 

100 . 0 . 

17.519, 

136.693. 

.0009090, . 0005666 . 
.0026571, .0037722. 


50.0. 

110 . 0 , 
25.636, 
163.767, 


60.0, 

120 . 0 , 

36.907, 

299.006, 


70.0, 

1 30 . 0 / 
52.301 . 
320.900/ 


60.0, 


73.051 , 


.0021361 

REAL NU.MFA.MFV.MMA.MUV.JF 
TAMB'TRECOV 
DO 10 N=1 ,9 

IF! TEMPI NI.GT.TAMBt GO TO 20 
10 CONTINUE 
20 T * TEMPIN) 

EPS * T - TAMB 
VAPM 
VAPL 
VEPS 
RHOH 
RHOL 
REPS 
RHOG 

RAHS9S. 32/26. 970 

PG « VAPL ♦ (10.0 - EPS l»VEPS/1 0 . 0 

PUNITS s 21 16.21/33. 932/12 . 

P*PAMD*2 1 16. 21/2 9. 92 13 ♦ PSTAT*PUNITS 
RHUM=RHUM/100. 


.0006266, .0011525, .0015603, 
.0099261, .0063625/ 


PSATI N) 

PSATI N- I » 

VAPM - VAPL 
RHOSATI N I 
RHOSAT I N- 1 1 
RHOH - RHOL 

RHOL ♦ (10.0 - EPS )«REPS/1 0. 


PVAP * RHUM-PG 
PA « P - PVAP 

RHOA * PA/I RA*( TAMO ♦ 959.671) 
RHOV * RHUM*RHOG 
W=RMOV/RHOA 
RHOM * RHOA ♦ RHOV 


MWA * 26.970 
MHV * 18.016 
MFV » RHOV/RHOM 
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jso nr a • i.o - nrv 

331 bm ■ i545.32"<mfa/mwa ♦ mfv.mwv> 

33? cp ■ mfa«o.? 40 ♦ nrv»o.445 

333 GC*32.1739 

334 JM778.26 

333 RCF«0.7"*0.S3SS3 

3)4 RHOG=( P/RM*POYN*PUNITS*RCF/1 CP*JP 1 1/1 TRECOV*459.47 1 

337 UINFiSGRTII.-GCTOYN-PUNITS/RHOG)/! 1 . 0»RAD«YNAX ) 

C OBTAIN STATIC TEMPERATURE 

C USING UINF FOR RECOVERY CORRECTION 

C UINF IS APPROX FREESTREAn VELOCITY AT FREESTREAM THERMOCOUPLE 

C RECOVERY FACTOR FOR TC HIRE PARALLEL TO FLOW 

338 RTC*0.66 

339 TINF *TRECOV-RTC*UlNF*UINF/( 2 . •GC*JF*CP ) 

340 VISC*I 11 .*0.0175*TINF )/! 1 . E06"RHOG )•( 1 .-.7«»W) 

341 PR* . 710*1 5)0. /( TINF *45 9. 47 I )••! . 1 >•( 1 . ♦ . 9«W > 

C NOTE FOR HIGH VELOCITY THIS ROUTINE SHOULD BE ITERATED 

342 K| TURN 

343 END 


344 SUBROUTINE STRESS) NMAX .REG, U.Y , UINF ,VISC ,CF2 .DJR , VPLO. YPHI , PSUMCF ) 
C 

C FRICTION COEFFICIENT FOUND BY SEARCHING DATA 

C FOR ALL POINTS SATISFYING FLAT PLATE LAW-OF-THE-WALL 

C BETWEEN THE RANGES OF YPLUS OF 75 AND 125 

C FOR THE CASE OF LOW RE USE THE INTERVAL OF YP BETWEEN 22 AND 

C 62. FOR CURVED AND RECOVERY REGIONS USE THE INTERVAL OF YP 

C BETWEEN 40 ANO 90 

C AN AVERAGE CF2 IS OBTAINEO FOR T HESE YPLUS VALUES 

C FOR LAMINAR BOUNDARY LAYER USE THE INTERVAL OF YP BETWEEN 0 AND 

C 20. 

C 

345 DIMENSION U( 11 , YM > .DUR! 11 .PSUMCF ( 60 1 , YPLUS! 60 ) , YLISTI 60 ) 

346 INTEGER REG 
C 

C • • • » a * NOTE' HIGH RE FOR UINF« 55-112 FT/SEC 

C • * • • * • NOTE' LOW RE FOR UINF .IT. 50 FT/SEC 

C 

347 3UMCF*0. 

348 NUMCF=0 

349 F YPO*0 . 

350 0CFR*0. 

351 EPS=5.E-3 

352 J*0 

353 YP*YPLO 

354 DUINFR=DUR(NMAX) 

C UNCERTAINTY IN Y POSITION. IN. 

355 DY*0 . 002 

356 IF (REG.EQ.O ) GO TO 602 

357 00 601 1=1 ,NMAX 

358 PSUMCF 1 1 )*1 2 . 0"UI I >"VI5C/« UINF*UINF»YI 1 1 1 

359 601 YPLUS! I )=S9RTI PSUMCF! I ) )»UINF •Y(I)/(VISC*12.0) 

360 602 DO 800 1*1 .NMAX 

361 PUY«U(I)*Y(I)/( 12. »VISC I 

C UNCERTAINTY IN Yt, RATIO 

362 DYPR*0 . 1 0 

363 IF (REG .NE. 0 I GO TO 603 

364 KONT *0 

365 600 KONT*KONT*1 

366 FYP*YP»! ALOG! YP1/0.4U5. 1-PUY 

367 IF 1YP.GT.YPMI) EPS«5.2-2 

368 IF ( ABS! FYP-FYP01.LT.EPS! GO TO 700 

369 EPS=5.E-3 

370 FOERV*! A LOG! YP ) ♦ 1 . I/O .41 *5. 

371 YP*YP-FYP/FOERV 

372 IF I KONT. LE. 200) GO TO 600 

373 WRITE (6,910) Y( I > 

374 910 FORMAT MOX, 'FAILURE TO CONVERGE FOR CF2 ' , 1 OX, * Y* ' ,F8. 3 ) 
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)7» RETURN 

37* 700 CONTINUE 

377 VPLUSIIUYR 

37* PSUNCFI I >■! VISC«YP«17./I Yl I H*UINF 1 )••* 

370 000 CONTINUE 

300 003 CONTINUE 

301 DO 075 IM.NNAX 

307 IF ( Y PI US I 1 1 . IT . TPIO I 00 TO 075 

303 IF IYPLUSI XI.OT.YPMI ) 00 TO 050 

304 SUNCF »SUNCF «PSUNCF I I ) 

305 DCFR*50RTI DCFR*»«7*4.6I DY/YI X I •••7*4.»0UINFR»«?*4 . •OYPR»«f 1 

304 NUNCF *NUNCF*1 

307 075 CONTINUE 

300 050 CF7»SUNCF/FU)ATI NUNCF ) 

304 DCF R*DCFR/F 10ATI NUNCF 1*100. 

3«0 WRITE 16.400) CF7.DCFR 

341 400 FORNAT I /I OX . * CF/7 * ' . E 1 7 . 5 . OX . 'UNCERTAINTY IN CF/E IPCTI ■ , ,F*.7> 

347 RETURN 

343 ENO 


344 SUBROUTINE DEL44 ( N.K1 .U.UPW.RAO ,T , Y ,TPL ,0E 1 
C SUBROUTINE WRITTEN BY NARCOS PINENTA 

345 DIMENSION Ul 1 I. Til 1.Y1 1 1.7(601 

344 ZBt.44 

347 00 101 Ml . K 1 

340 IFIM.EQ. 1 IZIKI'UIK >•( I .0«RAD*Y(K1 l/UPW 

344 IF(M.EQ.7)Z(K)*( TPl-TI K 1 I/I TPI-TCK1 I) 

400 101 CONTINUE 

401 00 10 Ml .K1 

407 IFI7IKI.GE.7RIG0 TO 7 

403 10 CONTINUE 

404 7 0 = Y(K-7 )«*2*Y(K-1 >♦ Yl K 1**7*Y( K-7 K-1 )»»C»Y( K l-YI K )»»2*Yt K- 1 »-Y( 

IK-1 >««»?«YC K-C )-Y(K-2 (*»7« Y t K ) 

405 F 1 *( I 7( K-1 1-71 K-7) MM YIK1-Y1 K-7 11-171 K 1-71 K-7 ))•( Y( K-1 I-YIK-7I11/D 

406 FC*( I7(K)-7!K-7))*(Y(K-1 1»*7-Y( K-7 >**7 »-( 7( K-1 )-7( K-7 ))•( Y(K >*«7-Y 
UK-7 )*»7 ) I/O 

407 F3*l YIK-7 1**7*Y(K-1 1*71 K >♦ Yl K 1»*2*Y( K-7 )*7( K-1 UYIK-1 )»»J*Y< K )»7( K 
1-71-Y1K)»*7»Y(K-1 )*ZI K-7 1-Y( K-1 >*»7*Y( K-7 l*Z< K >-Y( K-7 l**7*Y( K )*7( K 

7-1 ) I/O 

400 OE«VIK-l) 

404 DO 3 1^1 .10 

410 OE^DE*! ZB-F1 •DE«»7-Ff »0E-F3 i/I 7 .*F 1 •0E*F7 ) 

41 1 ERR = ( 7B-F 1 •DE»*2-F2»0E-F3 )/1 2 .*F 1 •0E*F7 1 

417 IF1ERR.LE.1 .E-03IGO TO 4 

413 3 CONTINUE 

414 4 CONTINUE 

415 RETURN 

416 ENO 
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417 

416 

416 


410 

421 

422 

423 

424 
426 

426 

427 
426 
424 

430 

431 

432 

433 

434 

435 

436 

437 
436 

439 

440 

441 

442 

443 

444 

445 

446 

447 
446 

4*9 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 


SUBROUTINE Dill 2 (Kt ,RAO.UPM,U.Y,T,Dt1T,DE2T,DUff.DOE2R> 
DIMENSION Ul 11. VM 1.TI 1 t.DURt I 1.7(601.61601 
DIMENSION 77160). 6GI60) 

C DUR. UNCERTAINTY IN VELOCITY. RATIO 
C DT. UNCERTAINTY IN V LOCATION, IN. 

0Y*0 . 002 
DO 10 1*1, Kl 
7( X )*U( 1 1/1 T1 1 ) *454. 64 1 
77 1 1 1*1 . -U( I l/UPM 
G6( X l*U( I )*E7( I l/UPW 
10 G(I)*7(X)*U(X) 

701-71 1 )*0.5 

G01-6I 1 )*0.5 

DY01-YI 1 1 

SUM1»701*DY01 

5UMP : 0 T 0 1 /( 1 . 0*RAD*Y( 1 )) 

SUM2-G0 1 *DY01 •( 1 .0«RAD*Y( 1 1 1 
DSUMUI 0.5*771 1 )*OT)**2 
05UM2*(0.5*GG( 1 )*DYI**2 
KLIMsKI-1 
DO 20 1*2. KLIM 
OY12«Y(I)-Y(I-1 1 
712*0. 5*(7( I I. 7(1-1 I) 

G12«0.5*(G( I )«CI X-1 )) 

SUM1 *5UM1 *71 2*DT 1 7 
YM*0 .5*1 Y( 1 1 ♦ Y 1 1-1 It 
SUMP*5UMP*DY1 2/( 1 .0*RAD*YH) 

sum:*sum:*gi2*dyi2*( i .o*rad*ymi 

D7 1 *1 ,414*DURI I )*U( I )/UPM 
0G1 *DF 1 •( 1 . -2 . *U( I )/UPW 1 

DSUM1 =DSUM1 ♦( 0 . 5*( Y( I ♦ 1 I— Y ( 1—1 1 )*D7 1 )**2 ♦ (0.5*177(1-11- 
1 7 7 ( 1*1 1 )*DT)**2 

DSUM2=DSUM2*(0.5*( Y( 1*1 l-YII-1 I l*DG1 )**2 ♦ ( 0 . 5*( G6( 1-1 1- 
1 GG( I * 1 1 )*0T 1**2 
27 DY01 -DY 1 2 

DYN*Y(K1 1 — YE K 1 -1 I 

5UM1 *SUM1 *0 . 5*( 7( XI )*7(K1-I 1 1*DYN 

YM*0 . 5*1 Y( K1 I ♦ Y( K 1 — 1 II 

SUMP: SUMP *0YN/( 1 . 0*RAD*YM ) 

SUM2*SU112*0.5*(G(K1 )*G(K1-1 11*0 YN*( 1 . 0*RAD*YM I 
D7N* 1 . 41 4*0UR1 K1 I 
DGN=D7N*( 1.-2. I 

DSUM1 =DSUM1 ♦! 0.5*DYN*D7N1**2 ♦ ( 0 . 5*( 77 < K1 ) *77 ( K1 -1 ) )*0Y )**2 
05UM2 : 0SUM2*( 0. 5*DYN*DGN 1**2 ♦ ( 0.5*1 GGIK1 )«GG(K1 -1 I )*OY )**2 
DSUM1 *SQRT ( 03UM1 1 
DSUM2*SQRT ( DSSJM2 ) 

DE 1 T*SUMP-SUM1 •( T ( K 1 )*459.69 >/UPW 
DE2T*SUMP-DE 1 T-SUM2*( TIKI 1*459. 69 1/( UPW*UPM I 
DDE2R-DSUM2/0E2T 
WRITE (6,30) OE1T.DSUM1 

33 FORMAT) /I OX, ’DEI (IN.) * ' , 7 7.4 , 5X, 'UNCERT IN DEI (IN.) *’,77.*) 
WRITE (6,40) 0E2T.DSUM2 

40 70RMAT(/10X« 'DE2 (IN.) * * , 7 7. 4 , 5X , 'UNCERT IN DEE (IN.) *',77.41 
RETURN 
ENO 
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470 

471 
47* 

473 

474 
47* 
4 7* 

477 

476 

479 

460 

461 
46* 
46) 
484 
46* 
464 

467 

468 

469 

490 

491 

49* 

49) 

494 

49* 

494 

497 

496 

499 

500 
*01 
* 0 * 
* 0 ) 
*04 


SUBROUTINE OHM* 1 K 1 .640 ,UPN ,U,T, Y , TPL .CPI .SUNK ,OUR .OSUTWRI 
C CAUTION THI* ROUT INI COMPUTES (NTHALPT THICKNESS BA5I0 
C UPON STAGNATION IN T MALPT PROFILES. THIS PRISUMS 

C THAT TH* TEMPERATURE PROFILES ARE STATIC TEMPI RATURE . 

DIMENSION UIII.TI1),ri1 t.OURl 1 I ,ZI 40 1 .021 40 I ,ZZ< 40 I 
C OUR. UNCERTAINTY IN VELOCITY. RATIO 
C OT, UNCERTAINTY IN TEMPERATURE, OEG f 
DTM.ES 

C DT, UNCERTAINTY IN V POSITION. IN. 

OY*0 . 00* 

00 10 IH.KI 

TB*I T( I I-TIK1 II/1TPL-TIK1 II 

0Z( II*UI XI*TB/U<K1 1*1 .414*5QRT(0UR( I )■•*«( OT/( TPL-TI Kl 1 1 )••* I 
Z*tII*UtII4TB/U<K1 I 

*( I »«CP1*< Tl I I-T1K1 I I ♦IlN 1 !••*-( UPM/I 1 .0*RAD*Y( 1 1 1 )••* 1/ 

1 (*.•)*. 1 74*776. I 
10 mi«UIII*2III/(T(II»4*9.49l 

2VI«CP1*(TPL-T<KI l)-UPW***/(*.*32. 174*776. I 
201 *21 1 1*0 . * 

DYOMYI 1 I 
SUMH*201*0YC1 
0SUMMM0.WZ1 1 )*OY !••* 

KLIM*K 1-1 
00 20 1*2 .KLIM 
0 Y 1 2* Y( I I-Y1 1-1 I 

C«(II*(I-I 1-21 1 1 I/I Yl 1-1 l-VIXIII-IIZIXI-ZIXH I I/I Y1 X l-YI !♦! 1111/ 
1IYII-I l-YI 1*1 II 

B*< 1 Z( X l-Z( I* 111/1 VII l-YI IH >ll-C*(Y(IltYIXH II 
A*Zl 1-1 l-B*Y1I-1 ) — C • Y I 1-1 ) • Y ( X-l I 

SUf1MsSUriH»A*lY(X*1 l-YI 1-1 ll/2.«.12S»B*<'V(I«1 I ♦ Yl 1 1 !••*-( Yl 1-1 !♦ 

i u 1 1 i,o. o>,i ;«c»i I VI i*i uti 1 1 >*•)-( yi i-i ,*u i . »»*)i 

D5UMH s 0SUMH*1 0 . 5*1 Y( I ♦ 1 I — V 1 1 — I 1 1 *021 X 1 1**2 ♦ 1 0 . **« 2Z1 I-I l-ZZl 1*1 I 
1 >*0Y 

*o oyoi*oyi2 

OYN*Y(K1 l-YI Kl -1 I 
5UMMiSUMH*(Z«K1 UZ1K1-1 I l*DYN*0.5 

0SUMN:DSUMN»I0.*«0Z(K1 I l«>2 ♦ < 0 . *•< 2*( K 1 )♦**< K 1 - 1 1 1*0 Y >••* 

DSUMM*S<JRT 1 OSUMN I 

SUM»l 5 5Ut1H*l TC Kl >«459 . 49 1/1 UPW*ZM I 

DSUMHRxDSUMM/SUMF 

IF ( RAD .LT. 0.001 I GO TO 2* 

5 UMH* SUMMER AD 

SUMH = ( 1 . 0/ RAD |*( EXP( SUMH 1-1.01 
2* WRITE (4.301 SUMH.OSUMH 

30 FORMAT1/10X, 'OEH2 1 IN. I* * ,F7.4 , *X , • UNCERT IN OEH2 ( IN. I* ' ,F7.4 1 
RETURN 
END 
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Appendix F 

Listing of the Stsnton Pr rgram 

The following listing Is thst of the Fortrsn program used to reduce 
the wsl 1-aessured hest flux dsts. 


F-l 


MATFIV 


STANTON NUMBER DATA REOUCUON PROGRAM 
CURVA'UBE RIB NASA - NAG- 1- I 

THIS PROGRAM CALCULATES STANTON NUMBERS AND (NTNALPT THICK. 
REVISED JULT 1900 

PRESENT VtRSION ALSO CALCULATCS 3EN3ITIVITT COEFFICIENTS TOR 
STANTON NUMBER AND UNCtRTAINTV IN STANTON NUMBER 


1 REAL K 

2 INTEGER PLATiO. PI ATM, 02LOC, PLATE2 

I COMMON ARGI611.ARG0 

A 01 Ml NS I ON SCTI62.A0.ll 

I 01 Ml NS I ON SMAXDRl AO I .SMINDR1 60 I . SMA/TSI AO I «SF*IN T St 601* 

1 SMAXRRt 60 I iSMINRRI 60 I 

6 OIMCNSION ST I 62 1 1 P0< 62 I tUNCERTI SS l .OSTI At ) • 
t USTRE It 62 1 .01 LOCI 62 1 .002 LOCI 62 I 

7 01 Ml NS I ON 01 SC I 6A I . XCMI 62 I 

S COMMON/ BIM /PANS, PST AT ■ TRICOV ,RMUM, POTN. TAMO 

* COMMON/ BLK2 /UINf .TINT .TAOIAB.RMOINF .VXSC.PR.CP.M 

10 COMMON/ HIM /TOO.TOI61I.MMI62I.QI 10I.TENDI Ml. 

1 Tf RAMI I 10l.OOOTt62l.OQOOTl62l.CMt Ml. 

2 OKKI 62 I .OSSt 61 1 .NTACt 61 1 .NTAOO .RCCNPI 62 I »STPt 62 I 

11 COMMON/ BLKA/OSTt 62 I .USTABSt 62 I .RUN! It I .DRCINt 62 1 .021 62 1 .0021 62 I 
1 .PC0C7M62 I.XI62 I 

It COMMON/ BLKS /SINC0I62. lAl.VMAX.VMIN.NA.NO 

II OIMINSION TTOPt JOI.TCINI M ) . TBOTt M I .R2XI 62 I , 

1 NRNI 6 I , KOMiVIT 1601, TSUPPt I 2 I .B( 66 I .MUM I 6l I . UB> 6} I 
16 OIMINSION TTOPMVt lAl.TCfNMVt JAl.TBOTMVt JAl.TSUPMVt 121, 

1 TOMVt 66 I . POST! 62 I . PXCHt 62 I 

15 CALL START(,i'&lNll**,0.) 

16 OATA X/ 25. 17. 26. 20. 27. 21. 20. 25. 26.20, 10. 11, 31. 36, 12. 17, II. 39, 

1 16. 62. 15. 65. 16. 60. 37. 51 .JO. 51. 19.56. 6059, 61 .62.62.66, 

2 ■* 1.67,66. 70,65.71,66.76,67.70,60.01 ,50.25,52.10,56.10, 

1 56.17,50. 17. 60. 16. 62. 16, 66. 15, 66. 15. 60. 15, 70. 16, 72. 16, 

6 76. IS, 76. 11,77. .,6, 70. 67, 79. 70, 00. 71, 01 .75,02.70,01.01 , 

5 06.06,05.06,06.09,07.92,00.95.09.90,91 .00,92.01,91.06, 

6 96 0 9,95 1 2, 96. 16,97. 1 7.90. 20,99. 21, 100. 26, 101 .20/ 

c 

17 00*0.01 


OATA 

6/ 

.0157, 

.0121, 

.0111, 

.0065, 

.0119, 

1 

.0120, 

.0167, 

.0166, 

.0000, 

.0097, 

.0105. 

1 

.0111, 

.0107, 







.00909, 

.01026, 

.00962, 

.01510, 

.01020, 


.00505, 

.00907, 

.00706, 

.00000, 

.00617, 

. 00 M2 • 


1 *»»0 

0, 






.00752. 

.00066, 

.00631, 

.01165, 

.00010, 

.02906, 


.01560, 

.00061, 

.01210, 

.00971 , 

.01562, 

.01210, 


.00001, 

.01207, 

.00010, 

.01077, 

. 01 M6, 

.01067, 


.01211, 

.00011, 

.00700, 

.0(655, 

.01579, 

.00776, 


19 IRUN*0 
C 

C CHANGE 06IGIN 07 X FROM START OF COPPER PLATES 

C TO START Or CURVATURE 

C 

20 UBO*. 05*00 

21 CO 110 1*1,62 

22 UOI XI*. 05*Bt 1 1 

21 110 XtXI * XIII - 69 . ]125 

26 USt 61 1*. 05*0161 1 

C 

25 00 120 1*1,26 

26 120 NUMtII * I ♦ 26 

77 DO 121 1*25. M 

20 121 NUMtII * I • 26 

29 00 122 1*19.61 

10 122 NUMtII « I - 30 

C 

C *1* READ PUN NUMBFR AND CONTROL PARAMETER 
C 

C NRN 0 DIGIT PUN NUMBER 

C I OUT PARAMETER TO TERMINATE PROGRAM 
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I OUT ■ 0 TO REAO OATA SET 

TOUT N( 0 TO T| BMINATl PROGRAM 

JOUT UNCERTAINTY OUTPUT INDEX 

JOUT >0 OUTPUTS AU UNCERTAINTIES WITH BASIC 

OUTPUT INFORMATION I INPUT DATA, CONVERTED 
DATA, STANTON NUMBERS, ETC.) 

PIUS MAXIMUM A NO MINIMUM SINSITIVITT 
COEFFICIENTS TOR EACH VARIABLE OR PARAMETER 
JOUT « I OUTPUTS ONLY RAM AfO CONVERTEO INPUT DATA 
JOUT«f OUTPUTS ONIT MAXIMUM AND MIM1MU, SENSITIVITY 
COEFFICIENTS FOR EACH VARIABLE OR PARAMETER 
JOUT*S OUTPUTS SENSITIVITY COEFFICIENTS AT ALL 
STATIONS PLUS MAX I MUMS AND MINIMUHS 
JOUT.R OUTPUTS UNCERTAINTIES MITH BASIC OUTPUT 

PLUS SENSITIVITY COEFFICIENTS AT ALL STATIONS 
PLUS MAX I MUMS AND MINIMUHS 
••• NOTE ••• 

JOUT * S AND JOUT *A TIEL9 LAPSE AMOUNTS OP 
OUTPUT. FOR SOME CONTROL OF OUTPUT SIZE 
MOO IF T SIZE OF MAIN UNCERTAINTY LOOP FIRST. 


31 

too 

REAO ( 5 . 1000 ) I NRN( I > , 1*1 ,4 ) , 1 OUT , JOUT 

32 

1000 

FORMAT! 4A2 ,12,13) 

S3 


IF (IOUT.NE.O) GO TO CIO 


c 

C 

REAO OATA RUN DESCRIPTION, A FORMAT 

3A 

w 

REAO 15.1010) IKOMMNTIX), 1*1,351 

35 

1010 

r 

FORMAT ( 35A2 ) 


C «3« 

REAO TEST CONDITIONS 


c 

c tamb AMBIENT TEMPERATURE tOEG FI 

c PAH6 AHBIENT PRESSURE (INCHES HO CORRECTED TO 3Z OEO F) 

C RHUM RELATIVE HUMIDITY (PERCENT I 

c 

3* READ (§,10201 TAMB.PAMB.RHUM 

37 1020 FORMAT I 7F 10.0, IS 1 

36 PLATI2*0 


SR 


C 

c 

c 

c 

c 

c 

c 

c 

p 

c 

c 

c 


READ TUTFUL CONDITIONS 

TRECMV TUNNEL AIR RECOVERY TEMPERATURE (I-C TC, MVI 
POTN TUNNEL AIR VELOCITY OYNAMXC PRESSURE (INCHES H20 I 
PSTAT TUNNEL CAGE STATIC PRESSURE ( INCHES M20I 
XVO VIRTUAL ORIGIN, TBL, FROM PGM PROFILE I INCHES I 
OXVO UNCERTAINTY IN XVO, FRCM PGM PROFILE (INCHES! 

DC INI T ENTHALPY THICKNESS. FROM PGM PROFILE (INCHES) 
002INT UNCERTAINTY IN 02LOC, FROM PGM PROFILE (INCHES) 
02LOC LOCATION OF PROFILE MEASUREMENTS (PLATE •) 

REAO (5,1020) TRECMV, POTN, PSTAT, XVO, DXV0,02XNIT,0D2INT ,02LOC 


A0 
A 1 

A.’ 

Hi 

AA 

A5 

A6 


C 

C *5* 

C 

c 

c 

c 

c 

c 


toso 


lost 

1 0S2 


READ DEVELOPMENT SECTION CONDITIONS 

PLATE 1 PLATE • WHERE DATA BEGINS 

PLATE2 PLATE R WHERE HEATING BEGINS 

TOMV(I) PLATE TEMPERATURE (MV) 

HMII) HEAT FLUX METER SIGNAL IMV) 

PCOEFFI 1 1 PRESSURE COEFFICIENT 

READ (5,1030) PLATE 1 . PLATE2 
FORMAT I 213 I 
PLATE0 » PLATE1 - 1 
REAO (5. 1031 ) TOMVO 

REAO (5,10311 ITOMVdl, HMIII, PCOEFFII), I«PLATE1 ,EA 1 
FORMAT! 3F10.0I 
FORMAT! F ( 0 . 0 I 



c 

C REAO TEST SECTION CONDITIONS 

C 

C TTOPMVI I I PLATE TOP TEMPERATURE (MV) 

C TCENMV(I) PLATE CENTER TEMPERATURE (MV) 
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67 


66 

SO 


52 


51 

5<* 

55 

M 

57 

55 

55 

60 

61 


C TBOTNVH) PLATE BOTTOM TEMPERATURE (NV) 

C (Ml 1 1 PLAT! POWER I WATTS I 

c pcoirfin pressure coefficient 

PI AO 15.10601 I TBOTNVI 1 1 , TCINMVI 1 1 > TTOPNV! 1 1 .QWI It. 

1 PCOEIIII). 1*25.551 
1060 7 OP NAT (5710.0) 

c 

c HEAD RECOVERY SCCTION CONDITIONS 

c 

C TOHVIII PI ATI TEMPERATURE I I-C TC, NV) 

C MH( 1 1 HI AT 7 LUX HI TCP SIGNAL (NV) 

C PCOCri(X) PSC55UP7 COC77ICXINT 

PC AO (5.1051) I TONV( I ),HH( I ).PCOC77( 1 1,1*56,6? ) 

PC AO (5.1052) TOMV1 6 5 I 

c 

C PC AO SUPPOPT TtNPCRATiJRCS 

c 

C TSUPNV(I) TIST SCCTION SUPPOPT TINPIPATUPI (I-C TC, NV ) 

c 

PCAO (5.1050) (TSUPNV(X), 1*1,12) 

1050 70PNAT 16710.0) 

C 

C WRIT! OUT ALL paw OATA •*••••• 

c 

III JOUT .10. 5. OR . JOUT .10.2 ) GO TO 252 
WRITK6. 50001 

WRITK4, 20001 INPN(I), 1*1,61 

2000 FORMAT IIOX.'RUN * ,6A2,' ••• CURVATURC PIG ••• NASA -NAG- 5- 5* . 

1 SX, 'INPUT OATA'/J 
WRIT! 1 6, 2010 I ( KONMNTI I ) , 1*1,55) 

2010 IORNAT I 10X.55A2/) 

C 

WRITE! 6 , 2020 > TAMB.PAriB.PMUN 

2020 FORMAT! 1 0X, AMB TEMP •',76.1,' I'.SX.'AMB PRESS t'.lA.t,' IN HO*. 
I 5*. 'REL HUN *',751.' /' ./) 

WRI T£ I 6 > 20 50 ) PST AT .POTN.TPECNV ,0?INIT .02 IOC ,002 1 NT .XVO.DXV* 


62 2050 FORMAT (lOX.'STAT PRESS * , ,75. 2,' IN H20 ufN PRESS *',76.5. 

1 ' IN M20 PICOV TEMP «',F6.3,' MV' //I OX, 'ENTHALPY THICKNESS ■ 

2 .77.6, • IN AT PLATE • ' . I 5.6X, 'UNCERTAINTY «',F5.S,' IN 1 ,//, 

5 I0X, 'VIRTUAL OPIGINi X*'.76.1,' IN' »6X, 'UNCERTAINTY *', 

6 75.2.' IN',/) 

C 

65 WRITEI6.2060) 

66 2060 FORMAT( /,1 OX,' PLATE ' ,6X,'TPLAT|* ,6X, 'HN' .OX, 'CP' ,/,22X, ' 1 NV ) ' ,6X, 

1 '(NV)*,/) 

65 I * 26 

66 WRITE! 6.2050 ) PLATCO. I, TONVO 

67 WRITEI6.2050) ( I .NUNI I ) ,TOMV( 1 1 ,MM( I ) , PCOE7III), X*PLATC1.26) 

65 2050 FORMAT! 10X. 12, 16,2111 .5.110.5) 

C 

66 WRITE! 6 ,2060 ) 

70 2060 IORNAT (/, 1 0X ,' PLATE ', 7X ,' TTOP TCEN TDOT ' ,6X, ' QOOT ' , 

1 5X , ' CP ' /22X , ' ( NV ) I MV l (NV) (WATTS)',/) 

71 WRITE! 6 , 2070 ) ( I ,NUH( I ) .TTOPMV1 1 ) , TCENMVI I ) , TBOTNVI I ) ,QW( I ) , 

1 PCOEIIII). 1*25.50) 

72 2070 FORMAT ( 1 0X. 12 , 16, IX. 5710. 5. 5X.7 7. 2 ,71 0. 5 ) 

C 

75 WRITE! 6. 5000) 

76 WRITEI6, 20001 (NRNiX), 1*1,6) 

75 WRITE! 6 .201 0 ) IKOMMNT(X), 1*1,55) 

C 

76 WRITEI6.2060) 

77 WRITE ( 6,2050 ) ( I , NUNI I ) , TOMVI I ) ,HH( I ) , PCOEFKI), 1*36,62) 

75 X * 65 

76 WRITE(6.2050) I, NUnilt, TOMVI X ) 

50 WRITE! 6 ,2050 ) (TSUPMVH), 1*1,12) 

51 2050 ICRMAT (///, 1 OX, 'SUPPORT TEMPERATURES 1 NV I • ,///,26X, 

1 'TOP CEN BOT' ,//,1 OX, 'UPSTREAM ' , 5FB. 5,/, 1 OX. 

2 'DOWNSTREAM- ,575.5, //,26X, 'UP CCN DOWN',//, 

5 lOX.'TOP ' ,5F6. 5./, 10X, 'BOTTOM ',570. 51 

02 252 CONTINUE 

C 
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perturbation block 


BA 


85 

BB 

B9 

90 

91 

92 

93 
99 

95 

96 

98 

99 
102 
103 
109 

105 

106 
» 

, C" 

109 

110 
111 
112 
113 
119 

115 

116 

117 

118 

119 

120 
121 
122 
123 
129 

125 

126 

127 

128 

129 

130 

131 


INPUTS VARIED HUIVIDUALLY BT IX 
PERTURBATION LOOP INDICES* 

J MAIN PERTURBATION LOOP 

JT SUBSIDIARY LOOP TO ALLOW CALCULATION OP PARTIAL 
DERIVATIVES MITM RESPECT TO TII*1> AND TII-11 
AS WELL AS Till 

JS SUBSIDIARY LOOP FOR SI 1*1 I AND Sill 
NTAGIU MARKS ARGUMENT IN SUBSIDIARY LOOPS. 



DISTINGUISHING AMONG I, 1*1 , AND 1-1 FOR 
IN COMPUTING SENSITIVITY COEFFICIENTS. 

LATER USE 

DATA 

OESC/ 1 TOI I ' , 

• » 

•TTOP’.MII '.'TCEN'.MIl 

, 

1 

•TBOT’ 

, Mil 

.•Bill*.' •.•Sill'.* 


2 

•HMIl' 

,' 1 

* *QI 1 1 ' » ' 


3 

'All)’ 

1 

.'PAMB',* ' , 'RHUM' i • 


9 

•TREC* 

, *ov 

• • POYN' , ' • , ' EMIS* , ' 


5 

•MCAL* 

* 

i 'KOI 1 1 • 1 1 * , ' RBOl ' » ' 1 1 


6 

'RLOO * 

(Mil 

,’RA ' , 'PSTA* . *T 


7 

•SPAN* 

• Mil 

• ' T5UP' » 'PI 0 ','RFW 

• . 
t 

8 

•TCAL* 

, 

,'PS * , •RHOS' • * 


9 

•XA • 

, 

,'RV * , *RB 


1 

'KUP 1 

i 

, 'KOOW' , *N V 


DATA 

UNCERT/900 . 0 , 

9*0.00 

, 0.003, 0.01, 0.020, 0.05, 0.02, 

1 

0.0005, 0 

.30, 1 

00. 0.003, 0.00500, 0.02, 

0 . C5 , 3* . 01 

2 

0.005, 

0.005 

20., 0.009, 1.09, 0.1, 2 

.0, 0.0001 . 

3 

0.003, 

100.0 

.005, 0.0066, 0.01790/ 



JS«1i JSC«11 JT*1 

AMO* 100.0 

00 160 I=PLATE1 .63 

ARGI 1 1*100.0 

CONTINUE 

00 500 J*1 ,33 

JCOUNT*J 

00 191 JS*1 .2 

JSC*JS 

NT*2 i MT*62 

TOMVI 69 1*0 . 0 

TTOPMVI 39 1*0.0 i TCENMVI 391*0.01 TBOTMVt 39 1*0 . 0 

00 19: JT* 1 ,3 

IFU.GT.6) GO TO 295 

IF1J.EQ.1 I GO TO 280 

JTC*JT 

IFUT.EQ.3I MT*63 
IF(NT.NE.O) GO TO 225 

IFIJ.EQ.2) CALL CHANGE! O.JTC.TOMVO, TOMVI 1 1 1 
1FIJ.EQ.6I CALL CHANGE! 0 . JTC.BO ,BI 1 1 1 
NT* 3 
CONTINUE 

00 193 I=WT ,MT ,3 
ICOUNT*I 

IF( I . LE . 38. AND . I .GE . 25 1 GO TO 220 

IFIJ.EQ.21 CALL CHANGE! ICOUNT ,JTC, TOMVI 1 1. TOMVI 1*1 1 1 

IF! J.EQ.3.ANO.JT.EQ.1 1 GO TO 223 

GO TO 229 

TOMVI 29 )=ARG< 291 

TOMVI 631*ARG(631 

TOMVO=ARGO 

IFII.EQ.90.0R.I.EQ.39I GO TO 229 

IFII. GT.21 TOMVI 1-2 »=ARG( 1-21 

IFU.EQ.61 CALL CHANGEl ICOUNT, JTC .B( 1 1 ,31 1*1 1 1 
GO TO 221 

IFU.EQ.31 CALL CHANGEl ICOUNT, JTC, TTOPMVI 1 1 .TTOPMVI 1*1 1 1 

IFI J. EQ. 9. ANO.JT.EQ.1 .ANO.I.GT.26 1 TTOPMVI 1-2 )=ARGI 1-2 1 
IFI J.EQ.9.AN0.JT.EQ.1 1 TTOPMVI 25 >= ARGI 25 1 

IFIJ.EQ.91 CALL CHANGEl ICOUNT , JTC,TC r NMVI 1 1 .TCENMVI 1*111 
IFIJ.EQ.5.AN0. JT.EQ.1 . ANO.I.GT.26 1 TCENMVI 1-2 1 *ARG( 1-2 1 
IF I J. EQ. 5. AND. JT. EQ. 1 1 TCENMVI 25 >*ARG( 25 1 
IFIJ.EQ.51 CALL CHANGEl ICOUNT , JTC , TBOTMVI 1 1 , TBOTMVI !♦ 1 1 1 
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1 t 
» 

51 
16 

17 

16 

19 

<*0 

91 

9E 

91 

96 

96 

97 
96 
99 
60 
61 

52 
61 
69 
65 

56 

57 
55 

59 

60 
61 
62 
61 
69 

65 

66 
67 
66 

69 

70 

71 

72 


175 

179 

175 

176 

177 
176 

179 

160 

161 

182 

161 

169 

165 

166 
167 
166 
169 

190 

191 

192 
191 
199 
195 
1 96 


1MJ.NE.6) CO TO 221 

IMJT.EQ.1ANO I. ST. 26 I TBOTMVI I-t >*ARG( I-f I 
1MJ1.1Q.1I 1601 MV I 25 I* ARGI 25 I 
ARGI I . 0 

Ifl I IQ. 01 GO 10 12 6 
(ft HI 6GI 1 I »• 

22 6 NIAGO'O 

in I. IQ. 611 GO 10 221 

HI Alt' 1*1 !■ - 1 
IMI.IQ.1 I GO 10 19) 

222 NTAGt 1-11*1 
1*1 CONTINUE 
295 inj.rq.7i BO* ARGO 
DO 199 1*1 ,95 
1C0UN1 * I 

in J.IQ. 7. AND JS.fQ. 1 I 6(1 I® ARGI 1 I 
IMI.Gt 25. AND. I. IE . 16) GO TO 29 7 

IMI. GT.62I GO TO 199 
iriJ.lQ.6) CALL CMANI 1COUNT.HMI I II 
irtJ.IQ.9l MM1 1 ) * ARGI I I 

GO 10 199 

297 IMJ. IQ 91 CALL CMANI t COUNT ,QM( I II 
inj.IQ.10l QWI I I* ARGI I I 
199 CONI I NUt 

II i J N.ltl CALL CMANI 1 . PAMB 1 

IMJ IQ HI (All l MANGIt 1 .9.RMUM.PAMBI 

IMJ fq 191 CALL CMANGM 1 .9.TRICMV.RMUM) 

.ii • l« 111 CALI CM ANSI 1 1 . . pd »n. tre cmv » 

in J.IQ. 16 1 POTN- ARGI I I 

iriJ.IQ.2CI CALL CMANI 1 .PSTAT I 

II I J I q mi rsi AT ARGI 1 I 

II IJ.Nl.21 . ANO . J.NC.29I GO TO 195 

DO 196 1-25.16 

ICOUNT«I 

II U.IQ.21I CALL CHANI ICOUNT .DILQCI I 1 1 
II 1 J.EQ.29) OILQCI 1 l*ARGI 1 1 
196 CONI INUl 
195 CONTtNLit 

IMJ. IQ. 291 CALL CHANI 1 , TSUPMVI 1 0 II 
IMJ. IQ. 251 TSUPMVI 10 !■ ARGI 1) 

260 CONI INUt 


C 

c ••••*>« DATA CONVIRMON BLOCK 

C 

in J.GT . 1 ) GO TO 250 

C AVOIDS CONVERTING PRl SSURI AND HUMIDITY MORI THAN ONCE 

RMUM * WMUM / 100. 

C CONVIRT ALL PRl 3SURCS TO LBI/SQH 

POTN » POYN • 5.195 
PSTAT ■ PSTAT • 5.195 
PAMB * PAMB * 70.726 


250 CONTINUr 

C CONVIRT Ail lirn’tWATURTS TROM MV TO DIG f 

IM J.tQ. 7. AND.J5.IQ.il GO TO 25) 

tri NT,TtICMV»t.t»TAMil 

ICC : TCI JCOUNT , TOMVO ,B0 .IAMB 1 
DO 100 I *1,1 1 

TSUPPt I 1 »TC( JCOUNT .TSUPMVI 11,0.0, T AMB I 

100 CONTINUE 

00 101 1 * Pi ATE 1 ,29 

TOl I > • TC I JCOUNT. TOMVI I), Bl I i.TAMB) 

101 CONTINUE 

DO 102 1*25,56 

T TOPI I 1 * TCI JCOUNT, TTOPMVI II. Bl I 1.TAMB1 
TCINI II TCI JCOUNT, TCENMV I I I .Bl I 1 . T AMO I 
TBOTl I I - TC I JCOUNT .TBOTMVI I I .Bl I I , T AMB I 

102 CONTINUE 

DO 101 1*39.6) 

TOl I )*TC( JCOUNT. TOMVI LI. Bl I I, TAMO) 

105 CONTINUE 
255 CONTINUE 
C 



197 

190 
1 99 

too 

SOI 

202 

20 ) 

209 

205 

20 * 


207 


too 

209 

210 
211 
212 
213 
21 '. 
2*4 
216 
217 
216 
219 


220 
221 
222 
22 ) 
KM 
22 5 
226 


227 

226 

229 

230 


231 

232 

233 

236 

235 

2 36 

237 

236 

239 

290 

291 

292 

293 
299 


C COM PUT! TM| Tf MPERATURES IN TNI FRAME NEAR THE ENDS OF TNE PLATES 
C IV A QUADRATIC FIT TO TNE TEMPS. MEASURED IN TNE PLATE SUPPORTS. 

C ALSO COF1PUTE AVERAGES OF END AND OVERAI L TEMPS. OF TNE PLATES. 

C 

At ■ 1 1 TSUPPI 7 )♦ TSUPPI 1 0 1-TSUPP1 6 1-TSUPP* 111 1/S. - ITSUP*»I6)» 

1 TSUPPI 1 1 I -TSUPPI 9 I -TSUPPI It 1 1/6. 1/22. 

At ■ I TSUPPI 61*1 SUPPI 1 1 1- TSUPPI 7 ) -TSUPPI 1 0 M/10. - A2 • 9. 

AO ■ I TSUPPI 7XT3UPPI 10 1 1/2. • A1 • 2. - A2 * 9. 

c 

00 109 1*25. 36 

TFRAMf I I I ■ A >i ♦ A1 « (1-29) • At « ( 1-29 )•( 1-29 I 
TOI 1 1 ■ Term 1 1 

TENOIII ■ I T TOPI I > ♦ TVOT (III /2. 

109 CONTINUE 

c 

TUP ■ I TSUPPI I) ♦ TSUPPI 2) ♦ TSUPPI SI) /I. 

TOOUN » (TSUPPI 9 I ♦ TSUPPI 5 I « TSUPPI 6 I ) /). 

C 

C COMPUTE NINO TUNNEL FLOW CONDITIONS 

CALL TUNNEL IJCOUNT ) 

C 

C COMPUTE NET ENERGY TRANSFER FROM DEVELOPMENT 

C REGION. TEST SECTION AND RECOVERY REGION 

CALL POWER I JSC • JCOUNT .TINF , TUP, TOOtRI. PLATE 1 I 
C CORRECTION FOR SPANWISE MEAT FLOW ANO UPOATE PLATE CENTER TEMP 
DO 106 I«1 ,62 
DOE LOCI X )»0.0 
106 CONTINUE 

DO 109 1*25.30 

TFC=0.0221*l TBOTI I J-TCENI I I )-0. 0127*1 TCENI I )-TTOP( 1 1 ) 

TFB*-0 .192*1 TCENI I l-TTOPI ID-9. 06*TFC 
TOI I l-TBOTI I >*5.S6*TFB-30.9*TFC 
DELQCI I 1*919. 0*TFC 
DDLLQCI I )-0. 10*DELQC( I 1 
QOOTI I XQOOT 1 1 XOELQCI X ) 

109 CONTINUE 
C 

C •■*•••• WRITE ALL CONVERTED DATA •*•*••• 

C 

IFIJ.GT.1 ) GO TO 251 

IF I JOUT . CQ . 3 . OR . JOUT . EQ . 2 ) GO TO 251 

WRITE! 6 , 3000 ) 

WRITE! 6 1 2 1 90 1 I NRNI I ) , 1*1,9) 

WRITEI6.2010 > (KOMMNTII), 1*1,35) 

WRITEI6.2100) 

2100 FORMAT! /, 1 0X, 'PLATE ' ,6X, ' TPLATE ' ,5X, ' MM ' , 0X , ' QF LUX ' ,/ , 

1 22X, '(F)' ,6X, ' I MV ) I BTU/NR/SQFT ) * ,/ ) 

1 * PLATEO ♦ 29 

WRITE! 6,2110) PLATEO. X. TOO 

WRITEI6.2110) 1 1 ,NUM( I ) , TOI I ) ,HMt I ) ,QOOT I X ) , X*PLATE1 .29) 

2110 FORMAT ( 1 0X , 12 , 19 ,F 1 0 . 2 ,F 1 0 . 3, F 1 1 . 2 ) 

C 

WRITE (6,2120) 

2120 FORMAT (/, 1 OX, ' PLATE ‘ ,6X, • TPLATE 1 ,5X, • QOOT • ,6X , 'QF LUX’ , 

1 6X, 'TrRAME * ./.22X. 'IF )' ,5X, * I WATTS ) (BTU/NR/SQFT) IF)’,/) 

WRITEI6.21 30 ) I I .NUMI I ) . TOI I ) ,Q< I ) .QOOTI I I , TFRAME I I 1 . 1 = 25,30) 
21)0 FORMAT! 1 0X , 12 , 19, 2F 1 0 . 2 , F 1 1 . 2 ,F 1 2 . 2 ) 

C 

WRITE! 6 , 3000 ) 

WRITEI6.2190) INRNIX), 1*1,9) 

2190 FORMAT MOX.'RUN * , 9A2 , • ••* CURVATURE RIG •*« NASA-NAG-S-3 *, 

1 5X, 'CONVERTED DATA'/) 

WRITEI6.2010) (KOMMNT(X), 1*1,35) 

C 

WRITEI6, 21001 

WRITE I 6.2110 I ( I .NLRK I ) , TOI I ) ,HM( I ) .QOOTI X), I* 39, 62) 

X * 63 

WRITE! 6.2110) I, NUMI I), TOI I ) 

WRITE! 6,2150) (TSUPPIX). 1*1,12) 

2150 FORMAT (///,1 OX, 'SUPPORT TEMPERATURES IDEG F)*,///,29X, 

1 'TOP CEN BOT' ,//,10X, 'UPSTREAM ' • 3F6 . 2 ,/, 1 0X , 

2 DOWNSTREAM' ,3F0. 2,//, 29X, 'UP CEN DOWN',//, 
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3 1 OX i ' TOP ’ . 3F6.2./.10X. 'BOTTOM * . V 8 . 21 

24 5 251 CONTINUE 

C 

c ■•••••• COMPUTE STANTON NUMBER ••••••• 

c 

C X REYNOLDS NUMBER BASED ON VIRTUAL ORIGIN TBL 

246 FACT*UINF/IVISC*1t. I 

247 DREX*FACT*DXVO 

246 00 HO 1*1.62 

249 110 REXIII«FACT«tXII) *V0 !•< 1 . 0-PCOEFFI 1 1 1»*0 . 5 

C COMPUTE STANTON NUMBERS 

250 OtNON=RHOINF*UINF*CP«3600. 

251 00 lit 1*1 .62 

252 ST I I )*QOOT< I 1/tOENOM*! T0( X l-TAOIAB )*l 1 .O-PCOEFFI X 1 >**0.5I 

253 111 CONTINUE 

C «••••• COMPUTE UNCERTAINTY IN STANTON NUMBER •••••• 

C 

c 

C OST I I I ORIGINAL ( UNPERTURBED 1 STANTON NUMBER 

C 05T1I1 CHANGE IN STANTON NUMBER OUE TO PERTURBATON 

C OARG CHANGE IN VARIABLE OR PARAMETER DUE TO PERTURBATION 

C POlII PARTIAL DERIVATIVE DST/DARG 

C SENCOI I . J ) SENSITIVITY COEFFICIENT P0( I )*ADSOLUTl UNCERTAINTY 

C IN VARIABLE OR PARAMETER 

C 

254 IF! J.EQ.1 1 GO TO 264 

255 DO 160 1*1 .62 

256 05TI I )*03T( 1 1-ST1 X I 

257 IF I J.UE.12.AN0.J.LE.16I GO TO 290 

256 IF< J.GE.21 .ANO.J.LE.25 .ANO. J.NE.23) GO TO 290 

259 IFU.GE.27I GO TO 290 

260 XFIJ.GE.33I GO TO 290 

261 GO TO 261 

C FOR VARIAB1ES OR PARAMETERS APPLYING TO ALL PLATES. 

C OARG IS MADE THE SAME AT EACH PLATE 

262 290 DO 161 IARG«1 .62 

263 161 ARG( IARG )*ARG( 1 1 

264 281 CONTINUE 

265 DARG=ARG( 1 1*0.01 

266 IFIJ.GT.71 GO TO 262 

267 IF ( NT AG( Xl.EQ.il DARG= ARG( 1*1 1*0.01 

268 IFINTAGI I I.EQ.-1 1 0ARG*ARG( 1-1 1*0 . 01 

C FOR J*6, B( 25 l*B( 26 I ALL EQUAL ZERO 

C FOR J--28, TM1 THERMOCOUPLE CALIBRATION EQUATION RESULT IS 

C PERTURBED BY 0.05F RATHER THAN V/. TO SIMPLIFY DETERMINATION OF ARC 

C FOR J DM 11 = 0.0 

269 262 IF I J. EQ.6. ANO . 1 . LE . 36. ANO . X.GE .25 I DARG*0.0001 

270 IF1 J. EQ.6. ANO. I. LE. 36. ANO. I. GE. 251 DARG*0.0001 

271 III J. EQ.6. ANO. JT.EQ.1 . ANO . I . EQ . 39 I OARG = 0.0001 

272 IF! J. EQ.6 . ANO . JT. EQ. 2 . ANO . I .EQ. 24 I DARG*0.0001 

273 IFIJ. EQ. 261 OARG=0.05 

274 IF IABSIDARG) .IT. t.OE-81 OARG=1.0E-6 

275 POI I)=DST(X)/OARG 

C UNCERTAINTIES IN Kill AND SHI CALCULATED IN SUBROUTINE POWER 

C USING RELATIVE (PERCENT) UNCERTAINTIES 

276 UNCERTI 1 0 )=DKK( 1 1 

277 UNCERT17I*DSS(II 

276 UNCERTI 6 )*UB( II 

279 UNCERT I 23 )=DDELQC( 1 1 

260 IF II GE. 381 UNCERTI 23 )*DOELQC( 36 1 

261 POU--POI I l*UNCCRT( Jl 

C COMPUTE SENSITIVITY COEFFICIENTS 

262 SENCOI I,JI=ABS(POU» 

263 445 CONTINUE 

C DETERMINE ARRAYS IN SUBSIDIARY LOOPS 

264 IFIJ.GT.6I GO TO 283 

285 SCTI I . J . JT 1* SENCOI X . J I 

286 283 IFIJ.NE.7) GO TO 160 

267 SCTlI,J,JSI=SENCniI,J» 

286 160 CONTINUE 

269 NT=NT-I 

290 MT=MT-1 
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291 

292 

295 

296 

295 

296 

297 
290 
299 

500 

501 

502 

505 

506 

505 

506 

507 
500 

509 

510 

511 

512 

515 

516 


515 

516 

517 
510 

519 

520 

521 

522 
525 


526 

525 

526 
52 7 

520 

529 

550 

551 

552 

555 

556 

555 

556 


557 

550 

559 

560 


561 

562 
565 
5h6 

565 

566 

567 
560 
569 



IFIJ.GT.6) GO 

TO 292 

162 

CONTINUE 


292 

IFIJ.NE.7) GT 

TO 291 

161 

CONTINUE 


291 

CONTINUE 
GO TO 500 


2B6 

00 165 1*1.62 


165 

OSTI I )*STl 1 ) 

NT*2 

MT*62 


500 

CONTINUE 



SCTI 26,2,1 )*SQRTI SENCOI 26 , 5 )*SENCOI 26 , 5 ) ♦SENCOI 26 ,6 )*SENCO( 26,6 ) 


1 < SfUCOi 26.5 t*SENCO( 26,5) 1 

5CT( 59,2.1 > * SOW T f SENCOI 59.5 l*SENCO( 59, 5 (♦SENCOI 59.6 )*SENCOI 59,6) 

1 ♦SENCOI 59,5) »SENCOI 59,5)) 

DO 102 J*2 ,6 
00 102 1*1 ,62 

102 SENCOI I. J)*(SCT(I,J,1 >«5CTII,J,1 •♦SCTI I,J,2 >*SCTI X.J.2 ) 

1 ♦SCTII,J,5)*SCTII,J,5>)**0.5 
00 105 1*1 ,62 

105 SENCOI I , J )*( SCTI X > J, 1 >*SCTII,J,1 MSCTI X , J,2 )*SCT( X.J.2) )**0 . f 
00 171 X*1 .62 
SUM*0 . 0 
OO 170 J*2 . 55 

170 SUH*SUN« SENCOI X , J )*SENCOI I • J ) 

USTABSI I )*SUM**0 . 5 

1 71 USTRE LI X >*USTABSI I l/STI 11*100.0 
C 

C COMPUTE OEL2 AND REOEL2 BASED ON ACTUAL ST-DATA 

CALL ENTHAL I FACT .02INIT .D02XNT , 02LOC ) 

00 172 1*1.62 
REENPI 1 1 = 0.0 
STPI I )*0.0 

IF IX ,LT. D2LOC ) 60 TO 172 
NP*I-D2L0Ca1 
REENP(NP)*REEN( X ) 

STPI NP )*OSTI X ) 

172 CONTINUE 
C 

c •••«•«• WRITE STANTON NUMBER DATA •*«***• 

C 

IFI JOUT.LE.5.AND.JOUT.GE.1 ) GO TO 205 
WRITE! 6 , 5000 ) 

WPITEI6.2200) (NRNtI), 1*1,6) 

2200 FORMAT 11 OX, 'RUN ' ,6A2 , ' «•* CURVATURE RI6 NASA-NAG-5-5 *, 

1 5X, ’STANTON NUMBER DATA’/) 

TA0BC=5 . *1 TAOIAB-52 . )/9. 

TINFC*5.*I TINF-52. )/9. 

UINFMS=UINF*0 . 5060 
XVOCM=XVO*2 . 56 
RHOKM5*RHOINF*16 . 02 
VISCI*VISC»0 . 0929 
CPJKGK*CP*6166. 

WRITEI 6,2210) TAOBC .UINFMS , TINFC .MHOKM5 , VISCI .XVOCM .CPJK6K , PR 
2210 FORMAT! 10X. ’TADB*’ ,F6.2, ’ DEG C UREF*' ,F12.2, ' M/S TINF* ' , 

1 F6 . 2 , ' DEG C VI OX, 'RHO* ' ,F7. 5, * K6/M5 VISC* ’ ,E 1 2 . 5 . ' M2/S’, 

2 6X, 'XVO* ' ,F7. 1 , ' CM’/IOX, ’CP*’ ,F6.0, ’ J/KGK PR*’,F16.5/) 

WRITEI 6 ,2220 ) I KOMMNTI I ) , 1*1 , 55 > 

2220 FORMAT! 1 OX, 55A2/) 

WRITEI 6 ,2250 ) 

2250 FORMATI 1 OX, ’PLATE* .SX.’X' ,5X, ’UPW* ,5X, 'K*,5X,6X, ’REXVO’ , 7X , * TO* ,6X 

1 , 'STANTON NO' ,6X, ’RE ENTH’ , 7X , *0ST ’ , 9X , ’OSTI '/. ) * , 7X , * OREEN* ,/ 

2 16X, ' I CM » *,* I M/S 1 * ,27X, ' DEG C',/) 

DO 115 1*1,62 

POSTI I )*05Tl X ) 

XCMI I )*XI X )*2 . 56 
PXCMI X )*XCM( X ) 

IF II .EQ. 25) WRITEI 6 ,6000 ) 

IF IX .EQ. 59) WRITEI 6,6100) 

6000 FORMATI /I OX, 'CURVE BEGINS’/) 

6100 FORMATI /I OX, 'RECOVERT BEGINS'/) 

IF (I.EQ.59) WRITEI 6 , 5000 ) 
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350 

351 
35? 
353 

356 

355 

356 

357 
155 

359 

360 

361 

36? 

363 

366 

365 

366 

367 


360 


369 

370 

371 
37? 


373 

176 

375 

379 

381 

38? 

386 

387 
391 
39? 
396 
398 

J99 

600 

601 

60? 

603 

606 


6^5 
6 6 

607 

608 

609 

610 
61? 
613 
616 


IF (I .67. II XCML*XI 1-1 1*2.56 

uioc*uiNrn5»( i . o-pcoeff i x i i**o.s 
TEMPC*5.*I TOI n-S?. 1/9. 

xr II .6T. II ACC*I 1 00 . *1 VISCX/1 ? . 0*ULOC I !• 

1 1 1 prot i f i i-i i-pcoeffiiii/ixchiii-xcfiliii/ii.o-pcoeffiiii 

IF 1 1 .EQ.1 I ACC *0.0 
XF II .6? . 0?10CI GO TO 116 
IF IX .61. PLATE? I CO TO 115 
WRITE I 6 > ??65 I I.XCHIII.ULOC.ACC.REXIII.TEMPC 
??65 F OPM AT I I OX. 13. F 7.1 ,F7. t . 1 X .C 1 0 . 3 . 1 X .E I ( . » .F 7. C I 
CO TO 113 

115 WPITEi 6 >??60 I I «XCMl 1 1 , UlOC . ACC > RE XI 1 I .TEMPC .OSTI X I .USTABSI I I .USTR 
tEU I ) 

2?60 FORMATl 10X,X3.F7.1,F7.?,1X.E10.1.1X.E1?.5,F7.?.?X.E1?.5,16X,E9.3, 

I 3X.f9.3l 
GO TO 113 

116 WRITEI6.22S0 I X , XCMI 1 1 ,ULOC , ACC ,PEXI 1 1 , TEMPC .OSTI X I .REENI 1 1 .USTABS 
I I 1 1 .USTREll 1 1 .DREENI 1 1 

??50 FORMATMOX.I3.F7.1 .F7. ? . 1X.E1 0 . 3. IX.E1 ? . 5.F7.? .21 EX.E1 ? . 5 I .?X. 

1 E9. 3.3X.F9. 3.6X.F7.0) 

113 CONTINUE 

MRITEI6.E260I OREX 

2260 rORMAT I/12X. 'UNCERTAINTY IN REX*',F6.0I 

r 

?85 CONTINUE 


C 

C ••••• ADDITIONAL UNCERTAINTY INFORMATION ••••• 

C 

C COMPUTE AND WRITE OUT MAXIMUM A NO MINIMUM SENSITIVITY 

C COEFFICIENTS FOR EACH VARIABLE OR PARAMETER 

c 


XFI JOUT.EQ.1 I GO TO 293 
WRITEI6.3000 I 

WRITEI6.2360 KNWNI I1,I«t .61 

2360 FORMAT! 10X, 'RUN' ,6A?, ' ••• CURVATURE RIG ••• NASA-NAG-3-3 '. 

1 5X, 'SENSITIVITY COEFFICIENTS IOST/OAPG • DELTA ARGI'//, 

? 50X. 'PREPLATE' 16X,*TC5T SECTION' , 1 5X, 1 AFTERPLATE * //SX, 

O' ARGUMENT ’.25X, 

9 'MAX MIN MAX MIN'.IOX. 

1 'MAX MIN'//) 

DO 186 J*?,33 
JCOUNT*J 

SMAXORI J 1*25.0 iSMINORIJ 1*25.0; SMAXTSI J 1*25.0 ISMINTSIJ I *25.0 
SMAXRRIJ 1*25.0 iSMINRRIJ 1*25.0 
CALL MAXMINM .26, JCOUNTI 

SMAXORI JlsVMAXJ SMINDRI J l«VMIN5 N1 *NA i N2 =NB 
CALL MAXMINI 25. 38. JCOUNTI 

SMAXTSI J l*VMAX; SMINTSI J UVMINi N3*NA i N6 = NB 
CALL MAXMIN1 39.62 > JCOUNT I 

SMAXRRIJ l*VMAXi SMINRRI J )*VMINS N5*NAi N6 = NB 

JOE5C1=2»J-3; J0ESC2=?*J-2 

WR ITEI 6 . 2300 IDESCl JOE SCI I .OESCI JDESC2 I .SMAXORI J I .SMINDRI J I . 

1 SMAXTSI J ) .SMINTSI J 1 . 5MAXRRI J I .SMINRRI J I 
2300 FORMAT! 6X, 2A6.26X .E 1 3. 5. 5E 1 3.5 I 

WRITEI6.2310 IJ.UNCERTI J I .N1 .N2.N3.N6.N5.N6 
2310 FORMAT l 6X, ' I J*',I2,' DELTA ARC* ' ,F9.5, ' )' . OX, '(•' .It,' )' • 
1 51 OX. ' I • ' ,12,' I' II 
106 CONTINUE 

WRITEI6.3000) 

293 CONTINUE 
C 

C WRITE OUT ALL SENSITIVITY COEFFICIENTS 

C 

IF I JOUT.LT. 3 I CO TO 286 
WR ITEI 6 > 3000 I 

MR ITEI 6. 2 350 1 1 NRNI 1 1, 1*1 ,61 
DO 190 J*2 • 33 
WRITEI6, 30001 

J0ESC1=2*J-3; JOESC2*2*J-2 

WR I TE 1 6 , 2351 I OESCI J0ESC1 I .DESCI JOESC2 I 

DO 190 1 = 1 .62 

190 WRITEI 6 .2352 I I .SENCOI I , J I ,ST( X I 
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<•15 


2)50 


FOBMATHOX, RUN’,AA2,* • •• CURVATLWE BIG ••• NASA NAG-)-) 
1 SENSITIVITY COEFFICIENTS’/) 


AH 

2351 

FORMAT! 1 OX, 'PLATE SCNSITIVITY 1 OST/O ' ,2AA, • » STANTON 

1 number’ ) 

AIT 

2)52 

FORMAT! 12X,I2,11X,112.A.12X,I12.A1 

<•18 

28 * 

CONTI NUC 

AI9 


NUIIP1 *A2-PLATE 1 

•.20 


XF (PLATE2.GT. PLATE) 1 NUMP1 *A2-PLATI2 

A2I 


IF IPLATI2 .IT. PLATED GO TO 287 

A22 


DO 288 I<PLATI2.A2 

AS) 


NOFF*I-PLATE2*1 

h2H 


POST1NOFF ) -OSTI I ) 

A25 


PXCM1NOFF I*XCMIX) 

A2A 

288 

CONTINUE 

417 

287 

CONTINUE 

<•28 


NUMP2*A2-D2LOC 

A29 


IF IIRUN .EG. 11 CALL PICTRGIO.Ot 

A)0 


CALL SETSMG1 21.0.01 

A ) 1 


CALL SETSMGI 2A.0.0I 

A)2 


CALL SUOJEGI-AO. .0 . 0 . 1 30 . 0.0. OOAO 1 

<♦33 


CALL 5ETSHGI 130.0.01 I 

A3A 


CALL SETSMG1 131 .0.01 I 

4 35 


CALL SETSMGI 132,1 .0) 

A 56 


CALL CSETGI ’OUPLEX#’ 1 

A)7 


CALL GBAPHGI ’•* .O.NUHP1 , PXCM, POST . ) . ’STREAHMISE DISTANCE (CM)#’, 
1 ’STANTON NUMBER# ’ . • STANTON NUMBER VERSUS STRIAMUISE DISTANCE#’) 

A 38 


CALL PICTRGIO.O) 

*•39 


CALL SETSMGI 23.1 .01 

AAO 


CALL SETSMGI 2A.1 .01 

AA1 


CALL CSETGI ’DUPLEX# • 1 

AA2 


CALL SUBJEGI 200. ,0.0005,5000. ,0.005 1 

Ah) 


CALL SETSMGI 133,2.0) 

AAA 


CALL SETSMGI DA. 2. 0 1 

AA5 


CALL GRAPHGI ’•’ ,0 .NUNP2 .REENP.STP, 3, 'ENTHALPY THICKNESS REYNOLDS N 
2UMOER#' , 'STANTON NUMBER# STANTON NUMBER VERSUS ENTHALPY THICKNES 
35 REYNOLDS NUMBER • ' 1 

AAA 


IRUN* 1 

AA7 


GO TO 200 

AAO 

210 

CONTINUE 

•.<•9 

3000 

FORMAT I 1M1 1 

A50 


CALL EXITG 

A51 


RETURN 

A52 


END 


A53 

FUNCTION TCI J.EMF .B.TAMB ) 


C 

FUNCTION CONVERTS TEMP 

FROM IRON-CONSTANTAN MV TO DEG F 

A5A 

T « 32.350 ♦ 35.300 • EMF 

- 0.A05 * EMF • EMF 

hSS 

1FI J.EQ.26 1 T=T«0. 05 


A5A 

TC«T*B*IT- TAMB) 


A57 

RETURN 


A5S 

END 



F-ll 



459 


<*60 

461 

462 


46 ) 

464 

465 

466 

467 

466 

469 

470 

471 
47t 
475 
* 7 * 

475 

476 

477 
476 
479 

460 

461 

462 
46) 
464 
46) 

| 

467 
466 
*59 

490 

491 


492 

49) 

*94 

49) 

496 

497 
496 


499 

)00 

501 

502 


SUBROUTINE TLAR4EL « J1 


THIS ROUTINE COM PUT I 3 THE HIND TUNNEL FLOW CONDITIONS 

UINF FREE STREAM VELOCITY (FT/SEC) 

TINF FREE STREAM STATIC TEMPERATURE (DEC FI 

RHOINF FREE STREAM DENSITY I IBM/FT)) 

VISC FREE STREAM KINEMATIC VISCOSITY (FTt/SECI 

CP FREE STREAM SPECIFIC HEAT (BTU/LBM/OEG R» 

PR FREE STREAM PR AND T l NUMBER 

M FREE STREAM ABSOLUTE HUMIDITY I LBM H20/LBM ORT AIR) 

COMMON/ BIK1 /PAMB.PSTAT.TRECOV.RHUH.POYN.TAMB 
COMMON/ BLK2 /UINF .TINF , TAOIAB .RHOINF , VISC . PR , CP, M 
COMMON ARC) 6)). ARGO 


OATA 


SATURATION DATA FROM K ANO K 1969 STEAM 
DIMENSION TEMPI lOMPSATI 101.RHOSATI 10) 
REAL MFA.MFV.MUA.MUV.JF 
40. . 

100 . 0 . 

17.519. 

1)6. 64). 

.0004090. 

.0026571 , 


TABLES 


DATA 


OATA 


50.0. 

110 . 0 . 
25.6)6, 
165.767. 
.0005666. 
.00)7722, 


60.0, 
120 . 0 , 
56.907, 
244.006. 
.0006266, 
.0049261 . 


70.0. 
1 ) 0 . 0 / 
52.501 . 
520.400/ 
.0011525, 
. 006)625/ 


10 

20 


TEMP/ 

90.0, 

PSAT/ 

100.627, 

RHOSAT/ 

.0021)61, 

DO 10 N*1 ,9 

IF ( TEMPI NI.GT.TAMB) GO TO 20 
CONTINUE 

OT * (TAMS - TEMPI N-1 M/10. 

PS * PSATIN-11 ♦ DT*( PSAT 1 N I - PSAT (N-1 II 

RHOS * RHOSAT I N-1 1 ♦ DT«1 RHOSATI N I - RMOSATI N - 1 1 1 

IFU.CQ.27I CALL CHANIt.PSI 

IFU.IQ.t6) CALL CHAN! 1 .RHOS I 

RA»1 545. 52/26. 970 

P « PAMB ♦ PSTAT 


60.0, 


75.051, 

.001560), 


RHOV * KHUM*RHOS 
PVAP * RHUM*PS 
PA * P - PVAP 

RHOA * PA/l RA*( TAMB ♦ 459.671) 

H*RHOV/RMOA 

RHOM * RHOA ♦ RHOV 

MWA ■ 26.970 

MUV a 16.016 

MFV * RHOV/RHOM 

MFA « 1 .0 - MFV 

RM ■ 1545. 52*<MFA/MWA ♦ MFV/MMV ) 


CP « MF A*0 . 240 ♦ HF V*0 . 445 
GC = 52 . 1 7)9 
JF*776. 26 

RFU: RECOVERY FACTOR FOR HIRE 

RFU : 0 .66 HIRE NORMAL TO FLOU 

RFU-O. 66 HIRE PARALLEL TO FLOU 

RFU-0.66 

IFU.EQ.25) CALL CHANI 1 ,RFH 1 

TINFsl TREC0V*459.67)/(1 . ♦ RFH*POYN*RM/(JF*CP*P I ) - 459.67 
RHOINF = P/RM/( TINF *459.67 ) 

UINF«SQR T( 2 . *GC*PDYN/RHOINF 1 

VISC*II1 .*0.0175»TINF I/I 1 .E06*RHOINF )*( 1 .-.7*H) 

PR=. 710*1 5)0. /I TINF *459. 67 I |RR| . 11*1 I ,*.9*HI 
CONVERT TO AOIABATIC HALL TEMPERATURE 
RCF:PR*»0. 5)555 

TADI AB = TINF *RCF*UINF*UINF/I 2 . *GC*JF»CP I 

RETURN 

END 
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SOS 


50 * 


505 

506 

507 
500 
509 


510 


SUBROUTINE POMCR I JSC. J,T INF, TUP, TDOWN. PLATE 1 ) 

THIS ROUTINE • 

Ml CORRECTS THE INOICATEO PLATE POWER READING FOR 

WATTMETER calibration ano circuit insertion losses 
III COMPUTES NET ENERGY LOST FROM PLATES BY FORCED 
CONVECTION HEAT TRANSFER 

IS) COMPUTES HEAT FLUX FROM RECOVERY REGION PLATES 

COMMON/ BLKS /TOO.TOI 6SI ,HM( 62 I ,QI SB I ,TENO( SB I , 

1 TFRAMEt 30I.Q00T162 I.DQOOTI62 1.QWI 301, 

2 OKKI 62 >.D3St63>.NTAG(63),NTAGO.RtENP<62>.STP(62) 

COMMON ARGI6SI.ARG0 

INTEGER PLATE! 

REAL KL.KR.KfP.K.KUP.KDOWN 
DIMENSION ROI t*l,RBO< 16I.RL00I 1*1 
DIMENSION RLI 1*I,K(62),SI6*I,SS163I, 

I AI62 ) • AAI 62 ) ,QR( 62 I 
HEAT FLUX METER CALIBRATION CONSTANTS 
DATA K/ *0.07. St. 59, SO. *2, 

I SI. OB, 29.06, SI. 50, 

SI. 50. *1.0*, 20.65, 

29.^6, SI. 15, 2 7 . *7 , 

END LOSS CONSTANTS FOR TEST SECTION I PLATES 25- SB I I BTU/HR/F I 
.MB' 291. .256, .22*. .2*1 , . 2BS, .227, 

,U7t .227, .230, .225. .196, .296. 2**., 

HEAT FLUX METER CALIBRATION CONSTANTS • 59-62 I BTU/MR/5QF T l/MV 


• 1-2* I BTU/HR/SQF T l/MV 
32.79, IS. 97, 55.72, 
35.70, 3*. 36. 5*. *7, 
30.51. 27.39, 32.62, 
32.11, 3*. 37, 35.00, 


33. *6, 
25.16, 
33.01 , 
32.3*, 


33.05, 
33.10, 
33.91 , 
30. *1 . 


35.21 , 
27.15, 
31 .*7. 
32.70, 


S3. 6*, 
30.71, 
29.60, 
22 . 00 , 


33.53, 
30. * 2 , 

22 . 6 * , 

36.50, 


31.92, 
31.39, 
31 .*1 , 
32.72/ 


COMPUTE UNCERTAINTIES IN K< X I 


511 


IF! J.GT.1 I GO TO 132 

512 


DO 130 I s 1 .62 

513 


OKKI I 1*0.01 

51* 


IFII.LE.2*) OKKI 1 )*0 . 03*K( 1 1 

515 


IFtI.GE.39) OKKI 11*0. 03*K( 1 1 

516 

130 

CONTINUE 

517 

132 

CONTINUE 

510 


00 131 1*1,62 

519 


IC0UNT*I 

520 


IFIJ.EQ.10) CALL CHAN! ICOUNT ,K( I 1 1 

521 


XFIJ.CQ.M ) KII)tARG(I) 

522 

131 

r 

CONTINUE 

523 


DATA A / 2**0.15706, 0.2320, 13*0. 

52* 


IF! J.NE . 1 1 1 GO TO 220 

525 


00 120 1*1 .62 

526 


ICOUNT*! 

527 


AAI I )*A( X ) 

520 


CALL CHANI ICOUNT, All)) 

529 

120 

CONTINUE 

530 

220 

JFI J.NE . 1 2 ) GO TO 222 

531 


DO 121 X*1.62 

532 

121 

A( I )*AA( I ) 

533 

222 

CONTINUE 


C 

AXIAL CONDUCTION LOSS CONSTANTS 

53* 


DATA S / 2**0.7, 0.0*7, 13*1.0 

1 .5*5, .65*. .609, .696, .690 

2 .761, .697, .777, .730, .710 

3 ’.775, .693, .755, .677, .772 

535 


00 1*2 X>1 ,63 

536 


IFtJ.CQ.O) SI I 1-551 I ) 

537 

r 

IF I JSC .NE . 1 ) GO TO 1*2 


C 

r 

COMPUTE UNCERTAINTIES IN SID 

530 

Ir 

0SS(I)*0.10*S(X) 

539 


SSI I >*31 1) 


(BTU/HR/F) 


.710, 

.730, 

.770, 


.720, 

.7**, 

.776, 


.666, .62*, 

.72*. .7*6, 

. 00 / 
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6*0 

1*2 

COMT1NUC 

5*1 


I7lJ.Nt.7l 60 TO 2*1 

5*2 


N3*2)MS*62 

5** 


IFI JSC.EQ.1 1 CO TO 2*0 

5*5 


N3*t 1115*65 

5*7 

2*0 

00 1*0 I*NS.MS,2 

5*6 


I COUNT *X 

5** 


CALL CHANGE! I COUNT, JSC, 51 1 1,51 !♦ Ill 

650 


NTA61 I 1*0 

551 


XFII.Nt.1 1 60 TO 2** 

552 


NT AGO *1 

555 


CO TO 1*0 

55* 

2** 

NTACI 1-1 »*1 

555 

1*0 

CONTINUt 

556 

2*1 

CONTINUE 


C 

WATTMETER CIRCUIT RESISTANCES 

557 


DATA RO / 6.16, 6.05, 8.1*, 6.12, 6.05, 6.11, 

1 6 11, 6.25, 6.0*, 6.07, 6.17, 7. *6, 

556 


DATA RBO / 6.06, 7. *2, 6.06. 601, 7.92, 7.**, 

1 7. *7, 6.06, 7. *5, 7. *6, 6.06, 7.66, 

55* 


OATA RLOO/ 7.69, 7.6*. 7.65, 7.66, 7.6*. 7.6*. 

1 7.6*. 7.67. 7.65. 7.6/. 7.66. 7.67, 

560 


OATA RA.XA.RV.RD/ 0.060, 0.065, 7560.0, 0.1760/ 

561 


OATA KUP. ►.DOWN / 0.066, 0.179/ 

562 


00 150 I* 1,1* 

565 


ICOUNT*! 

56* 


IMJ.EQ.18) CALL CHANI ICOUNT.ROI I )) 

565 


inj.EQ.19) CALL CHANGE1 ICOUNT, 6.RDOI I I.ROID) 

566 


IMJ.IQ.20) CALL CHANGE 1 ICOUNT >*.R,.00( X ) ,RDO( I ) ) 

567 


inj.EQ.21) RIOOI I )*ARG( I ) 

566 

150 

CONTINUE 

56* 


IFIJ.EQ.2I) CALL CHANI 1,RA) 

570 


inj.EQ.22) RA*0 . 06* 

571 


IFIJ.LT. 2*1 GO TO 250 

572 


inj.EQ.2*) CALL CHANI 1 ,XA ) 

575 


inj.EQ.50) CALL CHANGE 1 1 .6.RV.XA) 

5 7* 


inj.EQ.31) CALL CHANGE! 1 ,*.RP.RV) 

575 


iriJ.EQ.32l CALL CHANGE! 1 .*. KUP. R6 ) 

576 


IFIJ.EQ.33I CALL CHANGE 1 1 .* .KOOWN ,KUP ) 

577 

250 

CONTINUE 

576 

* 

I1IJ.EQ.3*) KDOWN*0.17* 


v 

c 

THIS BLOCK CORRECTS INDICATED WATTMETER READING USING 


c 

WATTMETER CALIBRATION EQUATION 

57* 


00 10 1*25.36 

560 


JUC* 1-2* 

581 


ICOUNT*I 

562 


QP-QWI 1 1/75. 

565 


QCOR : QP*l 0 . 0726»QP-0 . 0*2 7»QP*QP-0 . 02*2 ) 

56* 


QCOR*0 . ***QW( I )*QC0R<>75. 0 

5S5 


IE (QUID . *T. 5.0) GO TO It 

566 


QCOR-QCOR* 1 0 . 00633/QP )*SIN< QP**0 . 55 ) 

567 

n 

r 

Xnj.EQ.17) CALL CHANI ICOUNT .QCOR ) 


t» 

c 

THIS BLOCK CORRECTS FOR WATTMETER INSERTION LOSSES 

566 


SUMRO-ROI JWC ) 

56* 


5UMRB0-PB0I JWC ) 

5*0 


2ROSQ-SUMRO*SUMRO*l XA*XA ) 

5*1 


ZPPOSQ=SUHPBO«SUMRBO 

5*2 


RLIJWC) * RLOOIJWC) - RB 

5*5 


ZVALSQ*! RV*RA*RLODI JWC ? 1*1 RV*RA«RLOOI JUC ) )«XA*XA 

5** 


Q< I )*QCOR*l ZROSQ/ZRBOSQ >•( ZVALSQ/RV/RV ) 
1 •( RLI JWC )/( RA«RLOO( JWC ) ) ) 

5*5 

10 

c 

CONTINUE 


V 

c 

THIS BLOCK CORRECTS POUER DELIVERED TO PLATES IN 


c 

TEST SECTION TOR CONDUCTION ANO RADIATION LOSSES 

5*6 


sr«t . 

i* 7 


EMIS*0.05 

5*6 


Xru.EQ.16) CALL CHANI 1.EMIS) 

5** 


TAR*I TINF**60. )/»C0. 

600 


00 10* 1*25.36 


6 . 0 *, 

7. *7/ 
7. *6, 
7.66/ 
7.65, 
7.66/ 
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*01 

601 

601 

006 

000 


000 

007 


000 

000 

010 

on 

on 

on 

016 

010 

010 

017 

010 

019 


020 

621 

622 

621 


026 

625 


626 


627 

620 

629 


610 

611 

612 

611 

636 

635 

636 

637 
630 

639 

660 

661 


662 


TOR*l T011 1*600. 1/100. 

OCONO * hi I !•( TFNEM 1 1-TMUnn 1 1 1 « 51 1 1*1 T01 1 l-TOI 1-1 1 1 

i ♦ si in in toii moan n 

17 II .10. 251 OCONO * OCONO ♦ KUP*I TOI 1 1 - TUP I 

17 II .(0. 101 OCONO • OCONO ♦ KOOt*MMTOtI) - TDOMN 1 

or ad* At I >*sr *cmxs* . 1 71 on top*tor*tor*tor-tar*t arm ar*tar i 

c 

C ENERGY BALANCE IS APPLIEO TO PUT! 

QOOTII1 • 10111*1.9129 - OCONO - QRA01/AIX) 

109 CONTINUE 

c 

C THIS BLOCK COMPUTES HEAT 7 LUX FROM DEVELOPMENT REGION PLATES 

S7 « 1 .0 

CHISM. 05 

I7IJ.EQ.16) CALL CHANM.EHIS) 

TAR*I TIN7 *660. 1/100. 

c 

00 200 I*PLATE1.26 
TOR*ITOII1*960. 1/100. 

QRI 1 1 ■ 5F*EHIS*.t 719*1 TOR*TOR*TOR*TOR-TAR*TAR*TAR*TAR> 
I7II.6T.il CO TO 195 

oooti 1 1 • Ki i i*hmi i i*i i .neo.-Tot 111/700, » - o»< 1 ) 

1 - I SI 1 1*1 TOI 1 1-TOO I ♦ SI 2 1*1 TOI 1 l-TOI 2 1 1 1/AI 1 1 
CO TO 200 

195 OOOTI I >«KI I (•HMI 1 1*11. *160. -TOI 1 1 1/700. I - OR. I I 

1 - I SUNN TOI 1 1-TOI 1-1 II ♦ SI 1*1 MM TOI I l-TO< in 1 1 l/AI 1 1 
200 CONTINUE 

c 

C THIS BLOCK COMPUTES HEAT 7 LUX FROM RCCOVERT REGION PLATES 

DO 201 1*19.62 
TOR*l TOI 1 1*660. 1/100. 

QRH) « SF»EMIS*.1719*I TOR*TOR*TOR*TOR-TAR*TAR*TAR*TAR I 
201 OOOTI 1 )*K( X )*HMt I 1*1 1 .*160. -TOI 1 1 1/700. I - QRIII 

1 - I SI I >•( TOI X l-TOI X-1 II ♦ SI 1*1 ••( TOI 1 1-TOI 1*111 l/AI X I 
C 

c 

RETURN 

END 


SUBROUTINE ENTHAL I FACT .D2INIT .D02INT .02LOC I 

c 

C COMPUTE ENTHALPY THICKNESS BASEO ON CONTROL 
C VOLUME FOR ENERGY ADDITION 

C 

INTEGER D2LOC 

COMMON/ BLK2 /UIN7 ,TINF .TAOIAB.RHOINF ,VISC ,PR ,CP,H 
COMMON/ BLK3 /TOO , TOI 63 I .HMI 62 1 ,Q( 36 1 , TENOI 16 1 , 

1 TFRAMEI 36 ) .OOOTI 62 I .DQOOTI 62 1 .QWI 16 I . 

2 OKKI 62 I ,DS5I 63 I .NTAGI 63 I , NT AGO .REINPC 62 I .STPI 62 I 

COMMON/ BLK6 /STI 62 I ,DSTI 62 I .REENI 62 1 .DREENI 62 1 ,021 62 1 ,0021 62 1 , 

1 PCOEFFt 62 I .XI 62 I 
DIMENSION DXI62I 

DATA OX/ 29*1 ,026.1 .656.13*1 .996.26*1 .026/ 

C ASSUMES 02 AT BEGINNING OF PLATE • 02LOC 
OOX « 0.005 

02ID2LOCI > D2INTT*ST(D2LOC )*DX(D2LOC 1/2. 

PEENI 12LOC )*FACT*D2( 02L0C !•! 1.0- P.'OEFFI 02LOC I 1*«0 .5 
DD2I02L0C) * D02INT 

OREENI 02L0C )*FACT*DD2( 02L0C >•( 1 .O-PCOEFFI 02LOC I >**0.5 
12 « 02LOC ♦ 1 
C 

DO 100 1*12,62 

IF II .LE. 29 OR. I ,6C . 391 GO TO 20 

0211 ) * 0211-1 1 ♦ 3TII-1 1*1 1 .0*. 0569*021 1-1 ))*OX( 1-1 1/2. 

1 ♦ STI !)•( 1.0*. 0569*021 1-1 1 1*0X1 1 1/2. 

2 *1 0.5*1 1 1 .0*. 0569*021 1-1 1 1/. 0569 1*1 ALOGI 1 .0* 

3.0569*021 1-1 ) 1 1/1 1 .O-PCOEFFI 1 1 1 1*1 PC^EFFI 1 1 

9 -PCOEFFI X-1 11-2. *111.0*. 0569*02* 1-1 ))/. 05691* 

51 ALOGI 1 .0*. 0569*021 1-1 1 1 1*1 TOI 1 1-TOI X-1 1 1/1 TOa l*TOI 1-1 )-2.*TINf 1 
IF II .EQ.251 021 1 1=021 X 1-1 1 021 1 >/0 . 05691-11 1.0 
1*. 0569*021 II)/. 0031 66*1 ALOGI 1 .0*0.0569*021* 1 1 1 1 1*. 0569 
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441 40 TO 14 

444 tO CONTINUE 

441 Otlll • Dttl«1l ♦ STI 1-1 1*0X1 1-ll/t. • OTI X 1*0X1 X l/t. 

1 ♦io.i*o*tl-t i/i t o-pcoeffi i nm pcoirn i » 

I -PCOEFFII-1 M-t. *01(1-1 t*l TOIH-TOI 1-1 M/ITOI 1 1*7011-1 l-t.*TINF » 
444 It II .td. SOI OH I t*Otl 1 1*1 1 Oil I )/•. 0144 >•111.0 

1*. 0164*0*1 in/. 001164*1 ALOBI 1 . 0*0 . 0044*011 X III 1 1*. 0444 
447 10 CONTINUt 

44* 00*1 1 1 ■ •* 1TI 00*1 i-n*** ♦ 1ST! 1-1 l*OOX/t. »••* 

1 ♦ IOXIl-1 *0ST!I-1>/t. >••* ♦ I OTI X l*OOX/t . >••* ♦ 
t 1 0X1 X l*DSTI X l/«. )••* 1 

440 PC EMI X I • FACT • 021 X 1*1 1 .O-PCOEFFI X I )**0.S 

*so on* ini 1 1 * * act • oozix) *t t .o-pcocfm i > i**o t 

Oft too CONTINUt 

*1* RETURN 

*»1 I NO 

c 

c 


*14 SUB* OUT IN* CNANIN.XX) 

C 

C SUBPOUT I Nf CHAN FIRST STOPIS TNI VALUC OF XX AS A 

C MMBIP OF TH* DUNN I APS A f ARC AND THIN INCRIAStS 

C XX BT ONt PtRCENT. 

C 

*11 CONMON AR6I *1 1 ■ ARGO 

*14 IFl N.NC . 0 > CO TO 10 

*17 A»G0«XX 

*10 CO TO *0 

*14 10 AR6INXXX 

**0 *0 XX*XX»XX*0.0t 

**t RETURN 

**i L NO 

c 

c 


6*1 SUBROUTINE CHANGE! N.H.XX, TY I 

C 

C SUBROUTINE RESTORES THE VALUE OF TY, THE PREVIOUSLY 

C PERTURBED INPUT, TO ITS ORIGINAL VALUE, THEN STORES 

C XX IN THE ARRAY ARC, THEN INCREASES XX BY ONE PERCENT. 


4*4 

COMMON ARC! 6 11. ARGO 

4*1 

XFIM.E4.1 .0R.N.CE.63) CO TO 10 

*4* 

YY«ARGIN«. ) 

647 

IFIN.LT. 4) GO TO 10 

*64 

ZFIN.EQ.V* YY*ARG0 

6*4 

1FIN.NC.0I YY*ARG( N) 

*70 10 

IF(N.NE.O) ARGINKXX 

671 

IFIN.EQ.OI ARG0»XX 

67* 

XX*XX*XX*0.01 

*71 

RETURN 

*74 

ENO 


C 

C 
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67* SUMOUTINC MAXHINI N.M, JCOUNT I 

c 

C SUBBOUT INI MAXMIN 7INDS TNI MAXIMUM AMO MINIMUM VAIUIS 

C IN THI COMMONIO ARRAY OVIB TNI BANOI Of I IRON N TO M 

C ANO MARKS TNI I -LOCATIONS 01 TNI MAXIMUM AND MINIMUM 


6 76 

C 

COMMON/ BIAS /SINCOI62.S6I, 

VMAX 

•VMIN.NA.NB 

677 


VMAX ! 0 . 0 



676 


VMINMO.O 



676 


DO 610 X«N.M 



660 


IM SINCOC I .JCOUNT 1 . 11 .VMAX 1 

60 

TO 620 

661 


VhAX*SCNCOI I. JCOUNT) 



662 


HI I 



66 1 


60 TO 610 



66* 

620 

171 SeNCOt I . JCOUNT 1 . IT . VMIN 1 

60 

TO 640 

60S 


CO TO 610 



666 

640 

VMIN*S(HCO< I. JCOUNT I 



667 


NB* I 



666 

610 

CONTINUI 



666 


I7( VMAX. IQ. 0.0) NA*NB 



660 


R I TURN 



661 


KM) 
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Appendix C 
The Reduced Date 

The following pages show lists snd plots of the reduced data of the 
IS cases shown In Tsble 3-1. Stanton number data and starting profile 
data are presented for each case. Additional velocity and temperature 
profile date are listed according to the following cases: 

Case # s-dlstance In cm of profile data, In addition to 

the a - - 35 cm starting profiles. 

100779 10.0 (0 - 13*) 

4 1 • 0 (0 - 52*) 

61.0 (0 - 78*) 

119.0 

070280 0.0 

10.03 (0 - 13*) 

25.5 (C - 32*) 

40.5 (0 - 51*) 

61.3 (0 - 77*) 

72.3 

88.2 

104.0 

119.1 


022680 

10.0 

(0 - 13*) 

(Profiles 

25.0 

(0 - 32*) 

are 

41.0 

(0 - 52’) 

030180) 

61.5 

(0 - 78°) 


89.0 



104.0 

119.0 

060480 9.9 (0 • 12.6°), velocity only 

25.0 (0 - 32*), velocity only 

051080 40.8 (6 - 52.5*), velocity only 
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STKt AM*ISf distance 1CM) 

rig. 5-1. Case 070280: Stanton number varaua streanwlse 

dlatanca. 



enthalpy thickness Reynolds number 

Fig. C-2. Case 070280: Stanton number veraua enthalpy 

thlcknesa Reynolds number. 
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STKEAMWISE DISTANCE (CM) 

Fig. G-3 . Case 022680: Stanton nuaber versus st Teamwise 

dlatanca. 



ENTHALPY THICKNESS REYNOLDS NUMBER 

Fig. G-4 . Case 022680: Stanton number versus enthalpy 

thickness Reynolds ruaber. 
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Fig. C-5. Cant* 060480: Stanton number versus itruwlit 
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Fig. C-6. Case 060480: Stantcn number versus enthalpy 
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f.NTHAU’l THICKNESS REYNOLDS NUMBER 

Fig. G-10. Case 011380: Stanton number versus enthalpy 

thlckneaa Reynolds number. 
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STREAM WISE DISTANCE (CM) 

Fig. G-13. Cab* 030280: Stanton number versus atrt-funvlse 

distance. 



ENTHALPY THICKNESS REYNOLDS NUMBER 
Fig. C-14. Case 030280: Stanton number versus enthalpy 

thickness Reynolds number. 
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STANTON N'-IHWJI 


tot It 


a 



STRCAMWI.E PISTAN* C (CM) 

Fig. C-17. Caaa 0174R0: Stanton nuabar varaua atraaawlaa 

dlatanca. 



ENTMAU’Y THICKNESS RFYNOIDS NUNfirR 

Fig. C-1B. Caaa 012*80: Stanton nuabar varaua anthalpy 

thicknraa Kaynolda nuabar. 
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STMtAMWISC NS »N< l (CM ) 


Fig. C-19. Caae 06 2 SCO 1 Stanton nuaber varaua itrruvUf 
dlatanca. 
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ENTHALPY THICKNESS REYNOLDS NUMBER 

fig. G-20. Caat 062880: Stanton nuabar varaua enthalpy 

thlckneaa Reynold* number . 
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srRLAMwisc distan. r it*> 


fig. 0-2’, Use 0*2580: Stanton number versus etreamwlse 

distance. 



CNTMALP V THICKNESS REYNOLDS NUMSt'N 

^lg. C-22 Use 0*2580: Stanton number versus enthalpy 

thickness Reynolds number. 
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STREAMNISC DISTANCE (CM) 

Fig . G-27 . Casa 051080: Stanton number varaua atrranwlar 

dlaiance. 



ENTHALPY THICKNESS REYNOLDS NUMBER 

Fig. C-28. Case 051080: Stanton nuaiber varaua enthalpy 

thlckneaa Reynold* number. 
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Appendix H 

Modifications to STANS for Cbnvax Curvature Predictions 


The following listing shows the modifications made to STANS [1] to 
account for convex curvature effects. This version of the program was 
used to make the predictions of Pigs. 4-2 through 4-17. New programming 
Is listed, bracketed with the preceding and following lines of the orig- 
inal STANS program (lines with Identification on right side). 


H-l 


c ruiNOtio 

C THIS PR OCR AH HAS BUM HOOITHO TOR CONVCX CURVATURC 

C USING THt 01 LI IS MOOtl WITH THl SIMON CHANCES 

C Sl( THl TMO THISIS TOR OtTAIlS 

C LATCST REVISION SIP TIMBER 11(9 

1 INTI lit H CION, TLUIO.SOURCICSI'SPACI.BOOTOR. OUTPUT. TTPBC HAIN0140 


S/AOO/RBOMI S4 I ,OHOI 54 I .RONOl SC I . ITKl ,CBCC 1 00 I ,CITT I 54 I .CLAC.CUNV. HAINOIVO 
4CU.PLIS4I,SKIM( IOOI.STANTI 1 00 * .XPLOTII 00 I . AltNTHM 00 I , IPLOT .XSTART 
S 0 1 Ml NS ION AHCP(»),AM17I»I MAIN0S00 

4 01 Ml NS I ON T PPC tOO ) . ATPP1 1 00 1 .31 PPI tool »CTPPC 1 00 ) MAINO SI 0 

5 CALL STARTS! 'OlNll*' .0. ) 

C MAIN0S11 


7 

N1I40O 

HA I NO 140 

A 

IPLOT*0 


0 

NPLOT»i 


10 

S K(RROR>0 

HAINOSSO 

II 

00 7 1*1.100 


II 

SKIN! I 1*0.0 


IS 

STAtlTl 1 1*0.0 


14 

XPIOTI II ’.(i.O 


IS 

REIMTHI I >*0.0 


14 

7 CONTINUI 


17 

NINO'MINOM 

HAINUS40 


1 9 ITIKIRROR.CT.0I60 TO 1000 HA.N0J90 

10 XSTARTiXII) 

It 00 A H*l,NXBC 

|| IT I CBCI M I .LT. -1.11-0 .ANO. CACCH-I) .CT. -I.OI'AI XSTART«XCM> 

IS A CONTINUI 

c - STIPS HAIN0400 


SO 


AUXM2*AUXl(M-1 )*C AUX2I M I'AUXK M- 1 M*CXU*XCM-1 M/C XI H )-X< H*l M 

HAIN0470 


c 

c 

CURVATURE MODIFICATION , INTCRP. HALL RAOIUS 


SI 


CU*CBC( M- 1 IMCBCIMl'CtCCH-l 1 t*JXU-XCH-1 ) )/l XI H )-X( H-1 ) 1 


SI 

c 

IT (KIX.CQ.1 1 AHI*AMIE 

HA1N04A0 


H-2 


HAIN3*90 


1*0 H ‘01 1 1 /Ot It 

Ml B***orn 

M? Rtf1«0f L2*UGU*RNG/VIS6 


71AIN3SOO 


y > s cau out 

So* SKI III XPLOT l*C72 

3*7 STANTI 1PLOT KSTMI 

5*6 I» < ST* « > .IT. I. 01**1 ST ANT I I PLOT 1*1.0 

SO* XPLOTl XPLOT 1*1 XU'XSTARTIOO.** 

ISO REiHTttl IPiOTMRIH 

SSI Ift CONTINUE 

ssi xmvAR.ar.nco to 1000 


MAIN* I 10 


MAIN*! SO 


ss* Auxtii*Auxnn-i iMAuxiim-Auxiin*i n»'xo-X(n*i n/ix<m-xin*i n haxn*2so 

C 

C CURVATURE MOO 1 7 IC AT ION, INURP. WALL RAOIUS 

s*o cw*cbcim*i iMCBC<m*cociri*n>*ixu*xih-i n/ixim-xin-i n 

s*i spm*AesiCbCinn 

c 

C xr XT XS DESIRED TO INTRODUCE Tilt UALL BOUNDARY COtiOl T XOM AS MAIN*2*0 


SSI 

sss 

so* 

SOS 

so* 

507 

SOS 

569 

590 

591 

so: 

SvS 

59* 

39S 

59* 

397 

390 

399 

*00 

*01 

*02 

*03 

* 0 * 

*0S 


60 TO IS n/.XN*S90 

1000 CONTINUE 

17 INPIOT ,(Q. II CALL PXCTR6I0.0) 

CALI SET3MG12S.0.0I 

CALL SETSMGI I* > 0.3 I 

CALL SUBJEGI -40. .0.0,1 30.0,0.00*1 

CALL SETSMGI 130,0.01 I 

CALL SETSMGI 131 ,0.011 

CALL SETSMGMSS.t.O) 

CALL CSETGI ’DUPLEX* * i 

CALL CRAPMGI ' • ’ ,1 , XPlOT , XPLOT , ST ANT , 0 , 'STRE AMUISE DISTANCE ICN)*' 

I,* STANTON NUMBER , C7/2*'.'*'l 
CALL LIHESGI 'OASM*' .XPLOT, XPlOT. SKINI 
CALL PXCTRGIO.OI 
CALL SETSMGI 2’ . 1 . 0 I 
CALL SETSMGI 2* , 1 . 0 I 
CALL CSETGI 'DUPLEX*' I 
CALL SUBJEGI 200. ,0.0005,5000.0,0.005) 

CALL SETSMGI 133,2.01 
CALL SETSMGI I 3*. 2.0) 

CALL GPAPIIGl .1 ,JPLOT,REENTH,STANT,0. 'ENTHALPY THICKNESS REYNOLD 
IS NUMOEP*’ , 'STANTON NUMBER* I 
NPLOr«t 

I7INXNO.lt. NUMRUNIGO TO S HAIN*600 

CALL CXITG 

«tqp MAIN** 10 
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1/AOO/RBOMI S* I.OWI SAl.HOflQl SO I.ITKC .COCI 1 00 > .CC r r ( *<• I.CLA6.CUPV, ST t PO I so 
OCW.PLttOl.SKINI 001. ST anti 100 I .XPLOTI 1 00 1 .RICNTHI 1 0* I • l PLOT .XST ART 
C STIPOIAO 


l/AOO/RBONI $S I . OnC 1 SO I .ROMO I SO ) . ITKC .CBCI lOOl.CCFPISol .CLAG.CURV. AUXOOISO 
OCU.PLISOI.SKIMI 1001. STAMIMOOI. XPLOTI iooi.runthi iooi.ipiot.xstart 

?220 DIMENSION OVI SO I , SMOt SO I .STAB! BO ) ■ VI LUNOI So I.TPll So I AUXOOISO 


mt TPUT»RMH*UTAU/VISW AUX002O0 

mi KTMBU* 0 

1210 IP I INTO. OT. II GO TO 10 AUX00250 


1 2 36 IP IM00C.CQ.2I MOUNT* I AUX00270 

c 

C CURVATURE MODIFICATION — GUESS INITIAL VALUES OP PARAMETERS 
1217 SLTML 

12 IS OXSSHRiBXX 

1210 ALAM«SLT-01SSMR 

I ZAO PLAN*ALAN 

1201 PZLAtUPLAN 

1202 PlLAN*P2LAH 

1201 I9LT*N 

1200 R£M5*BEN 

c 

I20S 1 RAVG-BI 1 I AUXOOISO 


C AUXOOOt 0 

1270 10 KTMRU*KTMflO’ 1 

1271 00 SO 1*2. NP| AUX00620 


COMPUTE NIX IMG LENGTH AUX01000 

C 

c CURVATURE MODIFICATION TO HAKE MIXING LENGTH 

1301 A LAMA * ALAN 

t 302 IF ‘XU • LT . XSTART ) ALANA*0.65*YL 

H-4 


nos 

1FIINT8.LT. 271 ALAMA*.6S«VL 


1306 

SO AL»AIMG«ALAF1A«I .1765 

AUX01010 

ISOS 

ALHAX'AL 


ISOb 

IF IKA3E.EQ.1 .ANO.VH.LT.AL/AKI AL>AK«rn 

AUXO) 020 


1316 

35 EMUT*RH0AV«AL«AL»ABS((UCX»1 l*U( 1 1 I/I TC !♦! I'YI 1 1 1 1'OVI I )*0V( 1 1 

AUXO 1060 

ISIS 

PLI II*AL/VL 


l Si 6 

SPI2MBXX 


* SI 7 

SP1SMSLT 


■ SI 8 

IFIK2.NE.2.0R.KA5E.EQ.21GO TO 36 

AUXOt 090 


1332 

IS3S 

1111 X»1 )-T( 1 1 1 • /( GC*TAUU I 

TPU I >*AB5( EMUI X )a(U( 1*1 I’UI I ) I/I Yl 1*1 )-Y( 1 1 ) )/CGC«TAUWI 
VELCRO! I l*IUl 1*1 >-U(III/IYII*1 »- VI 1 1 1 

AUXO 1260 

1336 

IF (NPM.EQ.OI GO TO 69 

AUX01250 


1357 

60 PREFU,I>»0.5«<PPIJ,I*1 I'PRU.IM 

AUX01560 

1358 

66 CONTINUE 


1359 

69 CONTINUE 


1 3b0 

YI20.00 


1361 

YCRIT*0.0 


1362 

IF1CH.GE.-1 .OE-6) GO TO 91 


1363 

5CRIT*0.25«U-CH*YL >••0.331 


1366 

DO 90 1*2, NPI 


1365 

VELGROI I U(UIII*UII-1 ll/«V(X»-VII-1 II 


1366 

STAB! I )=CWU( I l/VELGROI 1 1 


1367 

IFISTABIII.GT.-SCRITI GO TO 90 


1366 

Yl 2*1 -SCRIT-STABl 1-1 I >/( STAB! I l-STABI 1-1)1 


1269 

YCRIT*Y12«(Y(II-Y(I-1 ))*Y(I-1 1 


1370 

GO TO 91 


1371 

90 CONTINUE 


1372 

91 P3LAM=P2LAM 


1373 

P2'.Af1*PlAM 


1376 

PLAflsSLT-OISSHR 


1375 

NPSHRsNPI-1 


1376 

TAUMAX=0 . 0 


1377 

00 92 1*3 iNPSHR 


1376 

92 IFITPLIII.GT.TAUMAXI TAUHAX*TPL( X 1 


1379 

CNTR*0.6S 


1360 

BANO*0 .25 


1381 

93 KLUNK*0 


1382 

NFIT*0 


1 383 

5UM1 *0.0 


1366 

SUM2*0.0 


1385 

SW13*0.0 


13e6 

5UM6*0 . 0 


1367 

00 96 1*1 ,NPSHR 


1368 

UPPER* 1 CNTR+BANO )»TFUHAX 
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i sat 
isto 

i »ti 
i st; 

i s«s 

i St* 

i sts 

IStb 
1 St7 

i sta 
i sot 
IA00 
IhOI 
i*.o: 
iaos 

1 hOh 

iaos 

IsOb 
I *«o / 
iaoo 
uot 

IhIO 

Ml 
loll 
101 S 
101'. 
1010 
I o 1 6 
101 7 
loia 
Ult 

io;o 

10J1 

ia:: 

io;s 


io;o 


io:s 
io:b 
io:j 
io;a 
io;o 
to so 

loSI 

ios; 

1055 
10 So 

loss 

1056 

IAS7 
10S6 
1oS9 
1 0oO 
1 oo I 
1 OoJ 
looS 
1 ooo 
looS 
loot 

1oo7 


1 008 

1 oot 

1 oSO 


A t OWM • I CNTR - B AND )• T AUF1AX 
1MKLUUK.CT.0I CO TO to 
FISRC'NPI -1 

17 1 I PH NSRC I . GE . I 0 . t S* T AUMAX I 1 KLUNKM 
IF I TPLI FISRC I ,CT . UPF’f H . OR . T PL I N5HC I . IT . A l OUR I CO TO to 
NF X T*Nt 1 T ♦ 1 
ft'JHI >SUF11 *11 NSRC I 
SUfi;*5UMJ»TPLINSRC I 
SUMS SUMS* TPH NSRC loilNSRCI 
SUM-. • Si M* ♦ 1 PL I NSHl I • 1 Pit NSRC I 
tO CONTINUE 

ir i fit it .to. oi co to ta 

AFiUM SUMS M sum •SUMJ/NF1T 1 
0 1 NUFI * SUflo * I SUM." • SUMi /NF l T I 

iMoiMuM.it .o.oi u-iiuift. 06 ) i i.ti n.ut n.TPii i », xh.npshri 
ifioano.gt.o.Soi CO TO too 
IMOIMUM.CT.O.O.ANO NFIT.GT.il GO TO 100 
t6 FORMAT ( {X.IS.nX.F R.S.aX.F 7. {.AX, 7 7.01 

ta hanc*dano«i .:•> 

CO TO tS 

ioo ti*i sum -sum:*! a^jm/oinuhi i,nf it 

lFIIN1C.GT.1t. AND . IN1C.Lt . SOI UWXTEI6.t7l 11,1(11 >U( 1 I . 

ITPtlSKVtl (Mi ii.i’tm. t ■ i .mi’smr i 

1 7 FiRMATi:X,IS,AX,Ft.5,AX,F 7 . f , oX . F 7 . 0 . OX , F 7 . 0 . AX . F 7 . 0 I 
10{ CONTINUE 
SLT«1E 

UHI1K6.H1* 11 .1CHIT .SLT.CH 

IF ICN .01. -1.01 b .OH. 1CR1T .IT. 1.01-61 GO TO 107 
SLT»1CRIT 

IT lit .11. 1 CHIT I GO TO 107 

ACO't . 00 

SL1R* 1 .0-ACO«1 I U/1CRIT 1-1 .01 
IF I SLIP .LI. O.SSI SITH-O.SS 
SlT'SLTRMCRXT 

WHIU 16.1111 1!,1CRIT,SLT,CW 
111 FORMAT MX.'OILSL f ROM SMI AR SIRrSS PROF ILl »‘,F10.5 

I.IOX.'OLtSL FROM SCHIT*',Ft0.5,'0tLSL IN MIXING UNGTM •’ 
{.M9.S. 1 CW- *,710.61 
107 I7ISLT.GT.1 .{•ILl tLTM.Wt 
C 7 I NO NEW OUT AST AR UASLO ON SHLaR lAHR THICKNESS 
C F I NO mi w I0M VI UX Hi 

KSLT*0 

00 10! 1*{.NP1 

1MKSLT.CQ O.ANO.MII.GT.SLTI 15LT*1 
111 II I I.GT.5LTI KSLT»1 
10S CONI1NUE 

IF1K5LT .EG.O I GO TO 10S 
USLT«UGU 

0IS5HR-I 1 .-U1 { l/USLT 1*( Yttl/Sl 
OMOMT !Ul{ I/USITI-OISSHR 
1SLT1 *ISLT-1 
00 1 Oh X > 1 , 1SLT I 

OISSMR tOlSSHRM I 1 . UCX-I l/USLT )M 1 ,‘UI X 1/USLT 1 1 
|t»| VI s )-Y« I III . 

DMT I - 1 U( l l/USLT !•( I U( X I/U5LT II 
OMT C • I U( X - 1 I/U1 XSLTI I • 1 1 .-IUI 1*1 l/UI XSLTIII 
ioa ohomt 'Jfiomtmomtmofit:)a(1(Xi-ki-i ii/{. 

OXSSHRiOISSHHM 1 . Ul 1SL11 l/USLT )•« SLT-11 XSLTI 1 1/{ . 

DM I Mi (UlISlTII USLTI*( 1 .-U( XSLTI l/USLT »•! *LT-V* XSLTI ll/f. 

0M0MT« OMOMT ’OMINC 

RIFIS OtlOMI»U3'.T»tn,Oi 1 l/VISCOl 1 I 

105 XMKSLT.tQ.OI 0IS5MH OXX 
ALAH*SLT-OXSSHR 

IF I INTO. GT.O. ANO. INTO . LT . 2tt I UR1TEI6.1061 INTG.SLT .0X55MR , 

1 1 L.rIF XT, RIMS. ALMAX.K THRU 

106 FORMAT! IX, 1 INTO **,!<•, SX.'SLT » 1 , F6 . 5 , SX , 'OISSMR « , ,7a.S,5X, 
1 ' 01 Ltt* ' .77 . a, SX, • KF XT ■ ' .Ii.hX, 'RIMS*' , 7 7 . t , AX, 1 LMAX* ’ ,77.6, 
tfX.'KTMRU*' ,XSI 

IF ( At’S! ALAM PILAF1 1 . LE . I 0 . 0S«1L 1 1 ALAMM PL AM‘AL AM »/{ . 

ioa CONTINUE 

oo no x*:,npi 
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AUX01600 


1559 

1560 

1561 


18' 7 


2/CN/AXX,BXX,CXX.0XX,EXX,K) ,K2 .KJ.SPI 56 1 , AUX1 1 1 00 I . AUX2M 00 1 , YPMAX OUT00160 
3/AOO/RBOMI 56 I UMOI 56 I ,ROMO< 56 I , ITKI .CBCC I 00 I ,CEFF I 56 1 .CL/G.CURV, 

6CM , PLI 56 I , SKINI 1 00 1 , ST ANT1 1 00 I .XPLOTI 1 00 1 .REENTMI 1 00 I , IPIOT ,X5T ART 
C OUTOOI SO 


C 


PROPERLY, IT MUST BE THE LAST 
IE 1 INTG.GT.29.A1I0.IMTC.LT.S9I 
IPIIMTt.tT.S9) SPACE s 1 0 
IF! IMTG.GT.2 9.AII0. INTG.LT. 19 I 
in lMTt.tT.S9) KSPACE 1 1 0 
IF) IMTG.GT.59.AMO.INTG.LT.701 
IP I IMTG.GT . 70 ) SPACE«10 
IF( IMTG.GT 59.AMO.INTG.LT.70I 
IFI IMTG.GT. 701 K3PACEM0 
IE I KIM, ME . 1 1G0 TO 600 


EQUATION SOLVED. 
SPACE*) 

KSPACE* I 

SPACE*) 

KSPACE*) 


OUT00260 


OUT00250 


WRITEI6,28h»MINTG.XU.UGU.CAY,FAH,R£M,CF 2,M,REH,STA,F( 1,1 i.CPL.AME OUT00660 
IPLOT*IPLOTM 

IE1K1 .GT.10 IWRITEI 6,286 M 5P1 1 1,1*1 .51.G.BTA QUT00650 


2/CM/AXX ,BXX ,CXX ,OXX , EXX ,K 1 ,K2 ,K3 ,SP( 56 1 , AUX1 (1 00 1 , AUX21 1 00 I .YPMAX 1NPU0160 
S/AOO/ROOMI 56 1 ,01101 5s I ,ROMO( 56 1.ITKE ,CBCt 100 1 .CEEF1 56 1 .CLAG.CURV , 

6CW, PLI 56 1 ,SKIN( I 00 ) , ST ANTI 100),XPLOT( lOOI.REENTMl 1 00 ) , IPLOT.XSTART 

IMPU01 50 


C I DECIMAL NUMBERS, IN THE FORM OF A TABLE.) INPU2160 

C * 

C CURVATURE MODIFICATION, RE AO HALL CURVATURE 

90 READ! 5,580 1 X( M 1 ,RH( M I , AUX1 ( M 1 , AUX21 M) ,CBC( M I INPU2150 

C * 

r XNPU2160 
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INPUISOO 


1054 


1000 

1001 

1002 


• 

C COECIMAL NUMQERS I 

REACH 5.560 I C LAG i CUR V 

c 

C CL AG IS THE CONSTANT IN THE EFFECTIVE CURVATURE LAG EQUATION. 

C RE COMMENDED VALUE IS J.J, CUHV IS THE MIXING LENGTH MODIFICATION 
C CONSTANT . RECOMMENDED VALUE IS 0. 

C 

C (DECIMAL MUM0ERS ) INPU1510 

REAO( 5 . S0O I GC.CJ.AXX.BXX.CXX.OXX.IXX INPU-.010 

CXX*CURV 
AXX*CLAG 

C INPUsOCO 
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